• <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不卡久久| 不卡一级毛片| 午夜精品一区二区三区免费看| 国产亚洲av嫩草精品影院| 亚洲国产欧洲综合997久久,| 少妇人妻一区二区三区视频| 一个人看视频在线观看www免费| 成熟少妇高潮喷水视频| 久久天躁狠狠躁夜夜2o2o| 精品一区二区免费观看| av.在线天堂| 午夜免费男女啪啪视频观看 | 97超碰精品成人国产| 久久精品国产99精品国产亚洲性色| 国产精品女同一区二区软件| 亚洲国产精品sss在线观看| 国产午夜精品论理片| 成年女人毛片免费观看观看9| 在线观看美女被高潮喷水网站| 色噜噜av男人的天堂激情| 欧美区成人在线视频| 亚洲三级黄色毛片| 精品人妻熟女av久视频| 久久综合国产亚洲精品| 国产一区二区三区av在线 | 老女人水多毛片| 国产欧美日韩一区二区精品| 九九热线精品视视频播放| av天堂在线播放| 伦理电影大哥的女人| 亚洲在线自拍视频| 亚洲国产精品成人久久小说 | 久久精品影院6| 久久人人精品亚洲av| 搡老熟女国产l中国老女人| 亚洲av免费高清在线观看| 亚洲熟妇熟女久久| 91av网一区二区| 欧美日韩综合久久久久久| 精品久久久噜噜| 夜夜夜夜夜久久久久| 亚洲成人久久爱视频| eeuss影院久久| 大又大粗又爽又黄少妇毛片口| 热99re8久久精品国产| 日韩欧美免费精品| 我的女老师完整版在线观看| 一区二区三区高清视频在线| 国产成人一区二区在线| 99久久成人亚洲精品观看| 国产精品一区www在线观看| 欧美日韩乱码在线| 国产成人91sexporn| 又粗又爽又猛毛片免费看| 99久久九九国产精品国产免费| 日韩大尺度精品在线看网址| 久久草成人影院| 日本欧美国产在线视频| 一夜夜www| 哪里可以看免费的av片| 亚洲精品日韩av片在线观看| 寂寞人妻少妇视频99o| 国产精品亚洲美女久久久| 日本一本二区三区精品| 欧美3d第一页| 日本一二三区视频观看| 日韩欧美免费精品| 国产精品人妻久久久影院| 亚洲国产精品sss在线观看| 色哟哟哟哟哟哟| 成人高潮视频无遮挡免费网站| 亚洲aⅴ乱码一区二区在线播放| 18禁黄网站禁片免费观看直播| 免费观看在线日韩| 欧美激情国产日韩精品一区| 又黄又爽又刺激的免费视频.| 国产精品一二三区在线看| 午夜福利视频1000在线观看| 在线观看美女被高潮喷水网站| 国产高清视频在线观看网站| 一级a爱片免费观看的视频| 成人毛片a级毛片在线播放| 国产伦一二天堂av在线观看| 国产色爽女视频免费观看| 成年免费大片在线观看| 美女cb高潮喷水在线观看| 亚洲aⅴ乱码一区二区在线播放| 亚洲av成人av| 日本 av在线| av天堂在线播放| 亚洲国产高清在线一区二区三| 国产精品99久久久久久久久| 中文在线观看免费www的网站| a级毛色黄片| 亚洲最大成人手机在线| 观看美女的网站| 长腿黑丝高跟| 国产精品av视频在线免费观看| 欧美区成人在线视频| 亚洲欧美日韩东京热| 成人一区二区视频在线观看| 晚上一个人看的免费电影| 最近2019中文字幕mv第一页| 99riav亚洲国产免费| 嫩草影院新地址| 久久精品综合一区二区三区| 香蕉av资源在线| 精华霜和精华液先用哪个| 亚洲精品日韩在线中文字幕 | 亚洲成人av在线免费| 一级毛片aaaaaa免费看小| 女的被弄到高潮叫床怎么办| 国产黄色视频一区二区在线观看 | 婷婷精品国产亚洲av| 日韩成人伦理影院| 99热精品在线国产| 亚洲精品乱码久久久v下载方式| 少妇熟女欧美另类| 简卡轻食公司| 禁无遮挡网站| 男女视频在线观看网站免费| 免费看美女性在线毛片视频| 国产av不卡久久| 国产精品一及| 少妇被粗大猛烈的视频| 2021天堂中文幕一二区在线观| 12—13女人毛片做爰片一| .国产精品久久| 国产免费男女视频| 欧美高清成人免费视频www| 国内精品久久久久精免费| 精品久久久久久久久久久久久| 亚洲欧美精品综合久久99| 久久精品夜色国产| 久久久国产成人免费| 男插女下体视频免费在线播放| 免费搜索国产男女视频| 国产黄色小视频在线观看| av在线天堂中文字幕| 91久久精品国产一区二区三区| 亚洲精品粉嫩美女一区| 亚洲精品色激情综合| 我的女老师完整版在线观看| 久久亚洲精品不卡| 午夜a级毛片| 久久精品国产自在天天线| 国产真实伦视频高清在线观看| 中文字幕av成人在线电影| 国产成人精品久久久久久| 国产亚洲精品久久久久久毛片| 成年版毛片免费区| 伊人久久精品亚洲午夜| 亚洲中文日韩欧美视频| 老女人水多毛片| 中文在线观看免费www的网站| 国产av一区在线观看免费| 亚洲国产精品久久男人天堂| 精品久久久久久成人av| 如何舔出高潮| 女的被弄到高潮叫床怎么办| 97超级碰碰碰精品色视频在线观看| 午夜激情福利司机影院| 亚洲国产欧美人成| 国产成人福利小说| 国产国拍精品亚洲av在线观看| 欧美精品国产亚洲| 国产精品一区二区三区四区久久| 亚洲av美国av| 亚洲av美国av| 日韩欧美三级三区| 国产精品国产三级国产av玫瑰| 搡女人真爽免费视频火全软件 | 国产av不卡久久| 免费高清视频大片| АⅤ资源中文在线天堂| 五月玫瑰六月丁香| 国产男人的电影天堂91| 亚州av有码| 日韩三级伦理在线观看| 国产男人的电影天堂91| 少妇的逼好多水| 国产精品久久久久久av不卡| 一边摸一边抽搐一进一小说| 亚洲精品影视一区二区三区av| av卡一久久| 亚洲婷婷狠狠爱综合网| 成人无遮挡网站| 一级毛片我不卡| 好男人在线观看高清免费视频| 日韩高清综合在线| 一级毛片我不卡| 午夜免费激情av| 99在线视频只有这里精品首页| av女优亚洲男人天堂| 欧美xxxx性猛交bbbb| 最近的中文字幕免费完整| 精品无人区乱码1区二区| 18禁裸乳无遮挡免费网站照片| 天天躁日日操中文字幕| 高清毛片免费观看视频网站| 美女cb高潮喷水在线观看| 一本一本综合久久| 亚洲欧美中文字幕日韩二区| 亚洲国产精品成人久久小说 | 美女内射精品一级片tv| 亚洲欧美清纯卡通| 亚洲国产色片| 少妇人妻精品综合一区二区 | 麻豆国产97在线/欧美| 噜噜噜噜噜久久久久久91| 一级毛片久久久久久久久女| av在线播放精品| 久久午夜福利片| 亚洲va在线va天堂va国产| 免费观看的影片在线观看| 久久九九热精品免费| 午夜爱爱视频在线播放| 久久久久性生活片| 久久久国产成人精品二区| 特大巨黑吊av在线直播| 国产久久久一区二区三区| 亚洲av成人精品一区久久| 日日摸夜夜添夜夜添小说| 中国美白少妇内射xxxbb| 国产精品国产高清国产av| 你懂的网址亚洲精品在线观看 | 18禁裸乳无遮挡免费网站照片| 男人和女人高潮做爰伦理| 热99re8久久精品国产| 亚洲欧美精品综合久久99| 老师上课跳d突然被开到最大视频| 国产成人a区在线观看| 18+在线观看网站| 亚洲av五月六月丁香网| 国产成人精品久久久久久| 99国产极品粉嫩在线观看| 国产不卡一卡二| 国产精品免费一区二区三区在线| 国产一区二区激情短视频| 男女之事视频高清在线观看| 在线观看午夜福利视频| 熟女人妻精品中文字幕| 成年免费大片在线观看| 国产精品永久免费网站| 精品国内亚洲2022精品成人| 国产精品亚洲美女久久久| 在线观看66精品国产| 日本免费一区二区三区高清不卡| 亚洲美女黄片视频| 日韩制服骚丝袜av| 有码 亚洲区| 国产高清有码在线观看视频| 亚洲成a人片在线一区二区| 国产在视频线在精品| 国内久久婷婷六月综合欲色啪| 人妻丰满熟妇av一区二区三区| 日韩欧美在线乱码| 天天一区二区日本电影三级| 你懂的网址亚洲精品在线观看 | av女优亚洲男人天堂| 免费看av在线观看网站| 又爽又黄a免费视频| 在线国产一区二区在线| 免费无遮挡裸体视频| 成人毛片a级毛片在线播放| 婷婷色综合大香蕉| 亚洲自偷自拍三级| 日本免费a在线| 国内精品宾馆在线| 又爽又黄无遮挡网站| 观看免费一级毛片| 久久亚洲国产成人精品v| 在线观看免费视频日本深夜| 午夜福利在线在线| 国产私拍福利视频在线观看| 淫秽高清视频在线观看| 欧美成人一区二区免费高清观看| 在线免费观看不下载黄p国产| 亚洲久久久久久中文字幕| 成人欧美大片| 亚洲av二区三区四区| 91久久精品电影网| 国产不卡一卡二| 亚洲内射少妇av| 91在线精品国自产拍蜜月| 日韩亚洲欧美综合| 晚上一个人看的免费电影| 大又大粗又爽又黄少妇毛片口| 亚洲人成网站在线播放欧美日韩| 久久久午夜欧美精品| 一区福利在线观看| 欧美一区二区国产精品久久精品| 亚洲aⅴ乱码一区二区在线播放| 乱系列少妇在线播放| 精品99又大又爽又粗少妇毛片| 美女被艹到高潮喷水动态| 亚洲av五月六月丁香网| 亚洲一区二区三区色噜噜| 看非洲黑人一级黄片| 麻豆国产97在线/欧美| 69av精品久久久久久| 一进一出抽搐动态| 亚洲国产精品国产精品| 欧美激情国产日韩精品一区| 1000部很黄的大片| 国产精品av视频在线免费观看| 嫩草影院精品99| 在现免费观看毛片| 亚洲不卡免费看| 国产伦在线观看视频一区| 国产高清视频在线播放一区| 久久精品夜夜夜夜夜久久蜜豆| 嫩草影院精品99| 日韩亚洲欧美综合| 在线播放国产精品三级| 亚洲国产日韩欧美精品在线观看| 日本黄大片高清| 九九久久精品国产亚洲av麻豆| 日韩欧美免费精品| 日韩 亚洲 欧美在线| 男女那种视频在线观看| 老熟妇乱子伦视频在线观看| 人人妻,人人澡人人爽秒播| 高清日韩中文字幕在线| 日日干狠狠操夜夜爽| 少妇人妻精品综合一区二区 | 国产成人a∨麻豆精品| 亚洲成a人片在线一区二区| 免费人成在线观看视频色| 亚洲中文字幕日韩| 日本 av在线| 日本爱情动作片www.在线观看 | av在线蜜桃| 亚洲综合色惰| 女人十人毛片免费观看3o分钟| 韩国av在线不卡| 亚洲丝袜综合中文字幕| 国产免费一级a男人的天堂| 亚洲精品色激情综合| 国产精品久久视频播放| 欧美成人免费av一区二区三区| 亚洲国产日韩欧美精品在线观看| 成人毛片a级毛片在线播放| 国产成年人精品一区二区| 国模一区二区三区四区视频| 欧美另类亚洲清纯唯美| 伦精品一区二区三区| 波多野结衣高清无吗| 精品久久久久久久久av| 亚洲av免费高清在线观看| 老师上课跳d突然被开到最大视频| 婷婷亚洲欧美| www日本黄色视频网| 观看美女的网站| 亚洲一区高清亚洲精品| 在线播放无遮挡| 波多野结衣高清无吗| 国产一区二区激情短视频| 最近在线观看免费完整版| 精品久久久久久久久久免费视频| 女同久久另类99精品国产91| 精品国产三级普通话版| 久久精品夜色国产| 成人一区二区视频在线观看| 色哟哟哟哟哟哟| 亚洲自偷自拍三级| 给我免费播放毛片高清在线观看| 级片在线观看| 插逼视频在线观看| 91麻豆精品激情在线观看国产| 国产精品久久视频播放| 亚洲不卡免费看| 亚洲经典国产精华液单| 久久精品久久久久久噜噜老黄 | 欧美+日韩+精品| 国产精品福利在线免费观看| 综合色av麻豆| aaaaa片日本免费| 欧美日本视频| 免费看a级黄色片| 深爱激情五月婷婷| 欧美色欧美亚洲另类二区| 人人妻人人看人人澡| 波多野结衣高清无吗| 中文字幕熟女人妻在线| 免费在线观看成人毛片| 尾随美女入室| 2021天堂中文幕一二区在线观| 深爱激情五月婷婷| 亚洲精品456在线播放app| 亚洲国产精品合色在线| 日本 av在线| 欧美成人一区二区免费高清观看| 尤物成人国产欧美一区二区三区| 在线a可以看的网站| 97超视频在线观看视频| 男人和女人高潮做爰伦理| 国产精品电影一区二区三区| 午夜福利高清视频| 狂野欧美激情性xxxx在线观看| 麻豆精品久久久久久蜜桃| 亚洲,欧美,日韩| 国产不卡一卡二| 欧美xxxx黑人xx丫x性爽| 午夜久久久久精精品| 日本黄色视频三级网站网址| 长腿黑丝高跟| 秋霞在线观看毛片| 国产精品亚洲美女久久久| 欧美日韩精品成人综合77777| 91久久精品国产一区二区成人| 亚洲人成网站在线播| 成人亚洲欧美一区二区av| 国产精品人妻久久久影院| 久99久视频精品免费| 尾随美女入室| 在线观看一区二区三区| 国产精品久久视频播放| 国产精品一区二区免费欧美| 国产高潮美女av| 91久久精品国产一区二区三区| 欧美丝袜亚洲另类| 亚洲,欧美,日韩| 国产精品av视频在线免费观看| 日本一二三区视频观看| 国产欧美日韩精品一区二区| 国产男靠女视频免费网站| а√天堂www在线а√下载| 欧美高清性xxxxhd video| 国产精品亚洲一级av第二区| 免费高清视频大片| 免费在线观看成人毛片| 久久久久国产精品人妻aⅴ院| 一进一出好大好爽视频| 亚洲国产色片| 欧美人与善性xxx| 日本 av在线| 久久久色成人| 亚洲国产日韩欧美精品在线观看| 日韩中字成人| 美女cb高潮喷水在线观看| 国产毛片a区久久久久| 久久久久久久久久黄片| 麻豆成人午夜福利视频| 亚洲欧美成人精品一区二区| 精品久久久久久久久亚洲| 亚洲精品国产av成人精品 | 国产精品,欧美在线| 最近2019中文字幕mv第一页| 国产精品乱码一区二三区的特点| 日韩精品有码人妻一区| 在线天堂最新版资源| 国产一区二区在线av高清观看| 久久精品国产99精品国产亚洲性色| 在线观看66精品国产| 精品久久国产蜜桃| 极品教师在线视频| 亚洲国产高清在线一区二区三| 九九热线精品视视频播放| 成人亚洲精品av一区二区| 国产极品精品免费视频能看的| 中文字幕熟女人妻在线| 高清毛片免费看| ponron亚洲| 欧美zozozo另类| 99热精品在线国产| 日本色播在线视频| videossex国产| 少妇的逼水好多| 搡女人真爽免费视频火全软件 | 国产 一区 欧美 日韩| 两个人的视频大全免费| 69人妻影院| 中文字幕久久专区| 精华霜和精华液先用哪个| 亚洲av免费高清在线观看| 国产精品99久久久久久久久| 一进一出好大好爽视频| 少妇人妻精品综合一区二区 | 国产精品一及| 国产精品一二三区在线看| 国产高清不卡午夜福利| 亚洲精品国产成人久久av| 直男gayav资源| 国产精品久久久久久亚洲av鲁大| 色播亚洲综合网| 精品国内亚洲2022精品成人| av国产免费在线观看| 1024手机看黄色片| 亚州av有码| 久久草成人影院| 两个人的视频大全免费| av视频在线观看入口| 国产av麻豆久久久久久久| 午夜a级毛片| 午夜激情福利司机影院| 亚洲专区国产一区二区| 亚洲最大成人av| 国产又黄又爽又无遮挡在线| 欧美日本亚洲视频在线播放| 日本一二三区视频观看| 欧美绝顶高潮抽搐喷水| 成人三级黄色视频| 午夜老司机福利剧场| 精品人妻一区二区三区麻豆 | 精品一区二区免费观看| 欧美色欧美亚洲另类二区| 你懂的网址亚洲精品在线观看 | 国产色爽女视频免费观看| 波多野结衣高清无吗| 日韩高清综合在线| 中文资源天堂在线| 亚洲精品色激情综合| 人妻久久中文字幕网| 欧美+亚洲+日韩+国产| 高清毛片免费观看视频网站| 熟妇人妻久久中文字幕3abv| 亚洲熟妇熟女久久| 欧洲精品卡2卡3卡4卡5卡区| 99国产极品粉嫩在线观看| 毛片一级片免费看久久久久| 麻豆一二三区av精品| 国产精品,欧美在线| 如何舔出高潮| 天天躁夜夜躁狠狠久久av| 99国产精品一区二区蜜桃av| 真人做人爱边吃奶动态| 国产精品国产三级国产av玫瑰| 特大巨黑吊av在线直播| 日韩中字成人| a级毛色黄片| 老熟妇仑乱视频hdxx| 久久九九热精品免费| 老熟妇乱子伦视频在线观看| 噜噜噜噜噜久久久久久91| 国产精品久久电影中文字幕| 日韩人妻高清精品专区| 午夜福利视频1000在线观看| 国产精品亚洲美女久久久| 久久久久久国产a免费观看| av黄色大香蕉| 最新中文字幕久久久久| 日本黄色片子视频| 大型黄色视频在线免费观看| 国产亚洲精品久久久com| 人人妻人人澡人人爽人人夜夜 | 亚洲美女黄片视频| 最近2019中文字幕mv第一页| 日本-黄色视频高清免费观看| 天天躁日日操中文字幕| 日韩国内少妇激情av| 久久久精品欧美日韩精品| 亚洲无线在线观看| 老女人水多毛片| 日韩制服骚丝袜av| 如何舔出高潮| 国产精品无大码| 亚洲成a人片在线一区二区| 日韩av不卡免费在线播放| 国产欧美日韩一区二区精品| 国产精品久久久久久久电影| 18禁在线无遮挡免费观看视频 | 最近视频中文字幕2019在线8| 一级毛片aaaaaa免费看小| 亚洲18禁久久av| 十八禁国产超污无遮挡网站| 国产精品av视频在线免费观看| 晚上一个人看的免费电影| 婷婷亚洲欧美| 69人妻影院| 熟妇人妻久久中文字幕3abv| 淫秽高清视频在线观看| 国产激情偷乱视频一区二区| 久久人人爽人人爽人人片va| 精品一区二区三区av网在线观看| 国产精品一及| 国产又黄又爽又无遮挡在线| 国产精品久久久久久精品电影| 午夜福利成人在线免费观看| 亚洲成人中文字幕在线播放| 国产单亲对白刺激| 亚洲欧美日韩无卡精品| 国产精品久久久久久精品电影| 夜夜夜夜夜久久久久| 欧美又色又爽又黄视频| 少妇人妻精品综合一区二区 | 国产成人a区在线观看| 黑人高潮一二区| 天堂网av新在线| 亚洲熟妇中文字幕五十中出| 淫妇啪啪啪对白视频| 欧美精品国产亚洲| 天堂网av新在线| 国产国拍精品亚洲av在线观看| 精品久久久噜噜| 久久人人爽人人爽人人片va| 十八禁国产超污无遮挡网站| 美女黄网站色视频| 亚洲人成网站高清观看| 国产精品人妻久久久影院| 久久久久精品国产欧美久久久| 免费av毛片视频| 老女人水多毛片| 国产成年人精品一区二区| 亚洲欧美日韩高清在线视频| 久久精品影院6| av专区在线播放| 变态另类丝袜制服| 国产激情偷乱视频一区二区| 国产一区二区亚洲精品在线观看| 亚洲人成网站在线观看播放| 综合色av麻豆| 国产亚洲精品综合一区在线观看| 变态另类丝袜制服|