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

    Response of Growing Season Gross Primary Production to El Ni?o in Different Phases of the Pacific Decadal Oscillation over Eastern China Based on Bayesian Model Averaging

    2021-07-26 14:38:42YueyueLILiDANJingPENGJunbangWANGFuqiangYANGDongdongGAOXiujingYANGandQiangYU
    Advances in Atmospheric Sciences 2021年9期

    Yueyue LI ,Li DAN ,Jing PENG ,Junbang WANG ,Fuqiang YANG ,Dongdong GAO ,Xiujing YANG,and Qiang YU

    1Key Laboratory of Water Cycle and Related Land Surface Processes,Institute of Geographic Sciences and Natural Resources Research,Chinese Academy of Sciences,Beijing 100101,China

    2College of Resources and Environment,University of Chinese Academy of Science,Beijing 100049,China

    3Key Laboratory of Regional Climate-Environment for Temperate East Asia,Institute of Atmospheric Physics,Chinese Academy of Sciences,Beijing 100029,China

    4National Ecosystem Science Data Center,Key Laboratory of Ecosystem Network Observation and Modeling,Institute of Geographic Sciences and Natural Resources Research,Chinese Academy of Sciences,Beijing 100101,China

    5State Key Laboratory of Soil Erosion and Dryland Farming on the Loess Plateau,Northwest A & F University,Shaanxi 712100,China

    6School of Life Sciences,University of Technology Sydney,Sydney 2007,NSW,Australia

    ABSTRACT Gross primary production (GPP) plays a crucial part in the carbon cycle of terrestrial ecosystems.A set of validated monthly GPP data from 1957 to 2010 in 0.5° × 0.5° grids of China was weighted from the Multi-scale Terrestrial Model Intercomparison Project using Bayesian model averaging (BMA).The spatial anomalies of detrended BMA GPP during the growing seasons of typical El Ni?o years indicated that GPP response to El Ni?o varies with Pacific Decadal Oscillation(PDO) phases:when the PDO was in the cool phase,it was likely that GPP was greater in northern China (32°–38°N,111°–122°E) and less in the Yangtze River valley (28°–32°N,111°–122°E);in contrast,when PDO was in the warm phase,the GPP anomalies were usually reversed in these two regions.The consistent spatiotemporal pattern and high partial correlation revealed that rainfall dominated this phenomenon.The previously published findings on how El Ni?o during different phases of PDO affecting rainfall in eastern China make the statistical relationship between GPP and El Ni?o in this study theoretically credible.This paper not only introduces an effective way to use BMA in grids that have mixed plant function types,but also makes it possible to evaluate the carbon cycle in eastern China based on the prediction of El Ni?o and PDO.

    Key words:East China,Bayesian model averaging,Gross primary production,El Ni?o,Pacific Decadal Oscillation,Monsoon rainfall

    1.Introduction

    Gross primary production (GPP) is the total amount of carbon dioxide "fixed" by land plants per unit time through the photosynthetic reduction of COinto organic compounds (

    Gough,2011

    ).The capacity of terrestrial GPP to partly offset anthropogenic COemissions is especially important for China,where the emissions of fossil fuel COinto the atmosphere are high (

    Gregg et al.,2008

    ).The subtropical forests in eastern China can take in a considerable amount of CO(

    Yu et al.,2014

    ),and contribute remarkably to the total carbon sequestration of Chinese ecosystems(

    Zhang et al.,2019a

    ).The climate in eastern China is dominated by the East Asian Summer Monsoon (EASM),which is one of the most active components in the global climate system and can significantly affect the movement of local rain belts (

    Ding et al.,2008

    ).Given the ecological importance of the monsoon climate in China and the high sensitivity of the carbon uptake of ecosystems in EASM region to rainfall(

    Dan et al.,2015

    ;

    Zhang et al.,2019a

    ),it is important to develop a deeper understanding of the response of GPP to the variability in rainfall during the EASM.Different modes of sea surface temperature (SST) can explain most of the inter-annual and inter-decadal variations of rainfall in China (

    Yang and Lau,2004

    ).Two modes of SST that reflect strong climatic teleconnections with the EASM rainfall are El Ni?o and the Pacific Decadal Oscillation (PDO) (

    Ma,2007

    ).An El Ni?o event is characterized by anomalous warming in the central and eastern equatorial Pacific Ocean and occurs,on average,every two to seven years.By contrast,the PDO (

    Mantua et al.,1997

    ;

    Qian and Zhou,2014

    ) is marked by an inter-decadal cycle,and is mainly present in the North Pacific (poleward of 20°N).During the warm phase of the PDO,SST tends to be anomalously cool in the central North Pacific,coincident with anomalously warm SST in the eastern North Pacific,whereas the spatial pattern of SST is reversed during the cool phase of the PDO.Because the inter-annual variability in the atmospheric concentrations of COinduced by El Ni?o is mainly attributed to tropical terrestrial ecosystems (

    Yang and Wang,2000

    ;

    Cox et al.,2013

    ;

    Wang et al.,2013

    ;

    Liu et al.,2017

    ),previous studies have primarily focused on how terrestrial carbon fluxes in the tropics responding to changes in the Earth’s climate (especially changes in temperature and rainfall)caused by El Ni?o (

    Gu and Adler,2011

    ;

    Bastos et al.,2013

    ;

    Wang et al.,2016

    ;

    Fang et al.,2017

    ).The most common way to investigate the decadal variability in the terrestrial carbon budget caused by the PDO is to use the long-term teleconnection between the PDO indices and specific carbon cycle variables (

    Wharton and Falk,2016

    ;

    Zhang et al.,2018

    ).However,this method is unable to reflect the impact of the variability of local chaotic weather events (

    Ito,2011

    ),especially when the responses of terrestrial carbon cycle to different El Ni?o events are contrasting (

    Liu et al.,2017

    ;

    Wang et al.,2018

    ).Hence,the mixed effects of El Ni?o and the PDO on the carbon cycle in the extratropical monsoon regions of China are still poorly understood.Even though the spatiotemporal patterns of rainfall anomalies over East Asia associated with individual El Ni?o episodes are variable (

    Wang et al.,2017

    ),

    Feng et al.(2014)

    still found relatively stable relationships between the EASM rainfall and El Ni?o under different phases of the PDO.When the PDO is in a cool phase,with a high-pressure area in the subtropical western Pacific clearly experiencing two northward shifts,positive rainfall anomalies move from southern to northern China.By contrast,when the PDO is in a warm phase,a greater amount of rainfall is observed in central China,whereas less rainfall is observed in southern and northern China.How GPP responds to these specific rainfall patterns is still unclear.

    The purpose of this study was to determine the relationship between the GPP and the changes in rainfall associated with El Ni?o under different phases of the PDO in eastern China.In particular,we attempted to identify the regions driving this relationship and to demonstrate the response of GPP to El Ni?o and the PDO.The results may have implications for how terrestrial carbon dynamics responding to variabilities in the Earth’s climate.

    2.Materials and Methods

    2.1.MsTMIP GPP and CRU–NCEP Datasets

    The rainfall patterns described in the fourth paragraph of the Introduction were revealed and explained by data spanning from 1957 to 2010.Given the data consistency,this study also focused on the same period and selected the same El Ni?o years (see

    Table 1

    ) as

    Feng et al.(2014)

    .Compared with other approaches,the Terrestrial Biosphere Model (TMB) is a preferred method to estimate carbon fluxes over such a long period because of the availability of climate forcing data and process-based simulation.The monthly GPP products of 15 TMBs in the Multi-scale synthesis and Terrestrial Model Intercomparison Project (MsTMIP) (

    Huntzinger et al.,2013

    ;

    Wei et al.,2014

    ) were applied here.The MsTMIP version 1.0 dataset includes four global (0.5° × 0.5° resolution) sensitivity simulations (SG1,SG2,SG3 and BG1) which can systematically assess the impacts of different forcing factors.The model-driven data in MsTMIP version 1.0 are climate variables (Climate),land use/land cover changes (LUCC),the atmospheric COconcentration (CO) and nitrogen deposition (N).The timevarying drivers in the four global sensitivity simulations are Climate (SG1),Climate+LUCC (SG2),Climate+LUCC +CO(SG3),and Climate+LUCC+CO+N (BG1).BG1 only includes a limited number of models as a result of the lack of a nitrogen cycle in other MsTMIP members.SG3 and BG1 are presumed to be closer to reality than SG1 and SG2 and were therefore used in this study.The CLASS-CTEMN model was excluded because its GPP estimates for China are much lower than is reasonable (

    Shao et al.,2016

    ).

    Table 1.Characteristics of the 16 El Ni?o years in this study.

    The available climatology datasets from the Climate Research Unit (CRU) (

    Harris et al.,2014

    ) and the NationalCenters for Environmental Prediction (NCEP)/National Center for Atmospheric Research (NCAR) (

    Kalnay et al.,1996

    )cannot fully meet the spatial and temporal requirements of MsTMIP.Therefore,a new dataset named CRU-NCEP was produced by combining the strengths of the datasets from the CRU and NCEP/NCAR as driver data for MsTMIP.The CRU–NCEP dataset not only maintains a good match to the observed rainfall in China (

    Zhao and Fu,2006

    ),but also corrects the known biases in temperature and downward shortwave radiation in the NCEP/NCAR reanalysis products(

    Wei et al.,2014

    ).Therefore,the globally gridded (0.5° ×0.5°) and monthly CRU-NCEP dataset was chosen to determine the relationships between GPP and the particular meteorological factors (rainfall,temperature,and downward shortwave radiation).

    2.2.Bayesian Model Averaging

    The optimum of the likelihood function is:

    After implementation of training with data

    y

    ,βfrom the solutions of Eq.(2) should satisfy the conditions:

    where

    w

    =βfrac+βfrac+βfrac.The BMA GPP in the target grid cell was given as:

    2.3.Data Preprocessing

    To highlight the impact of the EASM on GPP,May and September were defined as the start and end months of the growing season corresponding to the onset and withdrawal months of the EASM.When analyzing the response characteristics of GPP to El Ni?o over China and the corresponding mechanism,GPP and the meteorological data were not only limited to the growing season in each year,but also were detrended to minimize the possible influence of long-term trends.Trends were calculated by simple linear regression.

    3.Results

    3.1.Validation of the BMA GPP at the ChinaFLUX Sites

    The MsTMIP and BMA GPP data in the grid cells containing ChinaFLUX sites from 2009 to 2010 were chosen to evaluate the rationality of using BMA.The efficiency criteria used in this section are Eindices (see Appendix A.4).The Eindices for different models are listed in

    Table 2

    .BMA performed better than the MsTMIP models in the qualified grid cells corresponding to Changbaishan,Qianyanzhou,Yucheng,and Inner Mongolia,thereby proving its strength for the typical PFTs at these sites.Even though the Eindices calculated in the grid cells for Dinghushan and Haibei were not statistically significant,the relatively small absolute values still indicated that the magnitude of the BMA GPP was acceptable,especially in the subtropical region where Dinghushan is located.The evergreen broad leaf trees in Xishuangbanna and the grasses in Dangxiong are rainforests and alpine meadows,respectively,and the negative Eindexes suggested that the manner in which BMA was used in this study was unable to reduce the uncertainties from these two specific PFTs.

    Table 2.Nash–Sutcliffe efficiency index values for GPP simulated by 14 terrestrial biosphere models and Bayesian model averaging at eight ChinaFLUX sites.

    3.2.The Response of GPP to El Ni?o in Different Phases of the PDO

    Figure 1

    shows the average spatial anomalies of the BMA GPP and the CRU–NCEP rainfall for the selected El Ni?o years with different phases of the PDO.The GPP and rainfall with large anomalies were mainly distributed in northern China (32°–38°N,111°–122°E) and the Yangtze River valley (28°–32°N,111°–122°E).The similar spatial distribution of the GPP and the rainfall anomalies shown in

    Fig.1

    indicate that rainfall is the connection between the El Ni?o and GPP in eastern China (28°–38°N,111°–122°E).The spatial anomalies of the GPP averaged from multiple El Ni?o years (

    Fig.1a

    ,

    Fig.1b

    ) leads to the following hypotheses:(1) the zone in which GPP was significantly affected by El Ni?o was in eastern China (28°–38°N,111°–122°E);(2) when the El Ni?o years were in the cool phase of the PDO,GPP was higher in northern China(32°–38°N,111°–122°E) and lower in the Yangtze River valley (28°–32°N,111°–122°E);and (3) when the El Ni?o years were in the warm phase of the PDO,the spatial anomaly patterns of the GPP in these two regions were reversed.To check whether these hypotheses are valid,the response of GPP to each El Ni?o year was analyzed.

    Figure 2

    shows the spatial anomalies of GPP in eight typical El Ni?o years under the cool phase of the PDO.Except for 1968 and 1994,the graphs for other years in

    Fig.2

    all support these hypotheses and the most positive anomalies in northern China (32°–38°N,111°–122°E) were statistically significant.

    Figure 3

    shows the spatial anomalies of GPP in eight typical El Ni?o years under the warm phase of the PDO.Excluding 1957 and 1976,the positive anomalies of GPP for the remaining El Ni?o years were mainly distributed in the Yangtze River valley (28°–32°N,111°–122°E),which further verified the hypotheses from the averaged El Ni?o years.Less significant positive GPP anomalies were found in the Yangtze River valley.This difference suggests that plants in northern China are more sensitive to El Ni?o than those in the Yangtze River valley.The evolution of the zonal-time anomaly (May–September) of rainfall and GPP averaged from 111° to 122°E for typical El Ni?o years in the cool and warm phases of the PDO are shown in

    Fig.4

    and

    Fig.5

    .Even when there are some deflections,it is still clear that the large anomalies in the rainfall and GPP are synchronized both spatially and temporally.This confirms that the close relationship between GPP and El Ni?o in eastern China is related to the north–south movement of the rain belt.

    3.3.The Solid Relationship between GPP and Rainfall in Eastern China

    The inter-annual relationships between meteorological factors and vegetation activity can be changed in different dec-ades (

    Piao et al.,2014

    ).The EASM weakened from the late 1970s but has gradually recovered since the early 1990s(

    Liu et al.,2012

    ).Two abrupt climate change points are recognized in 1978 and 1992 (

    Ding et al.,2008

    ).To determine whether the relationships between GPP and the particular meteorological factors (rainfall,temperature,and the downward shortwave radiation) are consistent under different EASM intensities,they were compared in different time periods (1957–78,1979–92,and 1993–2010).

    Fig.1.Spatial anomalies of BMA GPP during the growing season (May–September) averaged from El Ni?o years under (a) the cool phase of the PDO and (b) the warm phase of the PDO in China.(c) and (d) are the same as (a) and (b) respectively,but for rainfall.The anomalies are relative to the mean of the time period 1957–2010.The dashed lines denote negative anomalies and the colored areas denote positive anomalies.The target zones (28°–32°N,111°–122°E;32°–38°N,111°–122°E) are shown by the rectangles.GPP and rainfall were detrended in advance.

    To confirm the role of key meteorological factors for GPP,partial correlations between GPP and rainfall,downward shortwave radiation,and temperature in China are given in

    Fig.6

    .This analysis shows that rainfall dominated the inter-annual changes of the GPP in eastern China from 1957 to 2010.By contrast,the roles of temperature and downward shortwave radiation in the same region were rather limited,although the effect of temperature was clearer in the time period 1957–78.The bar diagrams in

    Fig.7

    show the proportions of partial correlations between GPP and the particular meteorological factors in northern China and the Yangtze River valley,respectively.The highest proportions of positive partial correlations are between GPP and rainfall in all periods in

    Fig.7

    ,which clearly shows the importance of rainfall.The results in this section further enhance the credibility of the response characteristics of GPP to El Ni?o in eastern China being mostly due to its close relationship with rainfall.

    4.Discussion

    The mismatch in spatial scale between the observed GPP at the flux sites and the modeled GPP from MsTMIP in the corresponding 0.5° × 0.5° grids may lead to statistical errors.The monthly GPP from the MODerate-resolution Imaging Spectroradiometer (MODIS) in 0.05° × 0.05°grids and machine learning Model Tree Ensemble (MTE) in 0.5° × 0.5°grids (see Appendix A.5) were chosen to validate the quality of the BMA GPP with such errors.From the scatterplots averaged from northern China and the Yangtze River valley (see

    Fig.8

    ),we found that although the amplitude of the BMA GPP anomaly was less than the MODIS GPP anomaly,their Rvalues all passed 99% confidence tests.Compared with the single model in MsTMIP,BMA GPP is more similar to the multi-year averaged distribution of MTE GPP,especially in eastern China (see

    Fig.9

    ).This spatiotemporal consistency between BMA GPP and validation data in

    Fig.8

    and

    Fig.9

    not only proves the effectiveness of validating simulations in coarse grids with site observations,but also shows the rationality of simply classifying PFTs as trees,grasses,and crops when using BMA.There are no rainforests or alpine meadows in eastern China,so the uncertainties from these ecosystems can be neglected.We conclude that the quality of the BMA GPP in eastern China is reliable.

    Fig.2.Spatial anomalies of BMA GPP during the growing season (May–September) for eight selected El Ni?o years under the cool phase of the PDO in China.The anomalies are relative to the mean of the time period 1957–2010.The blank areas denote negative anomalies and the colored areas denote positive anomalies.The target zones (28°–32°N,111°–122°E;32°–38°N,111°–122°E) are shown by the rectangles.σ is the standard deviation of the BMA GPP anomalies,where ≥1.5σ means that the positive anomalies are significant.GPP was detrended in advance.

    Fig.3.Spatial anomalies of BMA GPP during the growing season (May–September) for eight selected El Ni?o years under the warm phase of the PDO in China.The anomalies are relative to the mean of the time period 1957–2010.The blank areas denote negative anomalies and the colored areas denote positive anomalies.The target zones (28°–32°N,111°–122°E;32°–38°N,111°–122°E) are shown by the rectangles.σ is the standard deviation of the BMA GPP anomalies,where ≥1.5σ means that the positive anomalies are significant.GPP were detrended in advance.

    Fig.4.Latitude–time cross-sections of the (a) rainfall (mm yr?1) and (b) BMA GPP (gC m?2 yr?1) anomalies during the growing season (May–September) for El Ni?o years under the cool phase of the PDO in eastern China(28°–38°N,111°–122°E).The dashed lines denote negative anomalies and the colored areas denote positive anomalies.GPP and rainfall were detrended in advance.

    Fig.5.Latitude–time cross-sections of the (a) rainfall (mm yr?1) and (b) BMA GPP (gC m?2 yr?1) anomalies during the growing season (May–September) for El Ni?o years under the warm phase of the PDO in eastern China(28°–38°N,111°–122°E).The dashed lines denote negative anomalies and the colored areas denote positive anomalies.GPP and rainfall were detrended in advance.

    Fig.6.Spatial patterns of partial correlations between the growing season (May–September) BMA GPP and (a) rainfall(RAIN),(b) downward shortwave radiation (DWSR),and (c) temperature (TEMP) in China in 1957–78.(d),(e),and (f)are the same as (a),(b),and (c),respectively,but in 1979–92.(g),(h),and (i) are also the same as (a),(b),and (c),respectively,but in 1993–2010.Gray regions indicate insignificant partial correlations (P >0.05).The target zones(28°–32°N,111°–122°E;32°–38°N,111°–122°E) are shown by rectangles.The GPP and meteorological factors were detrended in advance.

    Fig.7.Proportions of positive (+) and negative (–) partial correlations between the growing season(May–September) BMA GPP and rainfall (P),temperature (T),and downward shortwave radiation (R) in (a)northern China (32°–38°N,111°–122°E) and (b) the Yangtze River valley (28°–32°N,111°–122°E) in strong and weak East Asian Summer Monsoon periods (1957–78,1979–92,and 1993–2010).Gray bars indicate the proportions of insignificant partial correlations.All data were detrended in advance.

    Fig.8.Anomaly scatterplots of the monthly MODIS GPP against the monthly BMA GPP averaged in (a)northern China (32°–38°N,111°–122°E) and (b) the Yangtze River valley (28°–32°N,111°–122°E) in the time period 2000–10.All of the GPP data were detrended in advance.

    The general consensus is that GPP in the tropics is usually reduced due to the warm and dry conditions caused by El Ni?o (

    Cavaleri et al.,2017

    ;

    Yang et al.,2018

    ;

    Qian et al.,2019

    ;

    Yan et al.,2019

    ).By contrast,the current study was carried out in the temperate zone (28°–38°N,111°–122°E) of China.We found that GPP was not simply increased or decreased,but showed opposite north–south anomalies(divided by 32°N) when El Ni?o occurred under different phases of the PDO.However,the spatial anomalies of the BMA GPP in eastern China for 1957,1968,1976,and 1994 did not follow the general pattern seen in most selected El Ni?o years.This inconsistency suggests that using a combination of El Ni?o and the phase of the PDO to predict GPP is not completely reliable.The misplacement of the BMA GPP anomalies in eastern China for 1968 (minor El Ni?o)compared with 2006 (minor El Ni?o),1994 (moderate El Ni?o) compared with 2004 (moderate El Ni?o),and 1957(major El Ni?o) compared with 1991 (major El Ni?o) also indicates that a classification of the El Ni?o event as“strong” or “weak” is not sufficient to judge the interannual variations in the GPP.The diversities in how GPP responds to different El Ni?o episodes further complement and deepen our understanding of the nonlinear climatic effects of sea surface temperature on local plants.The inter-annual variabilities of the GPP derived from eddy covariance network observations are controlled by rainfall in temperate regions (

    Jung et al.,2011

    ),which is consistent with the results in this study.The different responses of GPP to El Ni?o in tropical and temperate regions are seen not only for the enzyme-driven land surface models in this paper,but are also reported for the model driven by remote sensing data (

    Zhang et al.,2019b

    ).Specific terrestrial biosphere models (e.g.,DLEM,LPJ,VEGAS,VISIT,CLM,ISAM,and ORCHIDEE) with a proven ability to reproduce the deviations in carbon flux caused by El Ni?o (

    Wang et al.,2016

    ;

    Chang et al.,2017

    ;

    Bastos et al.,2018

    ) also took part in MsTMIP and contributed to the BMA GPP.Therefore,our conclusions based on the BMA GPP are credible.The GPP in eastern China is difficult to describe quantitatively as a result of the uncertainties in the modeled GPP and the variety of climatic effects from different El Ni?o events (

    Schwalm et al.,2011

    ),therefore the responses of GPP to El Ni?o are only qualitatively shown here.This study only focused on the inter-annual relationships between different climate variables (rainfall,temperature,and downward shortwave radiation) and GPP,so the roles of specific climate changes (e.g.,an increase/decrease in the intensity or frequency of rainfall) were neglected.Given the complexities of the Earth’s climate system,the physical mechanisms in dry episodes associated with El Ni?o and the resulting state of the vegetation in eastern China still require further research.TBMs simulated the physiological processes of crops more poorly than other PFTs that were less affected by human activities (see Appendix A.2).A benefit from the advantages of the BMA method is that this defect in eastern China where crops occupy a large proportion of the area can be minimized by weighting different models according to their performance.Additionally,the intensity of agricultural management (especially irrigation) is difficult to evaluate at a regional level,so its impacts on the conclusions found in this study also need further investigation.

    Fig.9.Spatial patterns of the annual GPP in China averaged from MTE,BMA,and the various MsTMIP model contributors during 1982?2010.

    5.Conclusions

    We selected BMA,a method that can integrate the advantages of a number of different models into the final result,to produce a set of long time-series GPP data in China.When training BMA at ChinaFLUX sites and applying weight coefficients in spatial scale,PFTs were simply classified as trees,grasses,and crops.The validation of BMA GPP proved this kind of up-scaling method is reliable,especially in eastern China.

    In contrast to the tropics,where it is generally believed that local plant physiology can be significantly affected by El Ni?o,this paper focused on the neglected eastern monsoon region of China.The spatial anomalies of BMA GPP averaged from the multiple El Ni?o years in different phases of PDO produced the following hypothesis:El Ni?o during the cool phase of the PDO led to greater GPP in northern China (32°–38°N,111°–122°E) and less GPP in the Yangtze River valley (28°–32°N,111°–122°E);but El Ni?o during the warm phase of PDO reversed the GPP anomalies in these two regions.Even though there were some discrepancies,most individual El Ni?o years supported the hypothesis mentioned above,and therefore the conclusions have high applicability.It should also be noted that the effects of El Ni?o are more significant in northern China than in the Yangtze River valley.

    The synchronized spatiotemporal distribution in selected El Ni?o years and high partial correlation both in strong and weak periods of the EASM all demonstrated that rainfall linked El Ni?o and GPP in eastern China.Although EASM systems are complex (

    Ding and Chan,2005

    ) and the movements of the EASM rain belt are difficult to forecast(

    Gao et al.,2011

    ,

    2014

    ;

    Fan et al.,2012

    ),this study showed that we can still predict how GPP will change in eastern China (28°–38°N,111°–122°E) when El Ni?o occurs if we can identify the phase of the PDO.The results of this study also provide a reference for future research on the impact of El Ni?o on the carbon cycle of terrestrial ecosystems in China.

    .

    This study was jointly supported by the National Key Research and Development Program of China (Grant Nos.2016YFA0602501 and 2018YFA0606004) and the Strategic Priority Research Program of the Chinese Academy of Sciences(Grant Nos.XDA20040301 and XDA20020201).

    Data availability statement

    The MsTMIP outputs (

    https://nacp.ornl.gov/mstmipdata/

    ) and drivers (

    https://daac.ornl.gov/NACP/guides/NACP_MsTMIP_Mod el_Driver.html

    ) can be downloaded freely.The observed GPP from ChinaFLUX can be obtained from www.cnern.org.cn/data/init-DRsearch.The MOD17A2 products are available from

    http://files.ntsg.umt.edu/data/NTSG_Products/MOD17/GeoTIFF/Monthly_M OD17A2/GeoTIFF_0.05degree/

    .The BMA GPP for China are available in Zenodo with the identifier

    http://doi.org/10.5281/zenodo.3899605

    .

    APPENDIX A Supplementary Information for Materials and Methods

    APPENDIX A.1 Observed GPP from ChinaFLUX

    The monthly observed GPP data from eight typical flux measurement sites in the ChinaFLUX network (

    Yu et al.,2006

    ) were used to validate the modeled GPP.The selected eddy covariance tower sites cover all the main ecosystem types in China and have been successfully applied to validate the performance of different models (

    Li et al.,2013

    ;

    Wang et al.,2015

    ;

    Zhang et al.,2016

    ).Detailed information about these sites is given in

    Table A1

    .

    APPENDIX A.2 Uncertainties in the MsTMIP GPP

    The monthly MsTMIP GPP values in the grid cells corresponding to the location of each flux site were selected for comparison with the observed GPP.The statistical properties obtained in this way only had a credible significance when the dominant plant function types (PFTs) in the selected grid cells were the same as the observed PFTs at the flux sites and their proportions were >50% (

    Peng et al.,2015

    ).The land use and land cover change maps used by MsTMIP were prescribed by the SYNergetic land cover MAP (SYNMAP) (

    Jung et al.,2006

    ).

    Fig.A1

    displays the SYNMAP biome type fractions in the grid cells corresponding to the ChinaFLUX sites (see

    Table A1

    ).Even if weaccounted for the number of mixed trees,the proportion of the target PFT in the grid cell corresponding to Dinghushan was still <50% and did not meet the requirements mentioned above.The grid cell containing the Haibei site had similar problems.Evergreen broad leaf trees (EBT) and shrubs(SHB) are widely distributed in China.Unfortunately,Dinghushan and Haibei are the only two sites in China that share observations,respectively,for subtropical EBT and alpine SHB in the time period 2003–2010,and therefore cannot be replaced.

    Table A1. Information for the eight ChinaFLUX sites used in this study.

    Fig.A1. Biome type compositions in the SYNMAP grids corresponding to the ChinaFLUX sites.

    Taylor diagrams were used to indicate how well the modeled results matched the observed values in terms of three statistics:the Pearson correlation coefficient (R),the root-mean-square difference (RMSD),and the standard deviation (STD) (

    Taylor,2001

    ).The STD and RMSD are proportional to the radial distance from the origin to the observed point.R is defined as the cosine of the azimuth angle.In general,the closer a modeled point is to the observed point,the better the model performs.The Taylor diagrams in

    Fig.A2

    show how closely the GPP from the 14 MsTMIP models matched the observations from 2003 to 2010.The GPP simulations had large uncertainties,especially at YuCheng,where the ground cover was cropland.This means that simple equal-weighted averaging cannot correct the bias when most models have a poor performance for specific PFTs.The diversity in the multi-model results in

    Fig.A2

    shows the need to use a multi-model averaging method that can weight each model according to its simulation skill.

    APPENDIX A.3 Merging of PFTs when using BMA

    When use is made of the observed GPP from China-FLUX to train BMA,the dominant PFTs in the qualified grid cells cannot cover all of the PFTs in the SYNMAP legend.To reduce the errors caused by mixed PFTs in the grid cells and the limited number of ChinaFLUX sites,the PFTs at ChinaFLUX sites and used by MsTMIP were categorized as trees,grasses,and crops when BMA was applied.The detailed classifications of the PFTs at ChinaFLUX sites and SYNMAP legend are listed in

    Table A1

    and

    Table A2

    ,respectively.After being trained by the observed GPP at the grouped ChinaFLUX sites and the MsTMIP GPP in the corresponding grid cells during 2003–08,the weight coefficients of trees,grasses,and crops for different MsTMIP models were calculated by BMA (see

    Table A3

    ).

    Table A2. Summary of SYNMAP land cover types used for Bayesian model averaging.

    Table A3. Weight coefficients for 14 different models calculated by Bayesian model averaging.

    Fig.A2.Taylor diagrams for GPP simulations from 14 MsTMIP models at eight ChinaFLUX sites.The red dashed lines are the observed standard deviations at the sites.

    APPENDIX A.4 Nash–Sutcliffe efficiency index

    The Nash–Sutcliffe efficiency index (E;

    Nash and Sutcliffe,1970

    ) was used to evaluate the performance of the GPP from BMA and individual model in MsTMIP.Its mathematical expression is:

    where

    o

    is the ith observed data point,

    o

    is the mean of the observed data,and

    m

    is the ith modeled data point.The Evalues range between ?∞ and 1.Evalues <0 indicate that the mean value of the observed time series is a better predictor than the model.Evalues between 0 and 1 generally indicate that the performance of the model is at an acceptable level.Evalues equal to 1 indicate a perfect fit.

    APPENDIX A.5 MODIS GPP and MTE GPP

    The monthly MODerate-resolution Imaging Spectroradiometer (MODIS) Version MOD17A2 GPP product (

    Zhao et al.,2005

    ) was used to evaluate the temporal evolution of the BMA GPP in the target zones.To reduce the disadvantages caused by the difference in spatial resolution (MODIS:0.05° × 0.05°;MsTMIP:0.5° × 0.5°) and land cover classification system (MODIS:MOD12Q1;MsTMIP:SYNMAP),we only compared their spatially averaged anomalies from 2000 to 2010.The monthly Model Tree Ensemble (MTE) GPPupscaled from FLUXNET observations to global 0.5° × 0.5°grids during 1982?2010 (

    Jung et al.,2011

    ) was selected to estimate the magnitude of BMA GPP.Considering that the empirical relationships in machine learning algorithm are time-sensitive,we weakened the time dimension and only compared the spatial patterns of multi-year averaged GPP from MTE,BMA,and MsTMIP to show the rationality of BMA GPP distribution.

    又爽又黄无遮挡网站| 日本三级黄在线观看| av线在线观看网站| 赤兔流量卡办理| 日韩欧美国产在线观看| kizo精华| 我的老师免费观看完整版| 日韩中字成人| 免费观看人在逋| 午夜精品在线福利| 中文字幕人妻熟人妻熟丝袜美| 欧美潮喷喷水| 99视频精品全部免费 在线| 日本免费一区二区三区高清不卡| 色噜噜av男人的天堂激情| 九色成人免费人妻av| 在线观看美女被高潮喷水网站| 69av精品久久久久久| 国产老妇伦熟女老妇高清| 精品99又大又爽又粗少妇毛片| 成人午夜精彩视频在线观看| 波多野结衣巨乳人妻| 成年av动漫网址| 成人一区二区视频在线观看| 久久久久久久久久久免费av| 亚洲精品自拍成人| 长腿黑丝高跟| 大香蕉97超碰在线| 国产乱人偷精品视频| 亚洲成av人片在线播放无| 免费av不卡在线播放| 2022亚洲国产成人精品| 97在线视频观看| 一区二区三区免费毛片| 中文字幕免费在线视频6| 国产伦精品一区二区三区四那| 欧美最新免费一区二区三区| 国产三级中文精品| 高清av免费在线| 毛片一级片免费看久久久久| 极品教师在线视频| 国产午夜精品论理片| 日本-黄色视频高清免费观看| av免费观看日本| 少妇的逼好多水| 欧美一区二区亚洲| 秋霞在线观看毛片| 99久久精品一区二区三区| av福利片在线观看| 久久这里只有精品中国| 一夜夜www| 51国产日韩欧美| 久久久久久久久久久丰满| 91狼人影院| 网址你懂的国产日韩在线| 国产精品国产三级专区第一集| 七月丁香在线播放| 一边摸一边抽搐一进一小说| 欧美97在线视频| 国产免费一级a男人的天堂| 国产精品国产三级国产av玫瑰| 最后的刺客免费高清国语| 99久久成人亚洲精品观看| av福利片在线观看| 亚洲精品色激情综合| 国产亚洲av嫩草精品影院| 国产 一区 欧美 日韩| 最近中文字幕2019免费版| 国产精品久久久久久精品电影| 嫩草影院入口| 亚洲人成网站高清观看| 亚洲激情五月婷婷啪啪| 日本爱情动作片www.在线观看| 亚洲人与动物交配视频| 看片在线看免费视频| 白带黄色成豆腐渣| 五月伊人婷婷丁香| av视频在线观看入口| 建设人人有责人人尽责人人享有的 | 人人妻人人澡欧美一区二区| 久久久精品94久久精品| 午夜激情福利司机影院| 波多野结衣巨乳人妻| 日本-黄色视频高清免费观看| 成人国产麻豆网| 91久久精品国产一区二区三区| 美女cb高潮喷水在线观看| 中文字幕免费在线视频6| 欧美激情国产日韩精品一区| 婷婷色综合大香蕉| 国产精品女同一区二区软件| 最近视频中文字幕2019在线8| 亚洲成色77777| 免费人成在线观看视频色| 国产在线一区二区三区精 | 国产黄色小视频在线观看| 又粗又硬又长又爽又黄的视频| 成人综合一区亚洲| 99久久人妻综合| 99热这里只有精品一区| 久久精品国产鲁丝片午夜精品| videos熟女内射| 特大巨黑吊av在线直播| 三级毛片av免费| 51国产日韩欧美| 国产精品伦人一区二区| 亚洲图色成人| 国产精品麻豆人妻色哟哟久久 | 亚洲av免费在线观看| 69人妻影院| 国产亚洲一区二区精品| 不卡视频在线观看欧美| 亚洲欧美成人综合另类久久久 | 成人欧美大片| 午夜福利高清视频| 长腿黑丝高跟| 2021少妇久久久久久久久久久| 日本三级黄在线观看| 日韩国内少妇激情av| 亚洲精品久久久久久婷婷小说 | 中文字幕精品亚洲无线码一区| 国产综合懂色| 狂野欧美激情性xxxx在线观看| 久久99热这里只频精品6学生 | 看片在线看免费视频| 身体一侧抽搐| 成年女人永久免费观看视频| 免费黄色在线免费观看| 赤兔流量卡办理| 中文字幕av成人在线电影| 国产精品国产三级国产av玫瑰| 好男人视频免费观看在线| 亚洲不卡免费看| 国产淫语在线视频| 国产综合懂色| 搡老妇女老女人老熟妇| 久久6这里有精品| 少妇人妻精品综合一区二区| 色尼玛亚洲综合影院| 国产高清三级在线| 蜜桃久久精品国产亚洲av| 男女那种视频在线观看| 免费观看精品视频网站| 中国美白少妇内射xxxbb| 欧美不卡视频在线免费观看| 插逼视频在线观看| 中文天堂在线官网| 哪个播放器可以免费观看大片| 精品国产三级普通话版| 国产亚洲91精品色在线| 岛国在线免费视频观看| 桃色一区二区三区在线观看| 狠狠狠狠99中文字幕| 国产毛片a区久久久久| 国产色爽女视频免费观看| 高清毛片免费看| 成人漫画全彩无遮挡| 国产视频首页在线观看| 欧美精品国产亚洲| a级毛片免费高清观看在线播放| 国产成年人精品一区二区| 精品一区二区免费观看| 久久婷婷人人爽人人干人人爱| 最近最新中文字幕免费大全7| 日本猛色少妇xxxxx猛交久久| 丰满乱子伦码专区| 国产毛片a区久久久久| 纵有疾风起免费观看全集完整版 | 91av网一区二区| 亚洲av男天堂| 亚洲最大成人av| 久久久久网色| 国产精品久久电影中文字幕| 亚洲,欧美,日韩| 精品一区二区三区人妻视频| 免费无遮挡裸体视频| 中文在线观看免费www的网站| 看片在线看免费视频| 女人久久www免费人成看片 | 欧美性感艳星| 亚洲欧美一区二区三区国产| 国产 一区 欧美 日韩| 成人高潮视频无遮挡免费网站| 亚洲精品自拍成人| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 中文字幕制服av| 国产午夜精品论理片| 亚洲中文字幕一区二区三区有码在线看| 国内揄拍国产精品人妻在线| 国产精品女同一区二区软件| 亚洲国产精品成人久久小说| 午夜福利成人在线免费观看| 人妻少妇偷人精品九色| 麻豆精品久久久久久蜜桃| 别揉我奶头 嗯啊视频| 日本三级黄在线观看| 少妇裸体淫交视频免费看高清| 国产一区亚洲一区在线观看| 色噜噜av男人的天堂激情| 汤姆久久久久久久影院中文字幕 | 最近2019中文字幕mv第一页| a级毛色黄片| 成人亚洲精品av一区二区| 久久精品国产99精品国产亚洲性色| 久久久久久久久久成人| 精品人妻视频免费看| 精品久久国产蜜桃| 在现免费观看毛片| 欧美丝袜亚洲另类| 久99久视频精品免费| 亚洲av免费高清在线观看| 插阴视频在线观看视频| 国内精品宾馆在线| 久久久久网色| 建设人人有责人人尽责人人享有的 | 18禁裸乳无遮挡免费网站照片| 一级毛片久久久久久久久女| 日韩高清综合在线| 男人的好看免费观看在线视频| 亚洲第一区二区三区不卡| 国产成人a区在线观看| 亚洲在线自拍视频| 日本熟妇午夜| 亚洲性久久影院| 欧美最新免费一区二区三区| 久久热精品热| 2022亚洲国产成人精品| 国产精品1区2区在线观看.| 成人美女网站在线观看视频| 亚洲在线观看片| 欧美变态另类bdsm刘玥| 国产亚洲精品av在线| 国产精品99久久久久久久久| 小说图片视频综合网站| 亚洲欧美精品自产自拍| av播播在线观看一区| 色综合亚洲欧美另类图片| 国产欧美另类精品又又久久亚洲欧美| АⅤ资源中文在线天堂| av视频在线观看入口| 一级爰片在线观看| 少妇丰满av| 亚洲国产精品成人久久小说| av天堂中文字幕网| 日本免费a在线| 丝袜喷水一区| 男人舔女人下体高潮全视频| 99热这里只有是精品在线观看| 插逼视频在线观看| 午夜亚洲福利在线播放| av专区在线播放| 男人舔女人下体高潮全视频| 韩国高清视频一区二区三区| 久久精品91蜜桃| 美女xxoo啪啪120秒动态图| av黄色大香蕉| av国产免费在线观看| 2021少妇久久久久久久久久久| 国产精品不卡视频一区二区| 亚洲熟妇中文字幕五十中出| 久久精品国产亚洲av天美| 边亲边吃奶的免费视频| 成人午夜高清在线视频| 日韩在线高清观看一区二区三区| 国产国拍精品亚洲av在线观看| 干丝袜人妻中文字幕| 日本-黄色视频高清免费观看| 观看美女的网站| 成人综合一区亚洲| 亚洲av男天堂| 内地一区二区视频在线| 哪个播放器可以免费观看大片| 婷婷色综合大香蕉| 乱系列少妇在线播放| 夫妻性生交免费视频一级片| 久久精品久久精品一区二区三区| 日日摸夜夜添夜夜添av毛片| 国产精品三级大全| 亚洲国产精品国产精品| 久久精品久久久久久噜噜老黄 | 国产极品精品免费视频能看的| 日本爱情动作片www.在线观看| 成人欧美大片| 亚洲精品乱久久久久久| 国产在线一区二区三区精 | 亚洲国产色片| 九九在线视频观看精品| 亚洲五月天丁香| av在线亚洲专区| 国产国拍精品亚洲av在线观看| 亚洲图色成人| 午夜a级毛片| 久久99热6这里只有精品| 我的老师免费观看完整版| 欧美zozozo另类| kizo精华| 国产亚洲最大av| 久久人妻av系列| 看非洲黑人一级黄片| 精品欧美国产一区二区三| 国产伦一二天堂av在线观看| 舔av片在线| 成人综合一区亚洲| 国产高清视频在线观看网站| 真实男女啪啪啪动态图| 七月丁香在线播放| 高清在线视频一区二区三区 | 51国产日韩欧美| 欧美区成人在线视频| 岛国在线免费视频观看| 久久精品国产亚洲av涩爱| 有码 亚洲区| av在线老鸭窝| 国产私拍福利视频在线观看| 亚洲不卡免费看| 久久99蜜桃精品久久| 日日撸夜夜添| 久久精品影院6| 欧美一区二区国产精品久久精品| 波野结衣二区三区在线| 日韩在线高清观看一区二区三区| 中国美白少妇内射xxxbb| 国产精品综合久久久久久久免费| 深夜a级毛片| 美女大奶头视频| 亚洲av成人精品一二三区| 亚洲最大成人av| 国产精品熟女久久久久浪| 一级av片app| 国产探花在线观看一区二区| 一个人看的www免费观看视频| 神马国产精品三级电影在线观看| 亚洲18禁久久av| 午夜老司机福利剧场| 久久久久久久国产电影| 国产成人一区二区在线| 日本黄大片高清| 男女边吃奶边做爰视频| 秋霞伦理黄片| 久久久久久久久大av| 亚洲国产最新在线播放| 免费在线观看成人毛片| 联通29元200g的流量卡| 国产成人一区二区在线| 国产午夜精品一二区理论片| 久热久热在线精品观看| 99热全是精品| videos熟女内射| 久久这里有精品视频免费| 亚洲经典国产精华液单| 日韩精品青青久久久久久| 国产一区二区亚洲精品在线观看| 校园人妻丝袜中文字幕| 2021少妇久久久久久久久久久| 免费黄色在线免费观看| 又黄又爽又刺激的免费视频.| 伦理电影大哥的女人| 久久久亚洲精品成人影院| 在线免费观看的www视频| 国语对白做爰xxxⅹ性视频网站| 欧美高清性xxxxhd video| 国产成人一区二区在线| 国产极品天堂在线| a级毛片免费高清观看在线播放| 国产精品久久久久久精品电影小说 | 国产v大片淫在线免费观看| 日本与韩国留学比较| 秋霞在线观看毛片| 国产极品天堂在线| 女人被狂操c到高潮| 只有这里有精品99| 青春草国产在线视频| 日韩精品青青久久久久久| 水蜜桃什么品种好| 久久久午夜欧美精品| 中文字幕精品亚洲无线码一区| 纵有疾风起免费观看全集完整版 | 国产69精品久久久久777片| 日产精品乱码卡一卡2卡三| 最近的中文字幕免费完整| 日韩一本色道免费dvd| 国产精品美女特级片免费视频播放器| 午夜日本视频在线| 美女高潮的动态| 国产精品国产三级专区第一集| 99热全是精品| 亚洲国产精品sss在线观看| 久久精品国产99精品国产亚洲性色| av播播在线观看一区| 干丝袜人妻中文字幕| 国产真实伦视频高清在线观看| a级毛色黄片| 人人妻人人澡欧美一区二区| 嘟嘟电影网在线观看| 欧美日韩精品成人综合77777| 18禁动态无遮挡网站| 少妇人妻一区二区三区视频| 晚上一个人看的免费电影| 少妇被粗大猛烈的视频| 国产av一区在线观看免费| 七月丁香在线播放| 欧美97在线视频| 国产国拍精品亚洲av在线观看| 国产淫语在线视频| 国产高清国产精品国产三级 | 嫩草影院入口| 波野结衣二区三区在线| 国产亚洲精品久久久com| 身体一侧抽搐| 亚洲欧美一区二区三区国产| 国产三级在线视频| 亚洲欧美精品综合久久99| 26uuu在线亚洲综合色| 欧美高清性xxxxhd video| 精品一区二区免费观看| 国产成人a区在线观看| 少妇的逼好多水| 日韩一区二区视频免费看| 国产精品日韩av在线免费观看| 国产精品人妻久久久久久| 18禁在线播放成人免费| 99热网站在线观看| 国产免费又黄又爽又色| 有码 亚洲区| 欧美高清成人免费视频www| 深爱激情五月婷婷| 黄色配什么色好看| 亚洲综合精品二区| 两个人的视频大全免费| 国产精品国产三级国产专区5o | 69人妻影院| 亚洲国产欧美在线一区| 国产黄片视频在线免费观看| 国产乱人偷精品视频| 日本-黄色视频高清免费观看| 麻豆国产97在线/欧美| 大香蕉97超碰在线| 亚洲在线自拍视频| 男的添女的下面高潮视频| 国产av不卡久久| 成人二区视频| 色网站视频免费| 色综合站精品国产| 午夜精品国产一区二区电影 | 人人妻人人澡人人爽人人夜夜 | 欧美日韩综合久久久久久| 只有这里有精品99| 国产黄a三级三级三级人| 国产老妇伦熟女老妇高清| 一级黄片播放器| 中文字幕人妻熟人妻熟丝袜美| 精品熟女少妇av免费看| 五月伊人婷婷丁香| 两个人视频免费观看高清| 在线播放国产精品三级| 我的老师免费观看完整版| videossex国产| 大香蕉久久网| 青春草亚洲视频在线观看| 久久精品熟女亚洲av麻豆精品 | 国产乱人视频| 日本欧美国产在线视频| 久久久精品欧美日韩精品| 免费av观看视频| 国产高清有码在线观看视频| 日本黄色视频三级网站网址| 中文字幕精品亚洲无线码一区| 免费观看a级毛片全部| 久久久久久久久久久丰满| 国产国拍精品亚洲av在线观看| 中国美白少妇内射xxxbb| 色播亚洲综合网| 全区人妻精品视频| 韩国av在线不卡| 日韩精品有码人妻一区| 亚洲人成网站高清观看| 国产视频首页在线观看| 网址你懂的国产日韩在线| 男人和女人高潮做爰伦理| 亚洲精品自拍成人| 欧美成人午夜免费资源| 九色成人免费人妻av| 色5月婷婷丁香| 亚洲婷婷狠狠爱综合网| 天天躁夜夜躁狠狠久久av| 国产精品久久久久久久久免| 亚州av有码| 日韩在线高清观看一区二区三区| 久久精品夜色国产| 欧美一区二区精品小视频在线| 1024手机看黄色片| 国产免费一级a男人的天堂| 成人性生交大片免费视频hd| 日韩欧美 国产精品| 国产真实乱freesex| 熟女电影av网| 久久久久九九精品影院| 你懂的网址亚洲精品在线观看 | av天堂中文字幕网| 国模一区二区三区四区视频| 大香蕉久久网| 美女内射精品一级片tv| 一级黄色大片毛片| 亚洲欧美日韩无卡精品| av免费在线看不卡| 成年av动漫网址| 2021天堂中文幕一二区在线观| 中文字幕精品亚洲无线码一区| 午夜日本视频在线| 搡老妇女老女人老熟妇| 国产精品不卡视频一区二区| 久久久精品94久久精品| 亚洲三级黄色毛片| 99在线人妻在线中文字幕| 啦啦啦韩国在线观看视频| 国产三级中文精品| 日日撸夜夜添| 成年av动漫网址| 日韩制服骚丝袜av| 午夜福利成人在线免费观看| 国产午夜精品久久久久久一区二区三区| 非洲黑人性xxxx精品又粗又长| 天天躁夜夜躁狠狠久久av| 能在线免费观看的黄片| 国产极品精品免费视频能看的| 国产在视频线精品| av.在线天堂| 18禁裸乳无遮挡免费网站照片| 精品国产一区二区三区久久久樱花 | 我的老师免费观看完整版| 成人av在线播放网站| 久久久久久九九精品二区国产| 99久国产av精品国产电影| 高清午夜精品一区二区三区| 亚洲av免费高清在线观看| 国产高潮美女av| 熟女人妻精品中文字幕| 婷婷色麻豆天堂久久 | 久久久午夜欧美精品| 久久久久性生活片| 国产三级在线视频| 欧美成人精品欧美一级黄| 精品国产三级普通话版| av黄色大香蕉| 亚洲18禁久久av| 可以在线观看毛片的网站| 舔av片在线| 亚洲在久久综合| 免费观看性生交大片5| 亚洲国产最新在线播放| 亚洲在久久综合| 听说在线观看完整版免费高清| 欧美日本视频| 欧美性感艳星| 久久久久免费精品人妻一区二区| 综合色丁香网| 亚洲,欧美,日韩| 99国产精品一区二区蜜桃av| 国产黄片美女视频| 亚洲精品国产av成人精品| 大话2 男鬼变身卡| 三级毛片av免费| 人人妻人人澡欧美一区二区| 特大巨黑吊av在线直播| 久久婷婷人人爽人人干人人爱| 国产大屁股一区二区在线视频| 亚洲av免费在线观看| 亚洲丝袜综合中文字幕| 亚洲精品国产成人久久av| 亚洲性久久影院| 校园人妻丝袜中文字幕| 亚洲成色77777| 久久国内精品自在自线图片| 在线天堂最新版资源| 免费观看的影片在线观看| 国产免费一级a男人的天堂| 精品无人区乱码1区二区| 欧美日韩在线观看h| 国产91av在线免费观看| 热99在线观看视频| 久久久a久久爽久久v久久| 一级毛片电影观看 | ponron亚洲| 亚洲图色成人| 男女边吃奶边做爰视频| 亚洲色图av天堂| 成人性生交大片免费视频hd| 久久久成人免费电影| 日韩一区二区视频免费看| 亚洲国产最新在线播放| 成人三级黄色视频| 成人欧美大片| 一个人观看的视频www高清免费观看| 久久久午夜欧美精品| 国产成人a区在线观看| 一二三四中文在线观看免费高清| 久久久午夜欧美精品| 欧美成人午夜免费资源| 99九九线精品视频在线观看视频| av在线蜜桃| 久久久久久久久久久丰满| 视频中文字幕在线观看| 男的添女的下面高潮视频| 中文字幕人妻熟人妻熟丝袜美| 久久99热这里只频精品6学生 | 午夜精品在线福利| 在线观看一区二区三区| 日日摸夜夜添夜夜添av毛片| 欧美色视频一区免费| 亚洲伊人久久精品综合 | 狂野欧美激情性xxxx在线观看| 熟女电影av网| 午夜福利成人在线免费观看| 青春草亚洲视频在线观看| 两个人的视频大全免费|