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

    Flux Footprint Climatology Estimated by Three Analytical Models over a Subtropical Coniferous Plantation in Southeast China

    2015-01-05 02:02:06ZHANGHui張慧andWENXuefa溫學發(fā)
    Journal of Meteorological Research 2015年4期
    關(guān)鍵詞:張慧

    ZHANG Hui(張慧)and WEN Xuefa(溫學發(fā))

    1 Jinzhou Ecology and Agriculture Meteorological Center,Jinzhou 121001

    2 Key Laboratory of Ecosystem Network Observation and Modeling,Institute of Geographic Sciences and Natural Resources Research,Chinese Academy of Sciences,Beijing 100101

    Flux Footprint Climatology Estimated by Three Analytical Models over a Subtropical Coniferous Plantation in Southeast China

    ZHANG Hui1(張慧)and WEN Xuefa2?(溫學發(fā))

    1 Jinzhou Ecology and Agriculture Meteorological Center,Jinzhou 121001

    2 Key Laboratory of Ecosystem Network Observation and Modeling,Institute of Geographic Sciences and Natural Resources Research,Chinese Academy of Sciences,Beijing 100101

    Spatial heterogeneity poses a major challenge for the appropriate interpretation of eddy covariance data. The quantification offootprint climatology is fundamentalto improving our understanding ofcarbon budgets, assessing the quality of eddy covariance data,and upscaling the representativeness of a tower flux to regional or global scales.In this study,we elucidated the seasonal variation of flux footprint climatologies and the major factors that influence them using the analytical FSAM(Flux Source Area Model),KM(Kormann and Meixner,2001),and H(Hsieh et al.,2000)models based on eddy covariance measurements at two and three times the canopy height at the Qianyanzhou site of ChinaFLUX in 2003.The differences in footprints among the three models resulted from different underlying theories used to construct the models. A comparison demonstrated that atmospheric stability was the main factor leading to differences among the three models.In neutral and stable conditions,the KM and FSAM values agreed with each other,but they were both lower than the H values.In unstable conditions,the agreement among the three models for rough surfaces was better than that for smooth surfaces,and the models showed greater agreement for a low measurement height than for a high measurement height.The seasonal flux footprint climatologies were asymmetrically distributed around the tower and corresponded well to the prevailing wind direction, which was north-northwest in winter and south-southeast in summer.The average sizes of the 90%flux footprint climatologies were 0.36-0.74 and 1.5-3.2 km2at altitudes of two and three times the canopy height,respectively.The average sizes were ranked by season as follows:spring>summer>winter>autumn.The footprint climatology depended more on atmospheric stability on daily scale than on seasonal scale,and it increased with the increasing standard deviation of the lateral wind fluctuations.

    eddy covariance,flux footprint,flux footprint climatology,model comparison

    1.Introduction

    Spatial heterogeneity poses a major challenge for the appropriate interpretation ofeddy covariance(EC) data(Baldocchi,2008;Kljun,2010a;Biermann et al., 2011;Leclerc et al.,2014).Flux footprint climatology is a practical method to illustrate the portion of the sampled landscape that contributes the most to EC vertical turbulent flux at a given point over long-term periods(Amiro,1998;Kljun,2010a;Cai et al.,2011; Aubinet et al.,2012;Zhang et al.,2012).This method is used to understand carbon budgets,assess the quality of EC data,and upscale the representativeness of a tower flux to regional or global scales(Rebmann et al.,2005;G¨ockede et al.,2008;Chen et al.,2012).

    Flux footprint climatology,which is the long-term pattern of the flux footprint(Chen et al.,2009),provides a quantitative estimate of changes in the spatial representativeness of EC over time(Allen et al.,2011; Chen et al.,2011,2012).It is important to studywhether the factors affecting flux footprint climatology are the same as those affecting the flux footprint itself, such as wind direction,atmospheric stability,surface roughness length,and measurement height(Aubinet et al.,2012).It is not clear whether atmospheric stability,a parameter that fluctuates daily,strongly affects flux footprint climatology on monthly or seasonal scales.In addition,quantitative analysis to determine how much the size of the flux footprint climatology increases with measurement height is also necessary.Because trees grow taller over time,EC measurements are performed progressively closer to the canopy,which can impact the footprint climatology of flux measurements.

    Analytical models are often used for estimating flux footprints due to their ease of coding compared with the Lagrangian stochastic models(Leclerc and Thurtell,1990;Kljun,2010b),large eddy simulations(Cai et al.,2010),and closure models(Sogachev and Lloyd,2004).Currently,analyticalmodels such as the FSAM(Flux Source Area Model)developed by Schmid(1994),the KM model by Kormann and Meixner(2001),and the H model by Hsieh et al. (2000),are widely adopted.However,the strength, weakness,and applicability of each model have not been fully demonstrated.How the three models differ with variations in measurement height,surface roughness length,and atmospheric stability needs to be assessed.

    The objective of this study is to elucidate the seasonal variation of flux footprint climatologies and the major factors that influence them by using the analytical FSAM,KM,and H models based on eddy covariance measurements at two and three times the canopy height at the Qianyanzhou site of ChinaFLUX.ChinaFLUX is an observation and research network that uses eddy covariance and chamber methods to measure the exchanges of carbon dioxide,water vapor, and energy between the terrestrial ecosystem and the atmosphere in China(Yu et al.,2006).A comparison of the three models will be performed to evaluate their differences with changes in measurement height, surface roughness length,and atmospheric stability. Then,the three analytical footprint models are used to reveal the seasonal variations of the flux footprint climatologies.Finally,the major factors influencing the sizes of the flux footprint climatologies over the Qianyanzhou site are investigated.

    The paper is organized as follows.Materials and methods are described in Section 2.Comparison of the three flux footprint models and the simulated seasonal variations of the flux footprint climatologies are presented in Section 3.Effects of atmospheric stability,wind fluctuations,and measurement height on the flux footprint climatology are discussed in Section 4. A summary and conclusions are given in Section 5.

    2.Materials and methods

    2.1 Site description

    The Qianyanzhou(QYZ)flux site(26°44?52??N, 115°03?47??E)has been described in the previous studies such as Wen et al.(2006,2010)and Chen et al.(2010).This site is featured with a coniferous plantation forest in the subtropical continental monsoon region of China.The forest cover reaches 90% in the 1-km2area surrounding the tower and 70%in the 100-km2area surrounding the tower.Gently undulating terrain surrounds the site,with slopes of between 2.88°and 13.58°.Pinus elliottii,Pinus massoniana,and Cunninghamia lanceolata are the dominant tree species,with an average height of approximately 12 m in 2003.The prevailing wind directions are north-northwest in winter and south-southeast in summer.

    Eddy flux was measured at about two(23.6 m) and three(39.6 m)times the canopy height with two open-path eddy covariance systems consisting ofopenpath CO2/H2O gas analyzers(model LI-7500,Licor Inc.,Lincoln,Nebraska)and three-dimensional sonic anemometer/thermometers(model CSAT3,Campbell Scientific Inc.,Logan,Utah).The signals of these instruments were recorded at 10 Hz by CR5000 data loggers(model CR5000,Campbell Scientific Inc.,Logan, Utah)and block-averaged over 30 min for analysis and archiving.

    The eddy dataset was subject to a series of data quality control steps(Wen et al.,2010;Tang et al.,2014a,b).First,spurious data were detected in the datasets related to rainfall,water condensation,system failure,and insufficient turbulent mixing during the night.Second,planar fit rotation was applied, at monthly data intervals,to remove the effect of instrument tilt or irregularity on the airflow(Wilczak et al.,2001).Third,a Webb,Pearman&Leuning (WPL)correction was applied for removing the effect of fluctuation in air density on the fluxes of CO2and water vapor(Webb et al.,1980;Leuning,2005).Finally,to avoid possible underestimation of flux during stable conditions at night,the values of net ecosystem productivity(NEP)and evapotranspiration(ET)were excluded when the value of friction velocity,u?,was less than 0.17 m s?1.

    2.2 Flux footprint models

    Generally,the total flux footprint,f(x,y,zm) (m?2),is defined as the product of the crosswindintegrated footprint function,fy(x,zm)(m?1),and the crosswind dispersion function,Dy(x,y)(m?1);i.e.,

    where the field of f(x,y,zm)projected onto an x-y plane is called a source area.In practice,however,it is often desirable to obtain an estimate of the source area that is responsible for a given contribution,P, to the value of the measurement.Therefore,for the contribution P,the source area isΩP(Schmid,1994). In this study,to evaluate the agreement of the FSAM, KM,and H analytical models,the results for the 90% source area(P=90%)estimated by these three models are compared.

    A brief review of fyin the FSAM,KM,and H models is provided in the following section.

    2.2.1 The FSAM model

    The fyof the FSAM model is based on the twodimensional advection-diffusion equation and considers realistic velocity profiles and their atmospheric stability(Schmid,1994).The final function is

    2.2.2 The KM model

    The fyofthe KMmodeltakes into account atmospheric stability and uses the power law profile to calculate the verticalprofile of the eddy diffusivity,K(z), and vertical profile of the wind speed,u(z)(Kormann and Meixner,2001).The final function is

    2.2.3 The H model

    The fyof the H model uses two similarity constants,D and Q,to reflect atmospheric stability (Hsieh et al.,2000),according to the Gash(1986)footprint model for neutral atmospheric conditions.The final function is

    where L is Monin-Obukhov length,K=0.4 is von Karman’s constant,e=2.718,and zuis a length scale that combines zmand z0,i.e.,

    The two similarity constants,D and Q,are

    2.3 Parameterization

    Table 1 shows the initial parameters of the flux footprint for comparison of the three flux footprintmodels.The two measurement heights(zm)are 23.6 and 39.6 m.The measurement height should ensure that the EC is high enough to view the whole study area of interest.The surface distance over which wind has blown is termed fetch,and a rule of thumb,i.e.,the ratio of measurement height and fetch equals 1:100,is generally used.The two surface roughness lengths(z0)are 0.06 and 0.58 m.The Monin-Obukhov length,L,is-158 m for unstable conditions, 165 m for stable conditions,and 14048 m for neutral conditions.Near-neutral conditions are specified for |(zm-d)/L|<0.01 according to Hsieh et al.(2000). Therefore,(zm-d)/L<-0.01 is considered for unstable conditions and(zm-d)/L>0.01 is considered for stable conditions.The three friction velocities(u?)are 0.48,0.40,and 0.26 m s?1.The three standard deviations of lateral wind fluctuations(σv)are 0.99,0.75, and 0.58.The zero-plane displacement(d)is 8.4 m.

    Table 1.Initial parameters for the flux footprints:measurement height zm(m),surface roughness length z0(m), friction velocity u?(m s?1),standard deviation of lateral wind fluctuationsσv(m s?1),and Monin-Obukhov length L(m)

    The source area footprint results from the three models are compared.The variables a(near the end of the source area),i(the far end of the source area), j(the maximum lateralhalf-width of the source area), and Ar(the source area)are the characteristic dimensions of the source area.For the three models,the sensitivity agreement of the source area to(a)measurement height,(b)surface roughness,and(c)atmospheric condition can be assessed through comparing the four characteristic dimensions.The FSAM has been widely adopted but only applied during neutral and moderate atmospheric stability and with a limited range ofcrosswind turbulence intensity(Schmid,2002; Vesala et al.,2008).Therefore,the FSAM results are used as a reference and compared with the other two models.

    2.4 Flux footprint climatology

    To obtain the flux footprint climatology,the three analytical flux footprint models were run at a 30-min time step,and the data were accumulated to yield seasonal values for each pixel(x,y,zm-z0)separately. The accumulated values of each pixel were normalized by the cumulative seasonal values of the area of interestΩPto yield the flux footprint climatology(Chen et al.,2009),γ(x,y,zm-z0),as

    where i is the time step(i.e.,30 min),N is the total number of 30-min periods within a season,andΩPis the area of interest which gives a contribution,P,to the measurement.In implementation of the models, values forγ(x,y,zm-z0)were sorted in a descending order and accumulated from largest to smallest until a given fraction,P,was achieved.Finally,P-level profiles were produced.The calculated flux footprint climatology provides a map of the area around the tower that has contributed to the EC-measured flux.

    3.Results

    3.1 Comparison of the three flux footprint models

    3.1.1 Differences among the three models

    Figure 1 illustrates the 90%source areas predictedby the FSAM,KM,and H models at different measurement heights,surface roughness values,and atmospheric conditions.Table 2 shows the values of a,i,j, and Ar of the 90%source area as predicted by the three models at different measurement heights,surface roughness values,and atmospheric conditions.Figure 1 and Table 2 indicate that atmospheric stability is the main factor that has led to differences among the three models.

    In neutral and stable conditions,the KM results agree well with the FSAM results.The discrepancies in Ar between KM and FSAM are less than 3%.However,the H results are larger than the FSAM results. In neutralconditions,the discrepancies in Ar between H and FSAM are 4%-11%;whereas in stable conditions,the discrepancies increase to 18%-26%.When the measurement height is 23.6 m with a smooth surface,the discrepancies in Ar between H and FSAM are lowest in neutral conditions and highest in stable conditions.

    In unstable conditions,the agreement between the three models is better for a rough surface than for a smooth surface(Fig.1 and Table 2).Additionally,the agreement between the three models for the lower measurement height is better than for the upper measurement height.For smooth surfaces,the Ar magnitude follows the order:KM>FSAM>H;however,for rough surfaces,the order changes to KM>H>FSAM.When the measurement height increases to 39.6 m,the magnitude of results across the three models does not change.Thus,among the models, the differences due to roughness length are larger than those due to measurement height.

    3.1.2 Reasons for the differences among the three models

    Fig.1.The 90%source areas predicted by the FSAM,KM,and H models at different measurement heights,surface roughness values,and atmospheric conditions.

    The footprint differences among the three models originate from the different underlying theories used to construct each model.The FSAM and KM models are theoretical models,whereas the H model is an empiricalmodel.Thus,in Table 3,the FSAMand KMmodels have the same parameters,but these are not shared by the H model.The FSAM and KM models are based on a two-dimensional advection-diffusion equation,or K-theory(Schmid,2002;Foken,2008; Vesala et al.,2008),whereas the H model uses regression analysis to obtain three sets of empirical constants for D and Q in unstable,neutral,and stable conditions(Eq.(6))and forms a footprint model for thermally-stratified atmospheric flows(Hsieh et al., 2000).Therefore,in the present study,when atmospheric stability changes,changes in the H model are larger than those in the other two models.

    The discrepancy between FSAM and KM is caused by the same parameter being calculated by different methods(Horst and Weil,1992;Schmid,1994; Kormann and Meixner,2001).Table 3 shows the pa-rameterization methods of the FSAM and KM models.Monin-Obukhov similarity profiles are used in the FSAM model,and power law profiles are used in the KM model.Figure 2 illustrates the differences in the mean plume height for dispersion,effective speed of plume advection(U),gradient ofwith xand shape factor(r)between the FSAM and KM models in neutral conditions when measurement height is 23.6 m and roughness length is 0.58 m.If,U,and r were the same for FSAM and KM,then the results of FSAM and KM would be equal(i.e.,Eq.(2)is equivalent to Eq.(3)),while r is constant in the KM model but variable in FSAM.Differences in the other three parameters between the two models increase with the increasing upwind distance.

    Table 2.Values of the variables a,i,j,and Ar of the 90%source area predicted by the FSAM,KM,and H models at different measurement heights,surface roughness values,and atmospheric conditions

    Table 3.The parameterization methods for the mean plume height for dispersioneffective speed of plume advection(U),gradient ofwith xand shape factor(r)adopted in the FSAM and KM models

    Table 3.The parameterization methods for the mean plume height for dispersioneffective speed of plume advection(U),gradient ofwith xand shape factor(r)adopted in the FSAM and KM models

    a*:ˉz must be calculated numerically(Horst and Weil,1992).b*:K=0.4 is von Karman’s constant;p=1.55 is also a constant, andψ(pˉz/L)andφc(pˉz/L)are the diabatic integrations of the wind profile and stability function of heat.The specific calculations ofψ(pˉz/L)andφc(pˉz/L)were described by Horst and Weil(1992).Note:a and b are the constants of power laws for the eddy diffusivity and the vertical profile of the horizontal wind velocity,respectively(Kormann and Meixner,2001).The H model,which is empirical,does not use these parameters.

    KM H FSAM ?1 ˉz(x)a?ˉz(x)=∫∞0 ?r2 b z Cy(x,z)d z C y(x,z)d zˉz(x)=B ∫∞0 r x1r -U(x)U(x)=u?? ln 0.562ˉz ?? ?0.562ˉz U(x)=Γ(μ)?m?r2b a -ψ z0 k L Γ(1/r)araxmr-dˉz b(ˉz)1?rr r is variable r is constantb?d z dˉz d x d x=K2[ln(pˉz/z0)-ψ(pˉz/L)]φc(pˉz/L)d x=rBra

    Because the KMresults agree wellwith the FSAM results in stable conditions(Fig.1),the KM and FSAM models can be used equally at night(i.e.,under stable atmospheric conditions).It is considered that the H model should not be used immediately before and after dusk(dawn)when atmospheric stability changes because this causes the H results to change more significantly than the other two models(Fig.1). The discrepancies among the three models could be neglected if the dimension of the study area is sufficiently large.

    Based on previous studies,the strength,weakness,and applicability of the three models are reviewed in Table 4.The full version of FSAM is available on the website http://www.indiana.edu/~climate /SAM/SAM?FSAM.html,but some parameters,such asˉz,must be calculated numerically(Horst and Weil, 1992).Therefore,although the FSAM model is expressed by analytical formulae,it can only be evaluated numerically(Schmid,2002).The KM and H models provide analytical solutions because of their mathematical simplicity(Foken,2008).The KM model generally overestimates wind velocity near the ground,especially for unstable conditions and large roughness values(Kormann and Meixner,2001); while,in the present study,a change in atmospheric stability caused the H results to change in an inconsistent manner(Fig.1).The FSAM and KM models are restricted to surface layer scaling conditions due to the Gaussian distribution in the crosswind direction;whereas the H model is applicable to short-term changes in atmospheric stability.

    Fig.2.Differences of(a)mean plume height for dispersion,(b)effective speed of plume advection(U),(c)gradient with an upwind distance ofand(d)shape factor(r)between the FSAM and KM models.

    Table 4.Strength,weakness,and applicability of the FSAM,KM,and H models summarized from previous studies and from this study

    The conclusions presented in Table 4 are also validated by this study.For example,when atmospheric conditions changed from unstable to stable,the area given by the H model increased by more than double (much more than for the FSAM or KM models)(Fig. 1).When L<0,the size given by the KM model was the largest of the three models(Fig.1).This could be explained by the weakness of the KM model. For example,the KM modeloverestimated wind velocity near the ground,especially for unstable conditions. Because the wind velocity(u(z))was overestimated, the effective speed of plume advection(U)could also have been overestimated(U=u(z)/zm)(Kormann and Meixner,2001)and the EC sensor could access the further source according to the inverted plume assumption(Schmid,2002).In this scenario,the area of the footprint would be overestimated.

    3.2 Seasonal variations of the flux footprint climatology

    Figure 3 illustrates the seasonal variations of the flux footprint climatology as predicted by the three models for receptor locations at two and three times the canopy height.For clarity,the corresponding wind directions are also shown.There are similar seasonal patterns in the flux footprint climatologies assessed by the three models.

    The seasonal flux footprint climatologies are distributed asymmetrically around the tower and correspond well to the prevailing wind direction.In 2003, the 90%flux footprint climatologies of the three models at 23.6 m were 27%-29%from the north-northwest part of the receptor and 26%-28%from the southsoutheast part of the receptor.In summer,the flux footprint climatology in the south-southeast region increased to 51%.In autumn,the main distribution of the flux footprint climatology returned to the northnorthwest where 38%-52%contributed to the receptor.In winter,a 33%-50%contribution also originated from the north-northwest.The distributions of the flux footprint climatologies at 39.6 m were similar to those at 23.6 m.From spring to winter,the contribution of the north-northwest part at 39.6 m was 5% less than that at 23.6 m,but the contribution of the south-southeast part did not significantly change.

    The average sizes of the 90%flux footprint climatologies of the three models in 2003 for the receptor at 23.6 m were 0.61-0.74 km2in spring,0.45-0.67 km2in summer,0.38-0.42 km2in autumn,and 0.36-0.69 km2in winter.The seasonalaverage sizes followed the order:spring>summer>winter>autumn,whereas the seasonal changes in the sizes followed the order:winter>summer>autumn>spring.The sizes for the receptor at 39.6 m were on average 4.6 times larger than those at 23.6 m.The sizes were 2.69-3.19 km2in spring,2.07-3.28 km2in summer,1.74-2.36 km2in autumn,and 1.54-2.90 km2in winter.The difference between the two measurement heights did not change significantly with season.

    Fig.3.Seasonal variations of the flux footprint climatologies(left panels)simulated by FSAM(red),KM(black), and H(blue)models at 23-(solid)and 39-m(dotted)measurement heights,and the corresponding wind rose diagrams (right panels),for(a1,a2)spring,(b1,b2)summer,(c1,c2)autumn,and(d1,d2)winter over the subtropical coniferous plantation in Southeast China.The tower position was located at(0,0).

    4.Discussion

    4.1 Effect of atmospheric stability on flux footprint climatology

    Footprint climatologies depend on atmospheric stability on a daily scale.For example,Chen et al. (2008)revealed that the flux footprint area in the daytime was smaller than that at nighttime because more unstable conditions occur during daytime and more stable conditions occur at night.During stable conditions,the measurements can be affected by turbulence characteristics of the vertical wind velocity from further upwind(G¨ockede et al.,2006;Foken,2008).

    However,footprint climatologies less clearly depend on atmospheric stability on a seasonalscale.Figure 4a shows the linear regression(R2=0.38)between the size of flux climatology and the frequency of a stable atmosphere.Atmospheric stability varied for a short time;therefore,its influence on footprint climatologies was weaker for continuous observations over long periods.

    4.2 Effect of the standard deviation of lateral wind fluctuations on flux footprint climatology

    Figure 4b shows the linear regression(R2=0.73) between the size of flux climatology and the standard deviation of lateral wind fluctuations.The size of the flux footprint climatology increased with the increasing standard deviation of the lateralwind fluctuations. In Eq.(1),the crosswind dispersion,Dy,is generally assumed to be Gaussian(Horst and Weil,1992),

    Fig.4.Linear regression between the size of flux climatology and the(a)frequency of stable conditions and(b)standard deviation of the lateral wind fluctuations.

    4.3 Effect of measurement height on flux footprint climatology

    The results in this study showed that the size of the flux footprint climatology could be increased by three-fold when the measurement height increased to twice the canopy height(Fig.3).The majority of current footprint models rely,either implicitly or explicitly,on the inverted plume assumption(Schmid, 2002).The mean plume height(ˉz)and the effective speed of plume advection(U)which are characteristic parameters of plumes,depend on the measurement height according to the equations in Table 3(Horst and Weil,1992).Figure 5 illustrates the effect of measurement height onˉz and U.Variablesˉz and U were the results of the FSAM model for an unstable condition.When the measurement height increased,the mean plume height,ˉz,also increased and the effective speed of the plume advection,U,decreased;therefore, turbulence may have been received from a more distant source area(Aubinet et al.,2012).Therefore,the EC sensors should be high enough to access the environment through,and above,the main plant canopy of interest(Aubinet et al.,2012).

    5.Conclusions

    In this study,we elucidated the seasonal variations of flux footprint climatologies and the major factors that influence them using the analytical FSAM, KM,and H models based on eddy covariance measurements at two and three times the canopy height at the Qianyanzhou site of ChinaFLUX.

    A comparison of the three analytical models demonstrated that atmospheric stability was the main factor that led to differences among the three models on a daily scale.In neutral and stable conditions, the KM results agreed well with the FSAM results, whereas the H results were larger than the KM and FSAM results.In unstable conditions,the agreement among the three models for rough surfaces was better than that for smooth surfaces.Further,the lower measurement height(23.6 m)produced better agreement than the upper measurement height(39.6 m).The differences due to roughness length were larger than those due to measurement height.The differences in footprints among the three models were caused by the different underlying theories used to construct each model.The FSAM and KM models are theoreticalmodels,whereas the H model is empirical.

    Fig.5.The effect of measurement height on(a)mean plume height(m)and(b)effective speed of plume advection(m s?1).

    The FSAM model is available online,but some parameters,such asˉz,must be calculated numerically. The KM and H models provide analytical solutions because of their mathematical simplicity.The KM model generally overestimates wind velocity near the ground,especially for unstable conditions and large roughness values.When atmospheric stability varies, the H results change in an inconsistent manner.Because the FSAM and KM models both assume Gaussian distributions in the crosswind direction,and all the parameters ofthe two models are expressed for surface layer,these models are restricted to surface layer scaling conditions.The present study reveals that the H model is applicable to low-frequency changes in atmospheric stability.

    This study found that the seasonal flux footprint climatologies were distributed asymmetrically around the observation tower and corresponded well to the prevailing wind direction.The predominant distributions were north-northwest in winter and southsoutheast in summer.The average sizes of the 90% flux footprint climatologies were 0.36-0.74 and 1.5-3.2 km2at two and three times the canopy height,respectively.The seasonal average sizes followed the order: spring>summer>winter>autumn.

    The footprint climatologies clearly depended on the atmospheric stability on daily scale,but less clearly on seasonal scale.This was because atmospheric stability is only variable for a short time, and its influence on footprint climatologies is weaker for continuous observations over long periods.The size of the flux footprint climatology increased as the standard deviation of the lateral wind fluctuations increased and,therefore,crosswind dispersion increased.The flux footprint climatology increased three-fold when the measurement height increased to twice the canopy height.

    Acknowledgments.We thank Dr.Schmid for providing the source code of the FSAM model online.

    REFERENCES

    Allen,R.G.,L.S.Pereira,T.A.Howell,et al.,2011: Evapotranspiration information reporting.I:Factors governing measurement accuracy.Agr.Water Manage.,98,899-920.

    Amiro,B.D.,1998:Footprint climatologies for evapotranspiration in a boreal catchment.Agr.Forest. Meteor.,90,195-201.

    Aubinet,M.,T.Vesala,and D.Papale,2012:Eddy Covariance:A Practical Guide to Measurement and Data Analysis.Springer,New York,21-58 pp.

    Baldocchi,D.D.,2008:“Breathing”of the terrestrial biosphere:Lessons learned from a global network of carbon dioxide flux measurement systems.Aust.J. Bot.,56,1-26.

    Biermann,T.,W.Babel,E.Thiem,et al.,2011:Energy fluxes above Nam Co Lake and the surrounding grassland—The NamCo 2009 experiment.7th Sino-German Workshop on Tibetan Plateau Research,Hamburg,Germany,3-6 March,German TiP project and Institute of Tibetan Plateau Research.

    Cai,X.H.,J.Y.Chen,and R.L.Desjardins,2010: Flux footprints in the convective boundary layer: Large-eddy simulation and lagrangian stochastic modelling.Bound.-Layer Meteor.,137,31-47.

    Cai,X.H.,M.J.Zhu,S.M.Liu,et al.,2011:Flux footprint analysis and application for the large aperture scintillometer.Adv.Earth.Sci.,25,1166-1174.

    Chen,B.,Q.Ge,D.Fu,et al.,2010:A data-model fusion approach for upscaling gross ecosystem productivity to the landscape scale based on remote sensing and flux footprint modelling.Biogeosciences,7,2943-2958.

    Chen,B.Z.,J.M.Chen,G.Mo,et al.,2008:Comparison of regional carbon flux estimates from CO2concentration measurements and remote sensing based footprint integration.Global Biogeochem.Cycle,22,doi:10.1029/2007GB003024.

    Chen,B.Z.,N.C.Coops,D.J.Fu,et al.,2011:Assessing eddy-covariance flux tower location bias across the fluxnet-Canada research network based on remote sensing and footprint modelling.Agr.Forest Meteor.,151,87-100.

    Chen,B.Z.,N.C.Coops,D.J.Fu,et al.,2012:Characterizing spatial representativeness of flux tower eddy-covariance measurements across the Canadian carbon program network using remote sensing and footprint analysis.Remote Sens.Environ.,124, 742-755.

    Chen,B.Z.,T.A.Black,N.C.Coops,et al.,2009:Assessing tower flux footprint climatology and scaling between remotely sensed and eddy covariance measurements.Bound.-Layer Meteor.,130,137-167.

    Foken,T.,2008:Micrometeorology.Springer-Verlag, Berlin,Heidelberg,82-87.

    Gash,J.H.C.,1986:A note on estimating the effect of a limited fetch on micrometeorological evaporation measurements.Bound.-Layer Meteor.,35,409-413.

    Horst,T.W.,and J.C.Weil,1992:Footprint estimation for scalar flux measurements in the atmospheric surface layer.Bound.-Layer Meteor.,59,279-296.

    Hsieh,C.I.,G.Katul,and T.W.Chi,2000:An approximate analytical model for footprint estimation of scalar fluxes in thermally stratified atmospheric flows.Adv.Water Resour.,23,765-772.

    Kljun,N.,2010a:Attributing tall tower flux data to heterogeneous vegetation.19th Symposium on Boundary Layers and Turbulence,Keystone,Colorado,2-6 August,Amer.Meteor.Soc.

    Kljun,N.,2010b:BERMS sites revisited:Footprint climatology and 3D-LiDAR data.29th Conference on Agricultural and Forest Meteorology,Keystone, Colorado,2-6 August,Amer.Meteor.Soc.

    Kormann,R.,and F.X.Meixner,2001:An analytical footprint model for non-neutral stratification. Bound.-Layer Meteor.,99,207-224.

    Leclerc,M.Y.,and G.W.Thurtell,1990:Footprint prediction of scalar fluxes using a Markovian analysis. Bound.-Layer Meteor.,52,247-258.

    Leclerc,M.Y.,and T.Foken,2014:Footprints in Micrometeorology and Ecology.Springer-Verlag, Berlin,Heidelberg,71-98.

    Leuning,R.,2005:Measurements of trace gas fluxes in the atmosphere using eddy covariance:WPL corrections revisited.Handbook of Micrometeorology:A Guide for Surface Flux Measurement and Analysis. Springer,Netherlands,119-132.

    Rebmann,C.,M.Gockede,T.Foken,et al.,2005:Quality analysis applied on eddy covariance measurements at complex forest sites using footprint modeling. Theor.Appl.Climatol.,80,121-141.

    Schmid,H.P.,1994:Source areas for scalars and scalar fluxes.Bound.-Layer Meteor.,67,293-318.

    Schmid,H.P.,2002:Footprint modeling for vegetation atmosphere exchange studies:A review and perspective.Agr.Forest Meteor.,113,159-183.

    Sogachev,A.,and J.Lloyd,2004:Using a one-and-a-half order closure model of the atmospheric boundary layer for surface flux footprint estimation.Bound. -Layer Meteor.,112,467-502.

    Tang,Y.K.,X.F.Wen,X.M.Sun,et al.,2014a:The limiting effect of deep soilwater on evapotranspiration of a subtropical coniferous plantation subjected to seasonaldrought.Adv.Atmos.Sci.,31,385-395.

    Tang,Y.K.,X.F.Wen,X.M.Sun,et al.,2014b:Interannual variation of the bowen ratio in a subtropical coniferous plantation in Southeast China, 2003-2012.Plos One,9,e88267,doi:10.1371/journal.pone.0088267.

    Webb,E.K.,G.I.Pearman,and R.Leuning,1980:Correction of flux measurements for density effects due to heat and water vapour transfer.Quart.J.Roy. Meteor.Soc.,106,85-100.

    Wen,X.F.,G.R.Yu,X.M.Sun,et al.,2006:Soil moisture effect on the temperature dependence of ecosystem respiration in a subtropical Pinus plantation of southeastern China.Agr.Forest Meteor.,137,166-175.

    Wen,X.F.,H.M.Wang,J.L.Wang,et al.,2010: Ecosystem carbon exchanges of a subtropical evergreen coniferous plantation subjected to seasonal drought,2003-2007.Biogeosciences,7,357-369.

    Wilczak,J.M.,S.P.Oncley,and S.A.Stage,2001:Sonic anemometer tilt correction algorithms.Bound. -Layer Meteor.,99,127-150.

    Yu Guirui,Wen Xuefa,and Sun Xiaomin,et al.,2006: Overview of ChinaFLUX and evaluation of its eddy covariance measurement.Agr.Forest Meteor.,137, 125-137.

    Zhang Hui,Shen Shuanghe,Wen Xuefa,et al.,2012:Flux footprint of carbon dioxide and vapor exchange over the terrestrialecosystem:A review.Acta Ecol.Sin.,32,7622-7633.(in Chinese)

    Zhang Huiand Wen Xuefa,2015:Flux footprint climatology estimated by three analyticalmodels over a subtropical coniferous plantation in Southeast China.J.Meteor.Res.,29(4),654-666,

    10.1007/s13351-014-4090-7.

    Supported by the National Basic Research and Development(973)Program of China(2012CB416903),National Natural Science Foundation of China(31470500 and 31290221),and Knowledge Innovation Project of the Chinese Academy of Sciences (KZCX2-EW-QN305).

    ?wenxf@igsnrr.ac.cn.

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

    February 6,2015;in final form June 16,2015)

    猜你喜歡
    張慧
    A lightweight symmetric image encryption cryptosystem in wavelet domain based on an improved sine map
    張慧:二十年深耕細作,筑牢國際教育“金字招牌”
    華人時刊(2023年9期)2023-06-20 08:30:52
    幼兒合作意識培養(yǎng)“進行時”
    教育家(2022年17期)2022-04-23 22:21:35
    The Adverse Impact on Immigrant Women in the Workplace
    她放聲歌唱擊退病魔
    飲食科學(2018年3期)2018-06-01 09:28:59
    豆腐中的數(shù)學
    馬坤、張慧、楊帆、鄭玉秀作品
    豆腐中的數(shù)學
    13 Kidney and Urinary Tract
    王洪永?陳英俊?張慧
    免费av中文字幕在线| 久久中文字幕人妻熟女| av在线播放免费不卡| 欧美在线黄色| 国产不卡av网站在线观看| 91精品三级在线观看| 一级片免费观看大全| 久久久国产成人精品二区 | 亚洲精品中文字幕在线视频| 欧美午夜高清在线| 国产精品久久视频播放| 制服人妻中文乱码| 又紧又爽又黄一区二区| 18禁国产床啪视频网站| 一区二区三区精品91| 美女扒开内裤让男人捅视频| 777久久人妻少妇嫩草av网站| tube8黄色片| 久久国产精品大桥未久av| 69精品国产乱码久久久| 高潮久久久久久久久久久不卡| 女同久久另类99精品国产91| 欧美激情 高清一区二区三区| 人成视频在线观看免费观看| 国产亚洲欧美精品永久| 曰老女人黄片| 久久精品aⅴ一区二区三区四区| 国产伦人伦偷精品视频| 亚洲综合色网址| 超碰成人久久| 久久亚洲精品不卡| 亚洲一区二区三区欧美精品| 成人18禁高潮啪啪吃奶动态图| 国产精品98久久久久久宅男小说| 一进一出好大好爽视频| 欧美日韩中文字幕国产精品一区二区三区 | ponron亚洲| av天堂在线播放| 91成人精品电影| av欧美777| 国产亚洲欧美在线一区二区| 交换朋友夫妻互换小说| 国产av又大| 成人国语在线视频| 多毛熟女@视频| 精品第一国产精品| 手机成人av网站| 黑人巨大精品欧美一区二区蜜桃| 老汉色∧v一级毛片| 欧美亚洲日本最大视频资源| 久久精品国产a三级三级三级| 免费少妇av软件| 亚洲,欧美精品.| 国产精品自产拍在线观看55亚洲 | 老司机深夜福利视频在线观看| 久久香蕉精品热| 久久国产乱子伦精品免费另类| 国产主播在线观看一区二区| 国产精品.久久久| 亚洲aⅴ乱码一区二区在线播放 | 男男h啪啪无遮挡| 天天躁夜夜躁狠狠躁躁| 少妇裸体淫交视频免费看高清 | 欧美大码av| 精品少妇一区二区三区视频日本电影| 国产不卡av网站在线观看| 久久久久国内视频| 黄色丝袜av网址大全| 9热在线视频观看99| 麻豆成人av在线观看| 国产男靠女视频免费网站| 午夜福利在线免费观看网站| 少妇的丰满在线观看| 高清欧美精品videossex| 久久久国产欧美日韩av| 亚洲欧美一区二区三区黑人| 成人三级做爰电影| 色老头精品视频在线观看| 日韩 欧美 亚洲 中文字幕| 欧美精品人与动牲交sv欧美| 在线播放国产精品三级| 午夜两性在线视频| 国产亚洲欧美在线一区二区| 色婷婷av一区二区三区视频| 欧美日韩瑟瑟在线播放| 欧美黑人精品巨大| 免费看十八禁软件| 国产精品电影一区二区三区 | 老司机午夜福利在线观看视频| 91九色精品人成在线观看| 精品亚洲成a人片在线观看| 丝袜美腿诱惑在线| 国产成人欧美在线观看 | 国产色视频综合| 亚洲欧美日韩高清在线视频| 成人国产一区最新在线观看| 国产伦人伦偷精品视频| 校园春色视频在线观看| 国产精品久久久久久精品古装| 国产精品二区激情视频| 无人区码免费观看不卡| 9热在线视频观看99| 热99久久久久精品小说推荐| 18禁黄网站禁片午夜丰满| 香蕉久久夜色| 国产精品国产av在线观看| 大香蕉久久成人网| 亚洲 欧美一区二区三区| 亚洲专区中文字幕在线| 狠狠狠狠99中文字幕| 国产精品久久电影中文字幕 | 久久精品国产亚洲av高清一级| 午夜久久久在线观看| 水蜜桃什么品种好| 亚洲熟妇熟女久久| 国产单亲对白刺激| 动漫黄色视频在线观看| 久久久久久久午夜电影 | 午夜福利影视在线免费观看| 国产蜜桃级精品一区二区三区 | 12—13女人毛片做爰片一| 19禁男女啪啪无遮挡网站| 三级毛片av免费| 亚洲专区字幕在线| a级片在线免费高清观看视频| 久久午夜亚洲精品久久| 又黄又粗又硬又大视频| 亚洲成人国产一区在线观看| x7x7x7水蜜桃| 下体分泌物呈黄色| 亚洲中文av在线| 亚洲精品乱久久久久久| 黄网站色视频无遮挡免费观看| 十八禁网站免费在线| 黄色片一级片一级黄色片| 在线观看舔阴道视频| 巨乳人妻的诱惑在线观看| 国产精品久久视频播放| 精品卡一卡二卡四卡免费| 啦啦啦免费观看视频1| 18在线观看网站| 免费观看精品视频网站| 日韩成人在线观看一区二区三区| 国产色视频综合| 亚洲成人手机| 亚洲精品久久成人aⅴ小说| 嫁个100分男人电影在线观看| 精品少妇久久久久久888优播| 日韩制服丝袜自拍偷拍| 一级,二级,三级黄色视频| 淫妇啪啪啪对白视频| 丝瓜视频免费看黄片| 欧美日韩乱码在线| 免费观看精品视频网站| 日韩精品免费视频一区二区三区| 99精品久久久久人妻精品| 9191精品国产免费久久| 国产不卡av网站在线观看| 免费在线观看黄色视频的| 中文字幕精品免费在线观看视频| 捣出白浆h1v1| 久久久国产一区二区| 亚洲成人免费电影在线观看| 欧美日韩亚洲综合一区二区三区_| 欧美另类亚洲清纯唯美| 麻豆成人av在线观看| 丰满的人妻完整版| 香蕉丝袜av| 99riav亚洲国产免费| 中文字幕色久视频| netflix在线观看网站| 建设人人有责人人尽责人人享有的| 美女午夜性视频免费| 欧美激情高清一区二区三区| 久久这里只有精品19| 国精品久久久久久国模美| 精品一区二区三区视频在线观看免费 | 50天的宝宝边吃奶边哭怎么回事| 九色亚洲精品在线播放| 精品一区二区三卡| 亚洲一码二码三码区别大吗| 成人国语在线视频| 亚洲一区二区三区不卡视频| 成人国产一区最新在线观看| 国产一区二区激情短视频| 欧美人与性动交α欧美精品济南到| 亚洲色图综合在线观看| 日韩欧美一区视频在线观看| 亚洲免费av在线视频| 一二三四在线观看免费中文在| 捣出白浆h1v1| www.精华液| 少妇粗大呻吟视频| 久久精品91无色码中文字幕| 精品卡一卡二卡四卡免费| 久久精品aⅴ一区二区三区四区| 欧美日韩一级在线毛片| 黄色成人免费大全| 欧美日韩亚洲综合一区二区三区_| 亚洲欧洲精品一区二区精品久久久| 久久久久久人人人人人| 久久亚洲真实| 亚洲欧美一区二区三区久久| 日本黄色日本黄色录像| 正在播放国产对白刺激| 波多野结衣av一区二区av| 成人手机av| 高清欧美精品videossex| 一本一本久久a久久精品综合妖精| 在线十欧美十亚洲十日本专区| 免费在线观看完整版高清| 99久久精品国产亚洲精品| 在线观看免费视频日本深夜| 老司机影院毛片| 夫妻午夜视频| 亚洲美女黄片视频| 女警被强在线播放| 国产99白浆流出| 黄频高清免费视频| 99精品在免费线老司机午夜| 一区福利在线观看| 99热网站在线观看| 免费在线观看影片大全网站| 麻豆av在线久日| 好看av亚洲va欧美ⅴa在| 久久天堂一区二区三区四区| 99久久综合精品五月天人人| 国产深夜福利视频在线观看| 国产色视频综合| www日本在线高清视频| 捣出白浆h1v1| av天堂久久9| 国产日韩欧美亚洲二区| 久久香蕉激情| 日本黄色视频三级网站网址 | 亚洲成人免费电影在线观看| 日本a在线网址| 免费观看精品视频网站| 亚洲成人国产一区在线观看| 美女扒开内裤让男人捅视频| 日韩成人在线观看一区二区三区| 无人区码免费观看不卡| 亚洲精品中文字幕在线视频| 99riav亚洲国产免费| 一本大道久久a久久精品| 深夜精品福利| 成人18禁在线播放| 丁香欧美五月| 日韩一卡2卡3卡4卡2021年| 一区二区三区激情视频| 亚洲欧美一区二区三区黑人| 91九色精品人成在线观看| 天堂动漫精品| av天堂久久9| 搡老岳熟女国产| 精品熟女少妇八av免费久了| 国产一区二区三区视频了| 人成视频在线观看免费观看| av不卡在线播放| 国产精品.久久久| 制服人妻中文乱码| 80岁老熟妇乱子伦牲交| 91麻豆av在线| tube8黄色片| 在线观看日韩欧美| 久久婷婷成人综合色麻豆| 人成视频在线观看免费观看| 韩国av一区二区三区四区| 十八禁网站免费在线| 多毛熟女@视频| 村上凉子中文字幕在线| 很黄的视频免费| 动漫黄色视频在线观看| 国产成人影院久久av| www日本在线高清视频| 国产成人精品在线电影| 看片在线看免费视频| 黄片小视频在线播放| 高清av免费在线| 人人妻,人人澡人人爽秒播| 如日韩欧美国产精品一区二区三区| 手机成人av网站| 亚洲性夜色夜夜综合| 免费在线观看日本一区| www.999成人在线观看| 夫妻午夜视频| 国产高清激情床上av| av片东京热男人的天堂| 韩国精品一区二区三区| 国产国语露脸激情在线看| 亚洲熟女毛片儿| 在线看a的网站| 一区二区三区国产精品乱码| 成人免费观看视频高清| 国产成人免费无遮挡视频| 黄片小视频在线播放| 日韩欧美在线二视频 | 欧美国产精品一级二级三级| 亚洲精品久久成人aⅴ小说| 18禁裸乳无遮挡动漫免费视频| 久久人人97超碰香蕉20202| 制服诱惑二区| 热re99久久精品国产66热6| 免费人成视频x8x8入口观看| 亚洲午夜精品一区,二区,三区| 国内毛片毛片毛片毛片毛片| 搡老熟女国产l中国老女人| 五月开心婷婷网| 国产aⅴ精品一区二区三区波| av欧美777| 大码成人一级视频| 欧美av亚洲av综合av国产av| 色婷婷av一区二区三区视频| 亚洲少妇的诱惑av| 国产麻豆69| 一进一出抽搐gif免费好疼 | 久久国产亚洲av麻豆专区| 99精品欧美一区二区三区四区| 日韩欧美在线二视频 | 夜夜躁狠狠躁天天躁| 国产av一区二区精品久久| 精品久久蜜臀av无| 国产真人三级小视频在线观看| 操美女的视频在线观看| 精品亚洲成a人片在线观看| 视频区图区小说| 成人av一区二区三区在线看| 国产欧美日韩一区二区精品| 成人av一区二区三区在线看| 男女下面插进去视频免费观看| 久久午夜综合久久蜜桃| 超碰97精品在线观看| 久久国产精品大桥未久av| 国产男女内射视频| 动漫黄色视频在线观看| 大陆偷拍与自拍| 伦理电影免费视频| 国产精品综合久久久久久久免费 | 亚洲专区字幕在线| 欧美激情 高清一区二区三区| 老司机午夜十八禁免费视频| av免费在线观看网站| 成人亚洲精品一区在线观看| 亚洲国产中文字幕在线视频| 老熟妇仑乱视频hdxx| 免费观看人在逋| 精品少妇久久久久久888优播| 免费少妇av软件| 别揉我奶头~嗯~啊~动态视频| 妹子高潮喷水视频| tocl精华| 久久久久久人人人人人| 人人妻人人澡人人看| 18禁观看日本| 美女 人体艺术 gogo| 91字幕亚洲| 国产精品美女特级片免费视频播放器 | 国产精品一区二区免费欧美| 十八禁网站免费在线| 人妻丰满熟妇av一区二区三区 | 国产激情久久老熟女| 国产淫语在线视频| 法律面前人人平等表现在哪些方面| 精品一区二区三卡| 亚洲avbb在线观看| 黄色女人牲交| 黑人欧美特级aaaaaa片| 久久久水蜜桃国产精品网| 一级毛片女人18水好多| 亚洲av日韩精品久久久久久密| 啦啦啦视频在线资源免费观看| 啦啦啦 在线观看视频| 啦啦啦在线免费观看视频4| 啦啦啦 在线观看视频| 一级毛片精品| 午夜福利免费观看在线| 91九色精品人成在线观看| 亚洲九九香蕉| 色老头精品视频在线观看| 嫁个100分男人电影在线观看| 日韩中文字幕欧美一区二区| 国产熟女午夜一区二区三区| ponron亚洲| 国产欧美日韩综合在线一区二区| 欧美+亚洲+日韩+国产| 国产精品秋霞免费鲁丝片| av天堂在线播放| 国产成人免费观看mmmm| 啦啦啦免费观看视频1| av网站在线播放免费| 美女午夜性视频免费| 午夜久久久在线观看| av视频免费观看在线观看| 成人18禁在线播放| 深夜精品福利| 国产1区2区3区精品| 露出奶头的视频| 在线国产一区二区在线| 精品国产乱码久久久久久男人| 丝袜美足系列| 亚洲av第一区精品v没综合| www日本在线高清视频| 免费人成视频x8x8入口观看| 亚洲欧美激情在线| 中文字幕人妻熟女乱码| 女人高潮潮喷娇喘18禁视频| 成人三级做爰电影| 国产一卡二卡三卡精品| 看免费av毛片| 黄色视频,在线免费观看| 又黄又粗又硬又大视频| 日韩人妻精品一区2区三区| 麻豆乱淫一区二区| 久久香蕉激情| 国产野战对白在线观看| 成熟少妇高潮喷水视频| 一进一出好大好爽视频| 亚洲色图综合在线观看| 久久中文字幕人妻熟女| 亚洲色图av天堂| 99riav亚洲国产免费| 亚洲精品国产区一区二| 精品无人区乱码1区二区| 亚洲少妇的诱惑av| 熟女少妇亚洲综合色aaa.| 国产精品国产av在线观看| 亚洲在线自拍视频| 亚洲av电影在线进入| 国产亚洲精品一区二区www | 老熟妇仑乱视频hdxx| 老司机午夜十八禁免费视频| 久久狼人影院| 久久影院123| 亚洲精品美女久久av网站| 男男h啪啪无遮挡| 91av网站免费观看| 久热爱精品视频在线9| 中文欧美无线码| 午夜精品国产一区二区电影| 国产单亲对白刺激| 国产不卡av网站在线观看| 看免费av毛片| 91av网站免费观看| 一级黄色大片毛片| 日韩有码中文字幕| 亚洲美女黄片视频| 久久久久精品国产欧美久久久| 欧美黄色片欧美黄色片| 亚洲成人手机| 人成视频在线观看免费观看| 999久久久精品免费观看国产| 男女午夜视频在线观看| 国产麻豆69| 巨乳人妻的诱惑在线观看| 黄色丝袜av网址大全| 天堂√8在线中文| 一级毛片高清免费大全| 久久99一区二区三区| 欧美成人午夜精品| 国产国语露脸激情在线看| 久久国产乱子伦精品免费另类| 欧美中文综合在线视频| 十八禁网站免费在线| 亚洲午夜理论影院| 国产精品久久久久久人妻精品电影| 亚洲五月色婷婷综合| 免费观看a级毛片全部| 午夜免费成人在线视频| 亚洲精品国产一区二区精华液| 精品乱码久久久久久99久播| 亚洲午夜理论影院| a级毛片在线看网站| 国产99白浆流出| 亚洲欧美激情综合另类| 午夜福利乱码中文字幕| 黄色怎么调成土黄色| 黄色视频,在线免费观看| 欧美日韩精品网址| 欧美日韩一级在线毛片| 久久国产精品男人的天堂亚洲| 在线播放国产精品三级| 女人久久www免费人成看片| 久久精品亚洲熟妇少妇任你| 十八禁人妻一区二区| 欧美+亚洲+日韩+国产| 成人三级做爰电影| 一二三四社区在线视频社区8| 女人爽到高潮嗷嗷叫在线视频| a级毛片在线看网站| 女人被躁到高潮嗷嗷叫费观| 99久久综合精品五月天人人| 男女床上黄色一级片免费看| 国产在视频线精品| 免费看十八禁软件| 老熟女久久久| 色尼玛亚洲综合影院| 欧美激情极品国产一区二区三区| 欧美日韩瑟瑟在线播放| 成人黄色视频免费在线看| 女人高潮潮喷娇喘18禁视频| 国产人伦9x9x在线观看| 老司机影院毛片| 国产精品98久久久久久宅男小说| 国产极品粉嫩免费观看在线| 国产一区在线观看成人免费| 精品少妇久久久久久888优播| 日韩欧美免费精品| 国产视频一区二区在线看| 国产成人一区二区三区免费视频网站| 国产主播在线观看一区二区| 久久 成人 亚洲| 亚洲一码二码三码区别大吗| 欧美 日韩 精品 国产| 亚洲第一欧美日韩一区二区三区| 亚洲,欧美精品.| 色播在线永久视频| 亚洲精品自拍成人| 99国产极品粉嫩在线观看| 国产精品永久免费网站| 女性生殖器流出的白浆| 国产精品98久久久久久宅男小说| 黄色片一级片一级黄色片| 亚洲片人在线观看| 成人精品一区二区免费| 日韩视频一区二区在线观看| 天天躁日日躁夜夜躁夜夜| 久久中文字幕一级| 中文字幕高清在线视频| 国产精品九九99| 亚洲中文日韩欧美视频| 免费高清在线观看日韩| 午夜福利影视在线免费观看| 亚洲性夜色夜夜综合| 久久精品亚洲精品国产色婷小说| 国内毛片毛片毛片毛片毛片| 成人永久免费在线观看视频| 精品国产一区二区三区久久久樱花| 国产成+人综合+亚洲专区| 日韩欧美免费精品| 国产精品乱码一区二三区的特点 | 欧美日韩黄片免| 香蕉国产在线看| 欧美日韩亚洲高清精品| 热99re8久久精品国产| 婷婷精品国产亚洲av在线 | 国产在线观看jvid| 欧美成人免费av一区二区三区 | 中文欧美无线码| 久久久久久人人人人人| 三级毛片av免费| 精品亚洲成国产av| 成人三级做爰电影| 日韩欧美一区二区三区在线观看 | 精品人妻熟女毛片av久久网站| 免费日韩欧美在线观看| 午夜精品久久久久久毛片777| 久久精品熟女亚洲av麻豆精品| 国产1区2区3区精品| 午夜福利在线观看吧| 日本五十路高清| 久久ye,这里只有精品| 丝袜美腿诱惑在线| 婷婷成人精品国产| 自拍欧美九色日韩亚洲蝌蚪91| 在线天堂中文资源库| 欧美日韩av久久| 很黄的视频免费| 免费在线观看完整版高清| 两性午夜刺激爽爽歪歪视频在线观看 | 免费在线观看亚洲国产| 精品乱码久久久久久99久播| 欧美性长视频在线观看| 久久久国产成人免费| 国产高清激情床上av| netflix在线观看网站| 国产精品成人在线| 精品国产一区二区三区四区第35| 国精品久久久久久国模美| a级毛片在线看网站| 涩涩av久久男人的天堂| 人妻丰满熟妇av一区二区三区 | 一本大道久久a久久精品| 男人操女人黄网站| 一区二区日韩欧美中文字幕| 韩国av一区二区三区四区| 欧美成狂野欧美在线观看| 久久精品熟女亚洲av麻豆精品| 亚洲国产欧美网| 国产麻豆69| 黄片小视频在线播放| 欧美精品啪啪一区二区三区| 国产欧美日韩综合在线一区二区| 国产不卡av网站在线观看| 老司机深夜福利视频在线观看| 亚洲人成电影免费在线| 国产欧美日韩一区二区三| 超碰成人久久| 黄色片一级片一级黄色片| 亚洲第一青青草原| 又黄又爽又免费观看的视频| 亚洲人成伊人成综合网2020| 男女床上黄色一级片免费看| 老司机在亚洲福利影院| 国产区一区二久久| 亚洲精品粉嫩美女一区| 欧美成人午夜精品| 自拍欧美九色日韩亚洲蝌蚪91| 91成年电影在线观看| 妹子高潮喷水视频| 多毛熟女@视频| 欧美日韩乱码在线| 999久久久国产精品视频| 最新美女视频免费是黄的| 亚洲aⅴ乱码一区二区在线播放 | 免费观看人在逋|