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

    Improvement of Soil Moisture Simulation in Eurasia by the Beijing Climate Center Climate System Model from CMIP5 to CMIP6

    2021-02-26 08:22:22YinghanSANGHongLiRENXueliSHIXiaofengXUandHaishanCHEN
    Advances in Atmospheric Sciences 2021年2期

    Yinghan SANG, Hong-Li REN, Xueli SHI, Xiaofeng XU, and Haishan CHEN

    1School of Atmospheric Sciences, Nanjing University of Information Science &Technology, Nanjing 210044, China

    2State Key Laboratory of Severe Weather, Chinese Academy of Meteorological Sciences, Beijing 100081, China

    3Laboratory for Climate Studies &CMA–NJU Joint Laboratory for Climate Prediction Studies,National Climate Center, China Meteorological Administration, Beijing 100081, China

    4China Meteorological Administration, Beijing 100081, China

    ABSTRACT

    Key words:BCC-CSM,soil moisture,CMIP6,historical simulation,Eurasia

    1.Introduction

    Soil moisture (SM) plays a vital role in affecting climatic variability through influencing many physical processes in land–atmosphere interactions. By regulating surface water, evapotranspiration, and latent heat, as well as ground fluxes, SM variations feed back to the near-surface climate, extending to the boundary layer, impact vertical stability, and further affect precipitation (Seneviratne et al., 2010;Berg and Sheffield, 2018; Ruosteenoja et al., 2018). One salient character of SM is its long memory (Koster and Suarez,2001; Shinoda, 2001; Ruosteenoja et al., 2018), which can last 30–90 days for some climatic elements (Dirmeyer et al.,2009). The climatic response to the variation in SM changes with regions and seasons (Yeh et al., 1984; Koster et al.,2003; Koster et al., 2004). Especially in the Eurasian continent with its complex and diverse climate, SM not only affects ground fluxes and climate in the monsoon zone, but also impacts the evaporation and water vapor transport and further the precipitation and temperature inland (Ma et al.,2000; Zuo and Zhang, 2007; Dai and Zuo, 2010; Zhang et al., 2016a).

    Owing to the deficiency of reliable long-term and large-spatial-scale in-situ observational datasets, alternative data sources have been made available for researching SM variation. Represented by the Global Land Data Assimilation System (GLDAS), several assimilated SM data products have been developed. Research has found that data assimilation can significantly improve the accuracy of middle and deep SM estimates and the assimilated data products are reliable for describing the SM annual cycle and short-term variability (Renzullo et al., 2014). Simulations from offline land surface models and coupled general circulation models (GCMs) are another tool for investigating the continuity of SM spatially and temporally (Srinivasan et al.,2000; Koster et al., 2009). With the gradual improvements of model systems in the past several decades, their ability to reproduce SM has led to great progress, but there are still deficiencies. Guo and Dirmeyer (2006) pointed out that,although models can reproduce SM anomalies to a certain extent, they do not simulate absolute soil water content accurately. Xia et al. (2015) also showed that model simulations may overestimate the SM in the Northern Hemisphere.Cheng et al. (2015) analyzed the simulations of 20 CMIP5(phase 5 of the Coupled Model Intercomparison Project)GCMs and revealed that there has been an obvious decreasing trend in annual-mean near-surface SM over eastern Asia. Research also demonstrated that the majority of CMIP5 GCMs likewise simulate soil drying in the Northern Hemisphere, for nearly the whole continent of North America in summer and everywhere apart from the Arctic regions in spring (Dirmeyer et al., 2013). Thus, although there is a certain biases, simulations based on climate models are efficient tools to investigate the variations and effects of SM in the climate system.

    However, because so many models have been developed with different descriptions of land surface and SM-related schemes, under the same simulation framework,say the CMIP experiments, comparisons of multi-model simulations will provide an overall assessment of the model performance, but are inconvenient for investigating the reasons. On the contrary, comparisons of different model versions of a specified model system in the different phases of CMIP will provide a good opportunity to understand the SM simulation performance.

    The Beijing Climate Center Climate System Model(BCC-CSM) is one of the fully coupled climate system models that has participated in both the CMIP5 and CMIP6 experiments. Most of the CMIP6 simulations have now been finished, with datasets available via the ESGF website (Xin et al., 2019). Wu et al. (2019) summarized the general model improvements and performances of different components of models from CMIP5 to CMIP6. However, the SM has not yet been systematically evaluated. Therefore, in this paper,we assess the SM simulated by BCC-CSM through comparisons with assimilated data and observations as well as results of the different BCC model versions. Two medium resolutions are selected for comparison here, i.e., BCC-CSM2-MR for CMIP6 and BCC-CSM1.1m for CMIP5, in the historical experiments. Then, a horizontal comparison is made between BCC-CSM2-MR and 13 CMIP6 models. Finally,the coupling between precipitation and SM in model simulations is also discussed.

    2.Data and methods

    2.1.Data

    2.1.1.BCC-CSM historical simulations in CMIP5 and CMIP6

    Recently, BCC-CSM has been upgraded to its second generation (Wu et al., 2020). The medium-resolution version (BCC-CSM2-MR) has carried out most of the CMIP6 DECK and MIP experiments (Eyring et al., 2016; Xin et al.,2019). Its previous generation (BCC-CSM1.1m; Wu et al.,2013, 2014) participated in CMIP5 (Taylor et al., 2012; Xin et al., 2012). Among the experiments endorsed by CMIP5 and CMIP6, historical simulation is one of the entry cards for models to participate in the project. The historical period is defined as beginning in 1850 and extends to near to the present day (2012 for CMIP5 and 2014 for CMIP6). The CMIP historical simulation provides a good opportunity to assess model ability in simulating climatic variability and trends. In this study, monthly surface (0–10 cm) SM data from the historical simulation of BCC-CSM2-MR are evaluated by comparing with 13 other CMIP6 models (Table 1)and BCC-CSM1.1m in CMIP5. The monthly precipitation data of the BCC-CSMs are also used, to analyze the coupling between SM and precipitation.

    2.1.2.Reference data

    Representative of assimilated SM data products,GLDAS, which comprises ingested satellite- and groundbased observational data products and uses advanced land surface modeling and data assimilation techniques (Rodell et al., 2004), is widely used in studies of land–atmosphere interaction (Cheng et al., 2013; Zhang et al., 2016b). Therefore,we select the monthly data of near-surface (0–10 cm) SM from GLDAS v2.0 and v2.1 as the reference data in this study. In addition, we also use a station-based observational SM dataset in China produced by Wang and Shi(2019), to make the best use of the observations available at present. The SM measurements at a total of 1471 stations for the period January 1992 to September 2013 are from theNational Meteorological Information Center of the China Meteorological Administration. Among these, 732 stations with good spatial and temporal continuity from January 1992 to December 2012 are ultimately selected following quality control processes. We use the data of the first layer (0–10 cm)and transform the units into kg mfor comparison. In the analysis of the coupling between SM and precipitation, we use the monthly data of precipitation from the CPC (Climate Prediction Center of NOAA) Merged Analysis of Precipitation(CMAP; Huffman et al., 1997).

    Table 1. Details of the 14 CMIP6 models used in this study.

    2.2.Research region and diagnostic method

    The study region in this paper is Eurasia (10°–80°N,30°–160°E) and the study period ranges from 1979 to 2012,in which we use the annual mean and the four seasons defined as December–January–February (DJF) for winter,March–April–May (MAM) for spring, June–July–August(JJA) for summer, and September–October–November(SON) for autumn. For convenience of comparison, both the simulation data and reference data are regrided to a uniform resolution of 1° × 1° using bilinear interpolation.Although the dataset of the 732 stations is quality controlled, there are still plenty of missing values in the spatial and temporal range. Therefore, it is inaccurate and inconvenient to interpolate the station dataset to the grid. Accordingly, when using the station dataset from 1992 to 2012 as the reference to make comparisons, we interpolate the simulation data in China to the positions of the 732 stations. Pearson correlation (r), root-mean-square error (RMSE), standard deviation and linear regression, the most commonly used metrics (Legates and McCabe, 1999), are used to quantify the agreement between reference data and model simulations. Besides, empirical orthogonal function (EOF) analysis and spectral analysis are performed to identify the major modes of the surface SM anomalies in Eurasia and its periodic characteristics.

    3.Results

    3.1.Comparison between BCC-CSM2-MR and BCCCSM1.1m

    3.1.1.Spatial distribution in Eurasia The climatological mean surface SM in terms of annual mean and the four seasons of GLDAS and BCC-CSM simulations are firstly presented (Fig.1). Basically, except the less pronounced gradient in model simulations, both BCCCSM2-MR and BCC-CSM1.1m are able to capture the spatial patterns of the SM of GLDAS, but BCC-CSM2-MR is closer to GLDAS in terms of the higher pattern correlation.GLDAS shows a triple pattern in the annual mean, whereby the surface soil with higher soil water content is typically located in central and western Siberia and South China, while that with lower soil water content is mainly distributed in central Asia and the Arabian Peninsula (Fig.1a). From northwest to southeast, the SM in Eurasia exhibits a wet–dry–wet distribution, which is consistent with previous research(Nijssen et al., 2001; Guo et al., 2003; Erdenebat and Sato,2018). This pattern is also evident in the BCC-CSM2-MR and BCC-CSM1.1m model simulations (Figs. 1b and 1c).The pattern correlation coefficient between BCC-CSM2-MR and GLDAS is 0.97 (at the 0.01 significance level),which is higher than that between BCC-CSM1.1m and GLDAS (r = 0.93). In different seasons, the pattern correlation coefficients of BCC-CSM2-MR and GLDAS range from 0.96 to 0.97, which are higher than those of BCCCSM1.1m (between 0.92 and 0.94). The average pattern correlation coefficient between BCC-CSM2-MR and GLDAS is 0.97, which is a rise of 4.3% over the previous version.Therefore, BCC-CSM2-MR is more capable than BCCCSM1.1m in its reproduction of the distribution of the climate mean surface SM.

    Fig.1. Climate mean surface soil moisture (units: kg m?2) over Eurasia (1979–2012) on annual and seasonal time scales: (a,d, g, j, m) GLDAS as the reference data; (b, e, h, k, n) BCC-CSM2-MR and the pattern correlation coefficients with GLDAS; (c, f, i, l, o) BCC-CSM1.1m and the pattern correlation coefficients with GLDAS.

    The difference in patterns of the climate mean SM between the BCC model simulations and GLDAS in terms of their annual mean and the different seasons are shown in Fig.2. As the first two columns of the figure show, BCCCSM2-MR and BCC-CSM1.1m both tend to reproduce soil that is drier than in GLDAS over the whole of Eurasia, with a large coverage of negative values, especially in southern China. Underestimations of SM in climate models have been demonstrated in many previous studies. For example,Ramillien et al. (2003) indicated that Land Dynamics hydrological model tends to underestimate the absolute water storage in the soil and provide smoother values than in-situ measurements. However, in several regions (Siberia, Northeast China, the Yangtze–Huaihe River basin, etc.), the results of the BCC-CSMs tend to be significantly wetter than those of GLDAS. Comparing Figs. 2a and b, we can see that the areas and absolute values of most regions with large differences have reduced from BCC-CSM1.1m to BCC-CSM2-MR in terms of their annual mean. Taking the positive value in Siberia as an example, the SM difference between BCCCSM1.1m and GLDAS is greater than 20 kg m. However,it decreases to around 5 kg min BCC-CSM2-MR, with obvious improvement. This is also the case in the four seasons. Furthermore, compared to the previous-generation model, the RMSE values of BCC-CSM2-MR are lower both in terms of their annual mean and in the four seasons. The mean RMSE of the climate mean state declines by 7.5%from 9.37 (BCC-CSM1.1m) to 8.67 (BCC-CSM2-MR),which demonstrates that BCC-CSM2-MR is more capable of describing the actual distribution of the climate mean surface SM, with fewer biases. The differences between BCCCSM2-MR and BCC-CSM1.1m are shown in the rightmost panel of Fig.2. Clear improvements can be seen in Siberia,Northeast China, and the Yangtze–Huaihe River basin,which are basically the same as the areas of large BCC model–GLDAS differences mentioned above.

    Fig.2. Differences in climate mean surface soil moisture (units: kg m?2) between GLDAS and BCC simulations on annual and seasonal time scales: (a, d, g, j, m) differences and RMSEs between GLDAS and BCC-CSM2-MR; (b, e, h, k, n)differences and RMSEs between GLDAS and BCC-CSM1.1m; (c, f, i, l, o) differences between BCC-CSM2-MR and BCCCSM1.1m.

    Figure 3 compares the standard deviations of the climate mean surface SM over Eurasia in terms of their annual mean and in the different seasons between the BCC model simulations and the GLDAS data. GLDAS shows that the larger standard deviations tend to be found at high latitudes and in part of central Asia, and clearly during winter and spring (Figs. 3d and g). This illustrates that the surface SM over these regions varies greatly and is spread out over a wider range in winter and spring. Conversely, the variations of surface SM in the low–middle latitudes are quite small, as illustrated by the relatively lower standard deviations over these areas. Gu et al. (2019) indicated that significantly lower SM is generally found in Russia and northeastern Asia, which are similar to the areas with large variations in Figure 3. It appears that both BCC-CSM2-MR and BCC-CSM1.1m are able to capture the spatial patterns,except that the amplitudes of variation are relatively weaker than in GLDAS. For instance, the standard deviation over western Siberia during winter is supposed to exceed 6 kg m(Fig.3d), but the simulated values in the BCC models range between 4 kg mand 5 kg m. Despite the underestimations to a certain extent, BCC-CSM2-MR has made progress in terms of the standard deviation distributions, with higher pattern correlation coefficients than BCC-CSM1.1m in the annual mean and most seasons. The average pattern correlation coefficient increases by 4%, from 0.81 to 0.84. Therefore, BCC-CSM2-MR is better at describing the variations of surface SM.

    Fig.3. Standard deviations of climate mean surface soil moisture on annual and seasonal time scale: (a, d, g, j, m) GLDAS as the reference data; (b, e, h, k, n) BCC-CSM2-MR and the pattern correlation coefficients of standard deviations with GLDAS; (c, f, i, l, o) BCC-CSM1.1m and the pattern correlation coefficients of standard deviations with GLDAS.

    We also calculate the differences of the aforementioned standard deviations over Eurasia. Although the difference patterns between BCC-CSM2-MR and GLDAS are similar to those between BCC-CSM1.1m and GLDAS in terms of their annual mean (Figs. 4a and b), the discrepancies of BCC-CSM2-MR are less significant than those of BCCCSM1.1m in the four seasons, especially in western Siberia and central Asia during spring, summer and winter, because the areas and absolute negative values decrease. That is to say, there is a certain underestimation in BCC-CSM1.1m when describing the standard deviations of SM. In addition,BCC-CSM2-MR is more skillful than BCC-CSM1.1m in terms of RMSE. The RMSE values of BCC-CSM2-MR are distinctly lower than those of BCC-CSM1.1m in terms of their annual mean and all seasons. The average RMSE reduces by 15.6%, from 1.56 to 1.35. This demonstrates that BCC-CSM2-MR is more capable of describing the spatial distributions of surface SM standard deviations, with fewer biases. From the rightmost panel of Fig.4, we can see that the reason for the better performances in BCC-CSM2-MR is that the standard deviations described in BCC-CSM2-MR are systematically greater than those in BCC-CSM1.1m.This overcomes the defect of the underestimation in the previous generation, hence allowing BCC-CSM2-MR to present the degree of surface SM variation more accurately, especially in the middle and high latitudes.

    Fig.4. Differences in standard deviations of mean surface soil moisture between GLDAS and model simulations on annual and seasonal times scales: (a, d, g, j, m) differences and RMSEs between GLDAS and BCC-CSM2-MR; (b, e, h, k, n)differences and RMSEs between GLDAS and BCC-CSM1.1m; (c, f, i, l, o) differences between BCC-CSM2-MR and BCCCSM1.1m.

    To make the results of the comparison more robust, we also employ a station-based observational SM dataset in China as the reference data. Due to the relatively greater number of missing values in the cold seasons, Fig.5 only shows the difference patterns of the climate mean SM and the distributions of the standard deviations in summer; the results of the annual mean and other seasons are listed in Tables 2 and 3. From Figs. 5a and b, we can see that the SM over the Yangtze–Huaihe River basin and Northeast China simulated by BCC-CSM1.1m is obviously wetter than observed,but this situation improves in BCC-CSM2-MR. Figure 5c shows that through revising the overestimations over the regions mentioned above in the previous generation, the biases of the climate mean SM relative to observations are less severe in BCC-CSM2-MR. The reduction in RMSE(from 10.92 to 9.62) also demonstrates that the ability of the model simulation has been improved. It is worth noting that the patterns shown in Figs. 5a–c are consistent with the patterns in Fig.2, where the reference data are from GLDAS.Figures 5d and e show that both BCC-CSM1.1m and BCCCSM2-MR can capture the pattern of variation in SM over North China as being larger than in other areas, but the amplitudes are smaller than those of the observations to different degrees. However, BCC-CSM2-MR (standard deviation: 2.04) still performs better than BCC-CSM2-MR (standard deviation: 1.63), based on the closer average standard deviation to that of the observations (standard deviation:3.22). The simulations of BCC-CSM2-MR in terms of the annual mean and the other seasons are also improved(Tables 2 and 3).

    3.1.2.EOF analysis

    Fig.5. Differences in climate mean surface soil moisture between observations and BCC simulations over China in summer and the standard deviations of three datasets: (a–c) differences and RMSEs; (d–f) standard deviations.

    Table 2. RMSEs between observations and the BCC models in terms of climate-mean SM.

    Table 3. Mean standard deviations of the observed and BCCmodeled SM.

    The spatial features of surface SM anomalies over Eurasia from 1979 to 2012 are analyzed by using the EOF method. Figure 6 displays the first two principle components of the EOF analysis of GLDAS and the BCC model simulations. The EOF1 pattern in GLDAS, which explains 34.1% of the variance (Fig.6a), shows that there are significant negative signals in the western Siberia, central Asia and the Kamchatka region, while other areas in Eurasia are covered by weak positive signals. This pattern suggests that the surface soil tends to getting drier in western Siberia, central Asia and Kamchatka region, but slightly wetter in other areas. In the EOF2 pattern of GLDAS, with an explained variance of 8.9% (Fig.6b), there is a large-scale drying trend at high latitudes and mild wetting trend in other areas. The results are consistent with previous studies (Dong et al., 2007;Cheng et al., 2015; Gu et al., 2019). In BCC-CSM2-MR,the explained variance of the first two principle components reaches 20.0% in total (Figs. 6d and e). The EOF1 pattern of BCC-CSM2-MR that can capture negative values in western Siberia and central Asia has relatively high similarity with the EOF1 pattern of GLDAS. However, there are still some disagreements between BCC-CSM2-MR and GLDAS in the first principle component over the Arctic coastal region. This area is covered by positive values in BCC-CSM2-MR, which means the surface SM there tends to be wet. This trend is opposite to that of GLDAS. Therefore, BCC-CSM2-MR is able to capture the EOF1 spatial pattern of GLDAS, except at high latitudes in Arctic coastal regions. The EOF2 pattern of BCC-CSM2-MR, which shows that large-scale negative values cover the high latitudes, is similar to that of GLDAS, except that the gradient is not as obvious as shown in GLDAS. As for BCCCSM1.1m, EOF1 (Fig.6g) is unable to capture the spatial pattern of drying at high latitudes, as shown in the EOF1 of GLDAS and BCC-CSM2-MR. Neither western Siberia nor the Kamchatka region has the signs of a drying trend of surface SM. In the EOF2 pattern of BCC-CSM1.1m (Fig.6h),the high latitudes are covered by negative values, but their area and gradient are relatively smaller. In conclusion, compared to the previous-generation model, BCC-CSM2-MR is more skillful in describing the first two principle components of the EOF analysis of the surface SM anomaly over Eurasia from 1979 to 2012. The periodogram estimates of the spectra of the first principal component time series(PC1) of GLDAS and the BCC model simulations are shown in the rightmost column of Fig.6. The spectra above the red line are approved by the Markov “red noise” test. As shown in the panel, the corresponding period of GLDAS PC1 is 16.6 years. For the PC1s of BCC-CSM2-MR and BCC-CSM1.1m, the periods are 13.3 years and 8.3 years,respectively. Obviously, the period of PC1 in BCC-CSM2-MR is much closer to that in GLDAS. Therefore, BCCCSM2-MR is more capable of capturing the periodicity characteristics of the SM variation.

    Fig.6. The first two EOFs of annual mean soil moisture and periodogram estimates of the spectra of the PC1s: (a, d, g)EOF1 for GLDAS, BCC-CSM2-MR and BCC-CSM1.1m, respectively; (b, e, h) EOF2 for GLDAS, BCC-CSM2-MR and BCC-CSM1.1m, respectively; (c, f, i) periodogram estimates of the spectra of the PC1s of GLDAS, BCC-CSM2-MR and BCC-CSM1.1m (red lines: Markov “red noise” spectrum; green lines: upper confidence bound for Markov; blue lines: lower confidence bound for Markov).

    3.2.Evaluation of BCC-CSM2-MR among CMIP6 models

    We also compare the performance of BCC-CSM2-MR over Eurasia and China on annual and seasonal time scales with the other models participating in CMIP6, by using centered pattern correlation (computing anomalies from a central mean) and “amplitude of variations” (relative standard deviation). To make the comparison more comprehensive, BCC-CSM1.1m, the previous-generation model of BCC-CSM2-MR that participated in CMIP5, is also taken into account. Taylor diagrams (Taylor, 2001) are used to provide a visual representation of the aforementioned metrics.

    Centered pattern correlations between the surface SM simulated by the models and GLDAS are indicated by the azimuthal position of each dot in the Taylor diagrams. For the simulated annual mean surface SM (Fig.7a), the dots are relatively scattered. Correlations generally fall between 0.6 and 0.9 and tend to be clustered around 0.7. Compared with other CMIP6 models, BCC-CSM2-MR (r = 0.77), second only to EC-Earth3-Veg (r = 0.83) and MRI-ESM2-0 (r =0.78), performs at a high level. BCC-CSM2-MR also shows marked progress in correlation compared with BCCCSM1.1m (r = 0.62), with a rate of increase of 24.2%. In different seasons (Figs. 7b–e), although the correlations between the CMIP6 models and GLDAS are more variable,BCC-CSM2-MR is still well ahead of most models. Especially in winter (Fig.7b), BCC-CSM2-MR has the highest correlation coefficient with GLDAS (r = 0.77). Also, the coefficients of BCC-CSM2-MR are obviously higher than those of BCC-CSM1.1m in all seasons. The radial distance from the origin represents the standard deviation of the model simulation relative to the standard deviation of GLDAS (σ sim /σ obs). The closer to 1 the ratio is, the fewer biases in simulating SM variations the models have. In terms of the annual mean (Fig.7a), most models overestimate the standard deviation of surface SM considerably, with ratios above 1.25.BCC-CSM2-MR (σ sim / σ obs = 0.83) is one of the models able to give a comparatively accurate representation of the standard deviation. Meanwhile, compared with BCCCSM1.1m (σ sim / σ obs = 1.24), BCC-CSM2-MR also does a better job. During different seasons (Figs. 7b–e), the ratios of BCC-CSM2-MR are between 0.77 and 0.88 (the average ratio is 0.81), which places BCC-CSM2-MR at a better level among the CMIP6 models (the average ratio is 1.48) in representing similar standard deviations to GLDAS.However, compared to the previous-generation model, with slight overestimation (average ratio of 1.15), the improvements of BCC-CSM2-MR are not that obvious.

    Fig.7. Taylor diagrams for model-simulated surface soil moisture based on GLDAS over Eurasia: (a) annual mean; (b)winter; (c) spring; (d) summer; (e) autumn.

    Fig.8. Taylor diagrams for model-simulated surface soil moisture based on GLDAS over China: (a) annual mean; (b)winter; (c) spring; (d) summer; (e) autumn.

    For the area of China, the dots in Fig.8 are more clustered than in Fig.7, which means that the models perform relatively consistently over China. For the annual mean surface SM (Fig.8a), the correlation coefficients in the model simulations fall between 0.4 and 0.8. BCCCSM2-MR (r = 0.61) is at the mid-upper level among the CMIP6 models. Compared with BCC-CSM1.1m (r = 0.41),BCC-CSM2-MR has made an evident improvement. In the four seasons (Figs. 8b–e), the correlation coefficients of BCC-CSM2-MR always fall within the range of 0.5–0.8.The average is 0.64, which is higher than that of the CMIP6 models (r = 0.61), showing BCC-CSM2-MR to be more skillful in accurately and consistently representing the distribution of surface SM compared to most of the CMIP6 models.Meanwhile, BCC-CSM2-MR takes the lead in comparison with BCC-CSM1.1m (average correlation coefficient is 0.40). As for the ratios between the standard deviations of the model simulations and those of the reference data, the advantages of BCC-CSM2-MR are highlighted. The ratios of BCC-CSM2-MR are generally around 0.7–0.8 which are the closest to 1 among all the CMIP6 models. This indicates that BCC-CSM2-MR can capture the variability of the surface SM properly and maintain a minimum bias relative to the reference data. Compared with BCC-CSM1.1m,BCC-CSM2-MR carries forward the advantage of a similar standard deviation to that of GLDAS, and a slight improvement to the overestimation of the previous-generation model.

    In general, compared with other CMIP6 models and BCC-CSM1.1m, BCC-CSM2-MR is more competent in describing the distributions and variations of annual and seasonal SM, as shown by the relatively higher centered pattern correlation and standard deviation that is closer to that of the reference data, either in Eurasia or China.

    3.3.Role of precipitation in soil moisture simulation

    The above analysis shows that BCC-CSM2-MR has made certain progress in the simulation of surface SM compared to BCC-CSM1.1m. This might benefit from improvements to parameterization schemes in the component models, especially the land model implemented in BCC-CSM2-MR (Wu et al., 2019), such as the inclusion of a variable temperature threshold to determine soil water freeze–thaw,rather than a fixed temperature of 0 °C, a better calculation of snow cover fraction, and so on (Li et al., 2019). In the actual process of land–atmosphere coupling, there are many elements that have effects on the variations of SM, such as precipitation, temperature, wind, etc. Among these elements, precipitation is well known as the most vital and has thus been widely studied. The wetting of soil by precipitation, identified as the first part of the land–atmosphere feedback of water, is straightforward and intuitive and indisputably occurs in nature (Koster et al., 2003, 2004; Tawfik and Steiner, 2011). Numerous studies have shown that precipitation is the most direct and vital among the factors affecting SM. Zhang et al. (2008) used analysis data of precipitation and SM from GLDAS and pointed out that the strong land–atmosphere coupling lies mainly in semi-humid forest to grassland transition zones or in arid to semi-arid transition zones, including central Eurasia, northern China, etc.Research has also shown that the interaction between precipitation and SM exists in atmospheric general circulation models (Oglesby and Erickson III, 1989; Dirmeyer, 2000).

    The SM procedure adopted in BCC-CSM is almost the same as that in the NCAR Community Land Model (Oleson et al., 2004). SM is governed by infiltration, surface and sub-surface runoff, gradient diffusion, gravity, and root extraction through canopy transpiration. For one-dimensional vertical water flow in soils, the conservation of mass is stated as

    where θ is the volumetric soil water content, t is time and z is height above some datum in the soil column, q is the soil water flux, and e is the evapotranspiration loss. In the coupling between the land model and atmospheric model, liquid and solid precipitation from the atmospheric model will have an important effect on q and subsequently influence the simulation of SM. Therefore, the improvements in SM simulation in BCC-CSM2-MR may be attributable to a better simulation of precipitation. Influences of precipitation on SM in the models are thus discussed as follows.

    Differences in the climate mean of the surface SM (c.f.Fig.2) have shown that there are three regions where the difference values are significantly improved: Siberia(55°–64°N, 60°–87°E), Northeast China (40°–50°N, 120°–135°E), and the Yangtze–Huaihe River basin (28°–33°N,110°–121°E). The difference values in these three regions decrease by 32.7%, 30.0% and 20.6%, respectively. Hence,we select these three subregions (see the black sectors in Figs. 2a–c) and the whole of Eurasia as target areas to investigate the possible reasons behind the improvement in SM simulation, which has been linked with precipitation.

    Fig.9. Correlation coefficients between the annual mean time series of soil moisture and precipitation of BCC-CSM2-MR,BCC-CSM1.1m and reference data (GLDAS/CMAP) over Eurasia (EA), Siberia (SIB), Northeast China (NEC) and the Yangtze–Huaihe River basin (YH).

    Figure 9 shows the correlation coefficients between the time series of the annual mean surface SM and precipitation in the reference data (GLDAS for SM, CMAP for precipitation) and the BCC model simulations. In GLDAS/CMAP,the correlation coefficients range between 0.35 and 0.60,which are statistically significant at the 95% confidence level according to the Student’s t-test. As for the BCC model simulations, the correlation coefficients are generally significant in most areas except the whole of Eurasia.For Eurasia, which of course covers a wide range of longitudes and latitudes, the factors that influence SM in are more complex and diverse. Therefore, this complicated land–atmosphere interaction makes the proportion of the influence of precipitation on SM smaller and the correlation coefficients between them to be reduced further. Coupled climate models may not properly capture the complicated relationships between SM and factors of influence, and perhaps may even weaken the influence of precipitation excessively. Hence,the correlation coefficients in the BCC model simulations are quite low and even non-significant. As for the three subregions, the correlation coefficients of BCC-CSM2-MR in Siberia, Northeast China and the Yangtze–Huaihe River basin are 0.49, 0.67 and 0.87, respectively, while those of BCC-CSM1.1m are 0.60, 0.78 and 0.89, respectively. The former is apparently lower than the latter and much closer to the coefficients in GLDASCMAP (r = 0.52, 0.59 and 0.49, respectively). That is to say, the proportion of influence of precipitation on surface SM in BCC-CSM2-MR is not as much as in BCC-CSM1.1m. As seen on the annual time scale, the coefficients of BCC-CSM2-MR are more consistent with those of the reference data on the seasonal scale(not shown here). In other words, compared with the previous-generation model, BCC-CSM2-MR is able to represent,relatively realistically, the relationship between precipitation and surface SM in these three subregions.

    Fig.10. Scatterplots showing the anomalies of soil moisture in the reference data and BCC models, in conjunction with corresponding changes in precipitation, over (a) Eurasia (EA), (b) Siberia (SIB), (c) Northeast China (NEC) and (e) the Yangtze–Huaihe River basin (YH).

    The anomalies of surface SM in GLDAS and the BCC models, in conjunction with corresponding anomalies of precipitation in CMAP and the BCC models, over Eurasia and the three subregions, are shown in Fig.10. In Eurasia, the SM responses in GLDAS are positively correlated with the precipitation variations in CMAP (Fig.10a). By analyzing the slope of the regression line, we can conclude that BCCCSM2-MR is able to represent the correlation with a positive linear regression coefficient, but BCC-CSM1.1m shows an opposite correlation with a negative linear regression coefficient. In the three subregions (Figs. 10b–d), the dots are not as spread as they are in Eurasia. This means that precipitation anomalies play a more obvious and important role in increasing SM. In general, the regression lines in BCCCSM2-MR are closer to those in the reference data, as indicated by the similar spacing and trends of the dots. Taking Northeast China as an example, the linear regression coefficient of GLDASCMAP is 0.58, while the coefficients of BCC-CSM2-MR and BCC-CSM1.1m are 0.67 and 0.77,respectively. Obviously, the former is closer than the latter to the reference data. Therefore, compared to the previous version, BCC-CSM2-MR is more skillful in describing the relationship between precipitation and surface SM variations.

    We also evaluate the RMSE values of the mean precipitation between CMAP and the BCC model simulations on the annual time scale (Table 4). The RMSEs of BCC-CSM2-MR are lower than those of BCC-CSM1.1m to some extent.For Eurasia, the RMSE of the annual mean decreases by 9.6%, from 1.14 (BCC-CSM1.1m) to 1.03 (BCC-CSM2-MR). In the three subregions, the RMSEs of BCC-CSM2-MR tend to be obviously lower. The Yangtze–Huaihe River basin is the region with the most obvious improvement,where the RMSE is reduced by 40.4%. The progress is also significant on the seasonal time scale (not shown here).Hence, BCC-CSM2-MR is generally more skillful than BCC-CSM1.1m in realistically representing the variation of precipitation, as shown by the smaller deviation between the simulation and CMAP.

    From the above analyses, BCC-CSM2-MR is better able to properly describe the correlation between the surface SM and precipitation, the response of the surface SM variation to precipitation anomalies, and the variation of precipitation. These qualities may contribute to the better simulation of surface SM in BCC-CSM2-MR.

    4.Summary and discussion

    This study has evaluated the surface SM simulations of BCC-CSM2-MR over Eurasia on annual and seasonal time scales and compared them with its previous-generation version (BCC-CSM1.1m) and other CMIP6 models, as well as with GLDAS and station-based observations as reference data, and further discussed the possible reasons behind the improvement in surface SM simulations in relation to precipitation. The conclusions can be summarized as follows:

    BCC-CSM2-MR, which is participating in CMIP6, has clearly improved SM simulations compared with BCCCSM1.1m, part of CMIP5, not only in terms of spatial distributions of the climate mean surface SM, with higher pattern correlations with GLDAS, but also in terms of obviously reduced RMSEs. The standard deviations of BCCCSM2-MR also have closer spatial distributions to those of the reference data, and the RMSEs are generally lower. The first two principle components of BCC-CSM2-MR are able to capture the spatial patterns shown as EOF1 and EOF2 of GLDAS, while BCC-CSM1.1m struggles to describe the characteristics properly. Compared with the 8.3-year period of PC1 shown in BCC-CSM1.1m, the period in BCC-CSM2-MR (13.3 years) is closer to that in GLDAS (16.6 years).

    Among the other models participating in CMIP6, BCCCSM2-MR performs well in statistical diagnostics. The average centered correlation coefficient of BCC-CSM2-MR (r =0.77) is the third highest among the CMIP6 models. The average relative standard deviation is 0.83, closer to that of the reference data, which surmounts the overestimations in most of the CMIP6 models.

    The possible reasons behind the improvement in BCCCSM2-MR, as related to precipitation, have been discussed.BCC-CSM2-MR simulates correlation coefficients and regression coefficients that are closer to those of GLDASCMAP.The simulation of precipitation by BCC-CSM2-MR is also improved in terms of the RMSE.

    In addition to the above conclusions, we found that the improvement of BCC-CSM2-MR differs from season to season. Specifically, the performances in autumn and winter are better than in spring and summer. Taking the simulations of the climate mean SM for example, the pattern correlation coefficient between model and reference data increases by 4.3%, from 0.93 to 0.97, in autumn and winter, which is higher than that the rate of 3.7% in spring and summer. The reduction in RMSE is 8.0% in autumn and winter, while it is 6.6% in spring and summer. The different performances between the cold and warm seasons are possibly related to the optimization of parameterization schemes. In BCC_AVIM1.0 (the land component in BCC-CSM1.1m),liquid water freezes when the soil temperature decreases to 0°C, and the soil temperature will remain at 0°C until all the liquid water has frozen. However, liquid water can coexist with ice in the real world when the soil temperature is below 0°C. The relationship between soil water content and soil temperature is determined by the inherent characteristics of soil hydraulics and the thermodynamic equilibrium between soil water potential and soil temperature. Therefore, the method to calculate the soil freeze–thaw critical temperature used by Li and Sun (2008) is adopted in BCC_AVIM2.0 (the land component in BCC-CSM2-MR)to replace the unreasonable assumption used in BCC_AVIM1.0. In addition, the parameterization of snow cover fraction in BCC_AVIM2.0 is also adjusted for theabsence of fluctuating topography by taking into account the subgrid variability of topography in a model grid cell(Li et al., 2019). Thus, the optimization for the parameterization schemes of soil freeze–thaw processes and the snow cover fraction may result in BCC-CSM2-MR performing better in its autumn and winter season simulations.

    Table 4. RMSEs between the BCC models and CMAP data of annual mean precipitation over Eurasia and four subregions.

    Also of note is that the performance of BCC-CSM2-MR improves differently in different regions. Possible reasons related to parameterization schemes are discussed here.In Fig.2, the areas with significant improvement are mostly located in Eurasian mid–high latitudes (Siberian region and Northeast China). For the relatively higher latitudes, processes related to snow cover and soil water freeze–thaw are more significant and important than in the low latitudes.Therefore, the upgrading of related parameterization schemes mentioned above will make more obvious improvements in the high latitudes. It is worth noting that the Yangtze–Huaihe River basin is also an area showing significant improvement. This area is the main rice production region in China. The interactions between such agricultural land and the overlying atmosphere play an important role,where surface latent heat flux values are relatively large.However, the plant functional type of “crop ” in BCC_AVIM1.0 to represent rice paddies will underestimate the amount of surface evaporation, which is the important link between SM and the atmosphere. A new scheme for rice paddy fields was developed in BCC_AVIM2.0 to incorporate the addition of surface water above soil (Li et al.,2019). The essential difference in the calculation of latent heat flux between BCC_AVIM2.0 and the original crop scheme in BCC_AVIM1.0 lies in that there is no limit to the evaporation from a rice paddy in BCC_AVIM2.0. Therefore, the optimization of the parameterization over rice paddies may result in the model improvement in the Yangtze–Huaihe River basin.

    However, there are still some limitations in this study.First, we applied bilinear interpolation to regrid both the simulation data and reference data into a uniform resolution of 1° × 1° to make the evaluation more convenient, which might have affected the conclusions to a certain degree because of the inhomogeneity of SM. Although bilinear interpolation is a simple method and commonly used (Crow et al., 2012; Hsu et al., 2013; Yuan and Quiring, 2017), applying a more advanced downscaling or interpolation technique may provide better estimates of model-simulated SM.Second, this study has only discussed the correlation between SM and precipitation. In fact, the improvements in related land surface process schemes are also important, as well as the land–atmosphere coupling in other component models, and so precipitation is merely one of the mostrelated factors. Other factors need more and further investigations in the future.

    Acknowledgements.

    This work was supported by the National Key Research and Development Program of China (Grant Nos. 2018YFC1506004 and 2016YFA0602602).

    黄色 视频免费看| av天堂久久9| 日本vs欧美在线观看视频| 中国美女看黄片| 天天添夜夜摸| 99re6热这里在线精品视频| 久久久精品免费免费高清| 欧美精品高潮呻吟av久久| 飞空精品影院首页| 一区二区三区激情视频| 久久久久国产精品人妻aⅴ院 | 欧美不卡视频在线免费观看 | 亚洲av成人av| 人妻 亚洲 视频| 一级a爱片免费观看的视频| 日韩免费av在线播放| 欧美日韩黄片免| 激情在线观看视频在线高清 | 老汉色av国产亚洲站长工具| 亚洲第一av免费看| av线在线观看网站| 国产高清视频在线播放一区| 亚洲欧美激情在线| 黄色成人免费大全| 757午夜福利合集在线观看| 91在线观看av| 在线观看免费高清a一片| 久久人人爽av亚洲精品天堂| 日韩欧美在线二视频 | 国产精品综合久久久久久久免费 | 一区福利在线观看| 啦啦啦免费观看视频1| 国产又爽黄色视频| 啦啦啦免费观看视频1| 日韩成人在线观看一区二区三区| av线在线观看网站| 亚洲专区中文字幕在线| 丝袜在线中文字幕| 亚洲精品国产区一区二| 久久国产精品影院| 一边摸一边抽搐一进一出视频| 少妇猛男粗大的猛烈进出视频| 高潮久久久久久久久久久不卡| 在线播放国产精品三级| 欧美另类亚洲清纯唯美| 欧美一级毛片孕妇| 大片电影免费在线观看免费| 老司机福利观看| 美女扒开内裤让男人捅视频| 大码成人一级视频| 最新在线观看一区二区三区| 精品久久久久久电影网| 三级毛片av免费| 欧美日韩视频精品一区| 91国产中文字幕| 久久久久国内视频| 亚洲av欧美aⅴ国产| 亚洲国产欧美一区二区综合| 啪啪无遮挡十八禁网站| 在线看a的网站| 女性生殖器流出的白浆| 视频区欧美日本亚洲| 成人免费观看视频高清| 国产男女超爽视频在线观看| 人人妻人人澡人人爽人人夜夜| videosex国产| 五月开心婷婷网| videosex国产| 久久中文字幕人妻熟女| 久久久国产精品麻豆| 99国产精品免费福利视频| 啦啦啦 在线观看视频| 美女高潮到喷水免费观看| 久久久水蜜桃国产精品网| 在线天堂中文资源库| 美女高潮到喷水免费观看| 国产精品一区二区在线观看99| 视频区欧美日本亚洲| 国产男女超爽视频在线观看| 黑丝袜美女国产一区| 国产精品av久久久久免费| 中文字幕人妻丝袜制服| 亚洲七黄色美女视频| 亚洲aⅴ乱码一区二区在线播放 | 天天躁狠狠躁夜夜躁狠狠躁| 自线自在国产av| 国产精品影院久久| 男女午夜视频在线观看| 国产亚洲精品第一综合不卡| 亚洲全国av大片| 好男人电影高清在线观看| 亚洲av成人不卡在线观看播放网| 男人舔女人的私密视频| 久久久久久久国产电影| 丝瓜视频免费看黄片| 日韩成人在线观看一区二区三区| 精品一区二区三卡| 黄色女人牲交| 十八禁网站免费在线| 久久国产精品大桥未久av| а√天堂www在线а√下载 | 丝袜美足系列| 天天躁日日躁夜夜躁夜夜| 无限看片的www在线观看| 高清欧美精品videossex| 757午夜福利合集在线观看| 亚洲精品成人av观看孕妇| 美女视频免费永久观看网站| 精品电影一区二区在线| 日韩制服丝袜自拍偷拍| 国产男靠女视频免费网站| 涩涩av久久男人的天堂| 午夜福利欧美成人| 另类亚洲欧美激情| 人人妻人人爽人人添夜夜欢视频| 成年人免费黄色播放视频| 黑丝袜美女国产一区| 在线看a的网站| 国产黄色免费在线视频| 欧美日韩视频精品一区| 亚洲人成伊人成综合网2020| 精品一区二区三区四区五区乱码| 国产成人欧美在线观看 | av线在线观看网站| 在线观看一区二区三区激情| 精品视频人人做人人爽| 国产精品国产av在线观看| 亚洲久久久国产精品| 国产午夜精品久久久久久| 黑丝袜美女国产一区| 男女免费视频国产| 又大又爽又粗| 99热国产这里只有精品6| 国产成人免费观看mmmm| 欧美成人免费av一区二区三区 | 亚洲一区中文字幕在线| 女同久久另类99精品国产91| 午夜视频精品福利| 在线视频色国产色| 日韩欧美一区视频在线观看| 人妻一区二区av| 黄色视频,在线免费观看| 亚洲av熟女| 午夜精品在线福利| 国产区一区二久久| 国产精品欧美亚洲77777| 午夜免费鲁丝| 欧美乱码精品一区二区三区| 九色亚洲精品在线播放| 一进一出抽搐gif免费好疼 | 日韩免费av在线播放| 免费少妇av软件| 一本一本久久a久久精品综合妖精| 国产一区二区三区在线臀色熟女 | 中文字幕人妻丝袜制服| 男人的好看免费观看在线视频 | ponron亚洲| 亚洲av成人不卡在线观看播放网| 啦啦啦免费观看视频1| 亚洲av美国av| 久久精品人人爽人人爽视色| 免费人成视频x8x8入口观看| 久9热在线精品视频| 91麻豆av在线| 国产激情欧美一区二区| 亚洲自偷自拍图片 自拍| 精品国产超薄肉色丝袜足j| 亚洲美女黄片视频| 亚洲成人国产一区在线观看| 黑丝袜美女国产一区| 丝袜人妻中文字幕| 老司机午夜十八禁免费视频| 在线播放国产精品三级| 精品卡一卡二卡四卡免费| 亚洲五月色婷婷综合| 日本五十路高清| 国产精品国产av在线观看| videosex国产| 99riav亚洲国产免费| 国产成人一区二区三区免费视频网站| 成人三级做爰电影| 亚洲色图综合在线观看| 国产精品久久久久成人av| 老司机午夜十八禁免费视频| 国产精品一区二区免费欧美| 免费在线观看视频国产中文字幕亚洲| 精品一区二区三区视频在线观看免费 | 国产又爽黄色视频| 日韩三级视频一区二区三区| 母亲3免费完整高清在线观看| 中国美女看黄片| 交换朋友夫妻互换小说| 国产成人影院久久av| 国产精华一区二区三区| 欧美成人免费av一区二区三区| 最近最新中文字幕大全电影3| 欧美乱色亚洲激情| 欧美三级亚洲精品| 久久久久亚洲av毛片大全| 国产高清视频在线观看网站| 一夜夜www| 在线观看日韩欧美| 国内久久婷婷六月综合欲色啪| 99久久99久久久精品蜜桃| 精品国产亚洲在线| 一a级毛片在线观看| 久久香蕉国产精品| 在线观看午夜福利视频| 夜夜夜夜夜久久久久| 高潮久久久久久久久久久不卡| 少妇熟女aⅴ在线视频| 伊人久久大香线蕉亚洲五| 久久精品国产亚洲av香蕉五月| 黄色视频,在线免费观看| 午夜精品在线福利| 一卡2卡三卡四卡精品乱码亚洲| 国内精品久久久久精免费| 99久久九九国产精品国产免费| 国产乱人视频| 久久精品国产清高在天天线| 日本a在线网址| 欧美日韩福利视频一区二区| 男女那种视频在线观看| 亚洲最大成人中文| 神马国产精品三级电影在线观看| 高清日韩中文字幕在线| 亚洲七黄色美女视频| 少妇人妻精品综合一区二区 | 精品一区二区三区人妻视频| 国产主播在线观看一区二区| av在线天堂中文字幕| 精品一区二区三区视频在线 | 欧美av亚洲av综合av国产av| 国产黄a三级三级三级人| 18美女黄网站色大片免费观看| 中文字幕精品亚洲无线码一区| 日韩欧美国产在线观看| 国产精品99久久久久久久久| 日韩人妻高清精品专区| 在线免费观看的www视频| 久久久久久久精品吃奶| 亚洲欧美激情综合另类| 淫妇啪啪啪对白视频| 一个人看的www免费观看视频| 国产三级在线视频| 女生性感内裤真人,穿戴方法视频| 亚洲aⅴ乱码一区二区在线播放| 老司机午夜十八禁免费视频| 脱女人内裤的视频| av福利片在线观看| 国产午夜精品久久久久久一区二区三区 | 国产精品永久免费网站| 日韩国内少妇激情av| 久久久久亚洲av毛片大全| 亚洲内射少妇av| 亚洲成av人片免费观看| 中文字幕人妻熟人妻熟丝袜美 | 国产黄色小视频在线观看| 亚洲成人久久爱视频| 欧美一级a爱片免费观看看| 精品福利观看| 91av网一区二区| 国产精品精品国产色婷婷| 热99re8久久精品国产| 啦啦啦韩国在线观看视频| 日韩中文字幕欧美一区二区| 国产伦精品一区二区三区视频9 | 在线看三级毛片| 久久久久久久久中文| 每晚都被弄得嗷嗷叫到高潮| 丰满的人妻完整版| 啪啪无遮挡十八禁网站| 午夜福利成人在线免费观看| 欧美日韩瑟瑟在线播放| 亚洲国产精品999在线| 叶爱在线成人免费视频播放| 97超级碰碰碰精品色视频在线观看| 日韩大尺度精品在线看网址| 欧美日韩黄片免| 丰满人妻一区二区三区视频av | 天天一区二区日本电影三级| 精品人妻1区二区| 亚洲熟妇中文字幕五十中出| 免费看十八禁软件| 亚洲国产精品久久男人天堂| 中文资源天堂在线| 精品人妻1区二区| 亚洲无线在线观看| 欧美色视频一区免费| 成人18禁在线播放| 一个人看的www免费观看视频| 久久亚洲精品不卡| 亚洲黑人精品在线| 人妻久久中文字幕网| 特级一级黄色大片| 亚洲欧美日韩无卡精品| 国产黄色小视频在线观看| 成人鲁丝片一二三区免费| 在线播放无遮挡| 国产精品综合久久久久久久免费| 亚洲内射少妇av| 亚洲国产色片| 国产精品亚洲美女久久久| 日韩人妻高清精品专区| 成年女人永久免费观看视频| 99在线视频只有这里精品首页| 一级毛片女人18水好多| 悠悠久久av| 久久香蕉国产精品| 国产精品美女特级片免费视频播放器| 99热精品在线国产| 在线观看免费视频日本深夜| 18禁国产床啪视频网站| 国产亚洲av嫩草精品影院| 精品一区二区三区av网在线观看| 在线天堂最新版资源| 丰满的人妻完整版| 俺也久久电影网| 久久久久久久亚洲中文字幕 | 亚洲专区国产一区二区| 国产欧美日韩一区二区三| 一本久久中文字幕| 亚洲人成网站高清观看| av在线蜜桃| 亚洲国产欧洲综合997久久,| 1024手机看黄色片| 久久久精品大字幕| 日韩免费av在线播放| 成人特级av手机在线观看| xxxwww97欧美| 成人永久免费在线观看视频| 亚洲精品日韩av片在线观看 | 精品无人区乱码1区二区| 欧美乱码精品一区二区三区| 无遮挡黄片免费观看| 日日干狠狠操夜夜爽| 成人18禁在线播放| 亚洲人成网站在线播放欧美日韩| 美女高潮喷水抽搐中文字幕| 午夜福利视频1000在线观看| 天堂网av新在线| 久久久久久人人人人人| 久久久精品欧美日韩精品| 老司机午夜十八禁免费视频| 国产亚洲精品一区二区www| 搡老岳熟女国产| av视频在线观看入口| 亚洲第一电影网av| 亚洲中文字幕日韩| 天堂√8在线中文| 国产美女午夜福利| 韩国av一区二区三区四区| 在线视频色国产色| 色播亚洲综合网| 国产高清三级在线| 日韩免费av在线播放| 日韩欧美免费精品| 蜜桃久久精品国产亚洲av| 国产精品亚洲一级av第二区| 中文字幕人妻丝袜一区二区| 亚洲国产色片| 韩国av一区二区三区四区| eeuss影院久久| 欧美xxxx黑人xx丫x性爽| 日本一本二区三区精品| 一进一出抽搐gif免费好疼| 嫩草影院精品99| 成熟少妇高潮喷水视频| 麻豆成人av在线观看| 99精品欧美一区二区三区四区| 国产淫片久久久久久久久 | 精品一区二区三区视频在线 | 国产精品一区二区三区四区久久| 在线观看日韩欧美| 亚洲无线观看免费| 久久久成人免费电影| 成年人黄色毛片网站| 一级毛片高清免费大全| 色在线成人网| 精品久久久久久久末码| 亚洲在线自拍视频| 成人国产一区最新在线观看| 国产黄色小视频在线观看| 精品久久久久久久久久久久久| 午夜福利在线观看免费完整高清在 | 欧美乱妇无乱码| 亚洲国产欧美网| 少妇丰满av| 久久九九热精品免费| 中文字幕高清在线视频| 丁香欧美五月| 日本一二三区视频观看| 18+在线观看网站| 欧美在线黄色| 亚洲最大成人手机在线| 99久久99久久久精品蜜桃| 婷婷六月久久综合丁香| 最近在线观看免费完整版| 亚洲午夜理论影院| 男女下面进入的视频免费午夜| 久久久久国内视频| 99久久精品一区二区三区| 亚洲自拍偷在线| 十八禁人妻一区二区| a级一级毛片免费在线观看| 日本与韩国留学比较| 国产精品久久久久久人妻精品电影| 国产又黄又爽又无遮挡在线| av在线蜜桃| 性色av乱码一区二区三区2| 国产av一区在线观看免费| 国产精品亚洲一级av第二区| 女生性感内裤真人,穿戴方法视频| 99久久综合精品五月天人人| 国产视频一区二区在线看| 国产国拍精品亚洲av在线观看 | 国内精品久久久久精免费| 无遮挡黄片免费观看| 色精品久久人妻99蜜桃| 99久久精品一区二区三区| 18禁黄网站禁片免费观看直播| 久久性视频一级片| 免费av观看视频| 91麻豆精品激情在线观看国产| x7x7x7水蜜桃| 精品福利观看| 国产三级黄色录像| 欧洲精品卡2卡3卡4卡5卡区| 欧美日韩黄片免| 亚洲欧美日韩卡通动漫| 免费人成在线观看视频色| 亚洲精品成人久久久久久| 黄片大片在线免费观看| 内射极品少妇av片p| 欧美黄色片欧美黄色片| 亚洲第一欧美日韩一区二区三区| 中文亚洲av片在线观看爽| 天堂影院成人在线观看| 特级一级黄色大片| 搡老岳熟女国产| 久久久国产成人免费| 欧美日韩瑟瑟在线播放| 99在线人妻在线中文字幕| 老司机福利观看| 欧美一区二区精品小视频在线| 免费人成视频x8x8入口观看| 国产激情欧美一区二区| 久久这里只有精品中国| 久99久视频精品免费| 免费无遮挡裸体视频| 19禁男女啪啪无遮挡网站| 婷婷丁香在线五月| ponron亚洲| 亚洲黑人精品在线| 欧美高清成人免费视频www| 国产91精品成人一区二区三区| 又爽又黄无遮挡网站| 1000部很黄的大片| 成人午夜高清在线视频| 中文字幕久久专区| 少妇丰满av| x7x7x7水蜜桃| 亚洲成人免费电影在线观看| 国内毛片毛片毛片毛片毛片| 日韩大尺度精品在线看网址| 一a级毛片在线观看| 非洲黑人性xxxx精品又粗又长| 婷婷精品国产亚洲av| 叶爱在线成人免费视频播放| 国产在线精品亚洲第一网站| 男人舔女人下体高潮全视频| 在线观看免费午夜福利视频| 变态另类成人亚洲欧美熟女| 99国产综合亚洲精品| 亚洲av成人精品一区久久| 国产精品香港三级国产av潘金莲| 成人av在线播放网站| 亚洲最大成人中文| 成年女人永久免费观看视频| 级片在线观看| 亚洲国产精品999在线| 一进一出抽搐动态| 麻豆成人av在线观看| 18禁美女被吸乳视频| 国产精品电影一区二区三区| 国产伦精品一区二区三区视频9 | 怎么达到女性高潮| 国产不卡一卡二| 美女 人体艺术 gogo| 午夜两性在线视频| 亚洲av熟女| 欧美高清成人免费视频www| 欧美黑人欧美精品刺激| 狂野欧美激情性xxxx| 亚洲精品亚洲一区二区| 一本久久中文字幕| 国产精品一区二区三区四区免费观看 | 精品一区二区三区视频在线观看免费| 91在线精品国自产拍蜜月 | 国产精品免费一区二区三区在线| 成人亚洲精品av一区二区| 亚洲精品一区av在线观看| 久久久色成人| 熟女电影av网| 亚洲aⅴ乱码一区二区在线播放| 99久久99久久久精品蜜桃| 亚洲国产欧洲综合997久久,| 欧洲精品卡2卡3卡4卡5卡区| av国产免费在线观看| 99久久精品国产亚洲精品| 国产精品久久久久久久电影 | 波多野结衣巨乳人妻| 亚洲精品美女久久久久99蜜臀| 99热这里只有精品一区| 亚洲精品一区av在线观看| 亚洲人成电影免费在线| 国产精品亚洲一级av第二区| 男女视频在线观看网站免费| 韩国av一区二区三区四区| 国产三级黄色录像| 亚洲av成人精品一区久久| 男女做爰动态图高潮gif福利片| 欧美黄色淫秽网站| 国产精品影院久久| 熟女人妻精品中文字幕| 1000部很黄的大片| 免费看十八禁软件| 99久久九九国产精品国产免费| 免费观看人在逋| 午夜日韩欧美国产| 亚洲熟妇熟女久久| 精品久久久久久久久久免费视频| 久久99热这里只有精品18| 色av中文字幕| 村上凉子中文字幕在线| 51国产日韩欧美| 国内精品美女久久久久久| 99久久无色码亚洲精品果冻| 又黄又爽又免费观看的视频| 欧美一区二区亚洲| 麻豆一二三区av精品| 国产精品三级大全| 高潮久久久久久久久久久不卡| 色综合亚洲欧美另类图片| 精品99又大又爽又粗少妇毛片 | 成年免费大片在线观看| 1000部很黄的大片| 美女cb高潮喷水在线观看| 夜夜看夜夜爽夜夜摸| 在线观看午夜福利视频| 久久国产精品影院| 欧美黑人巨大hd| 国内精品一区二区在线观看| 日本a在线网址| 日本熟妇午夜| 亚洲av第一区精品v没综合| 成人国产一区最新在线观看| 亚洲最大成人手机在线| 国产精品1区2区在线观看.| 99久久无色码亚洲精品果冻| 欧美不卡视频在线免费观看| 搡老妇女老女人老熟妇| 国产伦在线观看视频一区| 亚洲狠狠婷婷综合久久图片| av专区在线播放| 波多野结衣高清作品| 国产 一区 欧美 日韩| 国产精品,欧美在线| 高清在线国产一区| 日韩欧美一区二区三区在线观看| 亚洲男人的天堂狠狠| 热99在线观看视频| 亚洲人成网站在线播| 3wmmmm亚洲av在线观看| 熟女少妇亚洲综合色aaa.| 亚洲精品成人久久久久久| 婷婷丁香在线五月| 国产真人三级小视频在线观看| 日本黄色片子视频| 国产野战对白在线观看| 亚洲欧美激情综合另类| av中文乱码字幕在线| 国产97色在线日韩免费| 日韩欧美精品v在线| 蜜桃亚洲精品一区二区三区| 亚洲成人精品中文字幕电影| 中出人妻视频一区二区| 日韩 欧美 亚洲 中文字幕| 午夜福利成人在线免费观看| 少妇丰满av| 国产精品 欧美亚洲| 免费av毛片视频| 两个人看的免费小视频| 国产精品免费一区二区三区在线| 真实男女啪啪啪动态图| 免费av毛片视频| 国产在线精品亚洲第一网站| 一二三四社区在线视频社区8| svipshipincom国产片| 国产欧美日韩一区二区三| 国产在视频线在精品| 国产精品久久久久久精品电影| av黄色大香蕉| 亚洲av免费在线观看| 午夜精品一区二区三区免费看| 九色成人免费人妻av| 人人妻人人看人人澡| 超碰av人人做人人爽久久 | 日韩精品青青久久久久久| 久久久久性生活片| 亚洲成人精品中文字幕电影| 亚洲欧美精品综合久久99| 亚洲片人在线观看| xxx96com| 国产精品三级大全|