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

    Skill Assessment of Copernicus Climate Change Service Seasonal Ensemble Precipitation Forecasts over Iran

    2021-04-20 04:02:18MasoudNOBAKHTBahramSAGHAFIANandSalehAMINYAVARI
    Advances in Atmospheric Sciences 2021年3期

    Masoud NOBAKHT, Bahram SAGHAFIAN*, and Saleh AMINYAVARI

    1Department of Civil Engineering, Science and Research Branch, Islamic Azad University, Tehran 1477893855, Iran

    2Department of Civil Engineering, Chalous Branch, Islamic Azad University, Chalous 46615/397, Iran

    (Received 2 February 2020; revised 15 October 2020; accepted 19 October 2020)

    ABSTRACT Medium to long-term precipitation forecasting plays a pivotal role in water resource management and development of warning systems. Recently, the Copernicus Climate Change Service (C3S) database has been releasing monthly forecasts for lead times of up to three months for public use. This study evaluated the ensemble forecasts of three C3S models over the period 1993-2017 in Iran’s eight classified precipitation clusters for one- to three-month lead times. Probabilistic and non-probabilistic criteria were used for evaluation. Furthermore, the skill of selected models was analyzed in dry and wet periods in different precipitation clusters. The results indicated that the models performed best in western precipitation clusters, while in the northern humid cluster the models had negative skill scores. All models were better at forecasting upper-tercile events in dry seasons and lower-tercile events in wet seasons. Moreover, with increasing lead time, the forecast skill of the models worsened. In terms of forecasting in dry and wet years, the forecasts of the models were generally close to observations, albeit they underestimated several severe dry periods and overestimated a few wet periods.Moreover, the multi-model forecasts generated via multivariate regression of the forecasts of the three models yielded better results compared with those of individual models. In general, the ECMWF and UKMO models were found to be appropriate for one-month-ahead precipitation forecasting in most clusters of Iran. For the clusters considered in Iran and for the long-range system versions considered, the Météo France model had lower skill than the other models.

    Key words: ensemble forecasts, Copernicus Climate Change Service, long-term forecasting, model evaluation, Iran

    1. Introduction

    Accurate precipitation forecasting is a key component in water resources decision making. In particular, accurate and timely forecasting of monthly and seasonal precipitation and streamflow improves the performance of meteorological and hydrological drought early warning systems.Moreover, meteorological forecasts increasingly rely on precipitation numerical models.

    Seasonal forecasts were initially based on estimates of signal-to-noise ratios that assumed full knowledge of ocean and/or land conditions, such that the variance of climate variables related to lower boundary forcing represented the signal. However, it has been shown that the interaction between the atmosphere, sea ice, land and ocean are also important (Doblas-Reyes et al., 2013). Seasonal climate forecasting centers around the world now routinely run coupled ocean-atmosphere general circulation models (GCMs).GCMs parameterize physics on coarse grids but involve coupling of ocean and atmosphere modules. The main aim of GCMs is to produce intra- to interseasonal forecasts driven by the slowly evolving boundary conditions, such as sea surface temperature (SST) (Duan et al., 2019).

    Ensemble forecasting is a technique used in numerical forecasts such that, instead of a single forecast, a series of forecasts depicting a wide range of future possible conditions in the atmosphere are produced. Nowadays, various centers produce ensemble forecasts of meteorological variables,such as precipitation and temperature, at different lead times from hourly to yearly, on a global scale, using numerical solutions of atmospheric hydrodynamic equations based on different initial conditions. Such forecasts have offered new opportunities in the water management sector. The Copernicus Climate Change Service (C3S) database has recently released the products of several significant European centers.

    The C3S database combines weather system observations and provides comprehensive information on the past,present and future weather conditions. The service, run by ECMWF on behalf of the European Union, provides seasonal forecasting protocols on its website. The C3S Seasonal Service is publicly available through the Climate Data Store. Since the C3S database has been only recently launched, not much research on its evaluation has been reported. Forecast skill assessment can provide valuable information for forecasters and developers such that decision makers may adopt appropriate strategies to mitigate climate-related risks.

    Manzanas et al. (2019) bias-corrected and calibrated C3S ensemble seasonal forecasts for Europe and Southeast Asia. The study adopted the ECMWF-SEAS5, UKMOGloSea5 and Meteo France System5 models with a onemonth lead time. For post-processing, simple bias adjustment (BA) and more sophisticated ensemble recalibration(RC) methods were employed, with the RC methods improving the reliability and outperforming the BA methods. The seasonal precipitation forecasts at the global scale were evaluated by Manzanas et al. (2014), who concluded that the best predictive skills were obtained in September-October and the poorest in March-May. MacLachlan et al. (2015) found that UKMO-GloSea5 had great forecast skill and reliability in predicting the North Atlantic Oscillation and Arctic Oscillation. Li and Robertson (2015) evaluated the ensemble forecast skill of the ECMWF, Japan Meteorological Agency(JMA) and CFSv2 models (seasonal center model of the NCEP) for summer and up to a four-week lead time. All three models provided good results for the first-week lead time. The forecast skill significantly declined in most clusters (except in some tropical clusters) during the second to fourth weeks. The precipitation forecast skill of the ECMWF model was significantly better than those of the other two models, especially for the third and fourth weeks.Shirvani and Landman (2016) studied the seasonal precipitation forecast skills of the North American Multi-Model Ensemble (NMME), two coupled ocean-atmosphere models and one two-tiered model, over Iran. They found low forecast skill for most models at all lead times, except in the October-December period at lead times of up to three months.

    Crochemore et al. (2017) evaluated the seasonal forecast skill of precipitation and river flow in 16 watersheds in France. The hindcasts of the ECMWF seasonal precipitation model (System 4) with a 90-day lead time were evaluated and post-processed using SAFRAN data (Météo France reanalysis product). They employed linear scaling (LS) and monthly and annual distribution mapping to post-process the raw precipitation data. The results indicated an increase in precipitation forecast skill after applying post-processing methods. Bett et al. (2017) assessed the skill and reliability of wind speed in GloSea5 seasonal forecasts corresponding to winter and summer seasons over China. The results showed that the winter mean wind speed was skillfully forecasted along the coast of the South China Sea. Lucatero et al. (2018) examined the skill of raw and post-processed ensemble seasonal meteorological forecasts in Denmark.They took advantage of the LS and quantile mapping (QM)techniques to bias-correct ECMWF precipitation, temperature and evapotranspiration ensemble forecasts on a daily basis. The focus was on clusters where seasonal forecasting was difficult. They concluded that the LS and QM techniques were able to remove the mean bias. Regarding the estimation of dry days and low precipitation amounts, the efficiency of QM was better than that of LS. Mishra et al.(2019) provided one of the most comprehensive assessments of seasonal temperature and precipitation ensemble forecasts of the EUROSIP multi-model forecasting system.One equally and two unequally weighted multi-models were also constructed from individual models, for both climate variables, and their respective forecasts were also assessed.They found that the simple equally weighted multi-model system performed better than both unequally weighted multimodel combination systems. Zhang et al. (2019) evaluated the ability of the seasonal temperature forecasts of the NMME over west coast areas of the United States. In general, the skill of the one-month lead time NMME forecasts was superior or similar to persistence forecasts over many continental clusters, while the skill was generally stronger over the ocean than over the continent. However, the forecast skill along most west coast clusters was markedly lower than in the adjacent ocean and interior, especially during the warm seasons.

    To the best of our knowledge, no research has yet been reported in which the precipitation forecasts of the C3S database have been comprehensively evaluated. In this study,the ECMWF, MF (Météo France) and UKMO ensemble precipitation forecasts were extracted from the C3S database in Iran’s geographical area for a period of approximately 24 years (1993-2017). The forecasts were compared with station data in different precipitation clusters. Evaluation of raw forecasts was performed in two stages: deterministic and probabilistic assessment.

    2. Methods and study area

    The C3S project was introduced in early 2017 and has since been routinely releasing forecast products. These products, which are taken from several European centers,are updated monthly (13th day of the lunar month at 1200 UTC) for up to a six-month lead time. Lead time refers to the period of time between the issue time of the forecast and the beginning of the forecast validity period. Long-range forecasts based on all data up to the beginning of the forecast validity period are said to be of lead zero. The period of time between the issue time and the beginning of the validity period will categorize the lead. For example, a March monthly forecast issued at the end of the preceding January is said to be of one-month lead time. The C3S climate data store is currently supported by the ECMWF, UKMO, MF,Centro Euro-Mediterraneo sui Cambiamenti Climatic, and Deutscher Wetterdienst centers. The NCEP, JMA and Bureau of Meteorology centers will be added in the near future.In the data store, a global-scale meteorological observed dataset has been used to obtain hindcasts since 1993 and may be used to improve forecast quality. The characteristics of the forecasting systems and their production methods in the C3S are presented in Table 1.

    In this study, the ECMWF, MF and UKMO monthly ensemble precipitation hindcasts of the C3S database in Iran’s geographical area (25°-40°N, 44°-64°E) at an approximate 25-km spatial resolution (0.25° × 0.25°) were extracted for a period of 24 years (1993-2017) at a three-month lead time with a monthly time step. Observed point data were extracted for 100 synoptic stations, operated by the Iranian Meteorological Organization, spread over eight different precipitation clusters as classified by Modarres (2006)based on the geography and climate using the Ward method. Figure 1 shows a map of Iran’s precipitation clusters, over which the locations of 100 stations used in this study are overlaid. Cluster G1 involves arid and semiarid clusters of central and eastern Iran, subject to high coefficient of variation and low precipitation. Cluster G2, spread in three distinct clusters (as shown in Fig. 1), mostly encompasses highland margins of G1. Cluster G3 involves cold clusters in northwestern Iran, while cluster G4 represents warm and temperate clusters along the Persian Gulf northern coast. Cluster G5 involves areas over the western border along the Zagros Mountains. Cluster G6, enjoying high precipitation, represents areas along the coast of the Caspian Sea in northern Iran. Cluster G7 is similar to G5 but receives more precipitation and encompasses two distinct clusters. In addition, cluster G8 is similar to G6, but receives more precipitation. Table 2 lists the number of stations and characteristics of the precipitation clusters.

    Using the inverse distance weighting method (Ozelkan et al., 2015), the forecast data were interpolated to the selected stations, and then monthly precipitation forecasts were compared with observed data. Evaluation criteria were employed in two stages, including deterministic and probabilistic assessment. It should be noted that the averages of the evaluation criteria in each precipitation cluster are reported in this study. All the evaluation criteria formulae used in this study are presented in Table 3.

    2.1. Evaluation criteria

    For deterministic evaluation, the Pearson correlation coefficients were used to compare the forecast values with those of the observed. The root-mean-square error skill score(RMSE) was also adopted, to calculate the error intensity of the forecasted monthly precipitation. Since each cluster has different mean precipitation, the skill score of this criterion was used for a fair assessment. The RMSEis 1 for a perfect forecast, and 0 when the forecast equals the climatology. It should be noted that in order to prevent dispersion,anomaly forecasts and observations were used for calculating the RMSE. A monthly anomaly is the actual monthly value minus the climatology of the same month (i.e., the average over 1993-2017).

    Table 1. Details of the selected numerical forecast models.

    For the probabilistic evaluation, the continuously ranked probability score (CRPS) and the relative operating characteristic (ROC) were used. To calculate the ROC, for each category (low/middle/upper tercile) a binary observation (0: the category occurred, 1: the category was not observed) and a probabilistic (between 0 and 1) forecast were created. The latter is derived based on the number of members predicting the category, out of the total number of available members. Then, for each tercile category, based on the values of hit and false alarm rates that were obtained based on a probability threshold varying from 0 to 1, the ROC curve was drawn and the area under the curve was calculated. Finally, the ROC skill score (ROC) was computed as 2

    A

    - 1, where

    A

    is the area under the ROC curve. The ROCrange is -1 to 1, where zero indicates no skill when compared to the climatological forecasts and 1 represents the perfect score (Manzanas et al., 2014). Moreover, using the CRPS, the agreement of the cumulative distribution function (CDF) forecasts with the observed CDF was studied.The perfect score for the CRPS criterion is 0. The CRPS criterion was calculated based on Ferro et al. (2008). As with the reason for using the RMSEss, the CRPS skill score(CRPSss) was used for a fair assessment between clusters.The CRPSss was determined based on the relationship presented in Table 3.

    The evaluation results were interpreted on a monthly,three-month average, and annual basis. On the monthly basis, the performance of each model was individually investigated in all 12 months; whereas for the 3-month average,the average precipitation in all four seasons was determined.Then, forecasts were evaluated in three (lower, middle and upper) tercile categories.

    Fig. 1. Layout of the observation stations and G1-G8 precipitation clusters [the boundaries of the eight groups are based on the research of Modarres (2006)].

    2.2. Multi-model forecasts

    Three selected models were combined to generate multi-model forecasts using the simple arithmetic mean(MMM), multivariate regression (MRMM), and biasremoved (BRMM) techniques. In the MMM, simply the mean of the three model forecasts was calculated in each time step.

    In the MRMM combination, a multivariate regression was developed to combine the forecasts of the three models based on the following relationship (Zhi et al., 2012):

    In the BRMM approach, a multi-model was generated based on the following relationship (Zhi et al., 2012):

    Table 2. Details of the precipitation clusters and the number of selected stations in each cluster (Kolachian and Saghafian, 2019).

    The correlation coefficient and relative root-meansquare error (RRMSE) were used to evaluate and compare the skills of individual forecast models with the constructed multi-models. These criteria were selected based on a similar study by Zhi et al. (2012). Since precipitation varies greatly among clusters, the RMSEs of each cluster were divided by the mean observed value of the same cluster to make a fair assessment. It should be noted that the evaluations in this part of the study were carried out separately for the four seasons.

    Table 3. Evaluation criteria used in this study (Tao et al., 2014; Aminyavari et al., 2018).

    2.3. Evaluation based on drought indices

    The Standardized Precipitation Index (SPI) was also adopted, to evaluate the skill of the models to forecast dry/wet years. The SPI represents the number of standard deviations (SDs) that observed cumulative precipitation deviates from the climatological average (Guttman, 1999). To calculate the SPI, the gamma distribution was fitted to 30 years(1987-2017) of monthly precipitation series and then converted into the standard normal distribution. The mean SPI of the stations in each of the eight precipitation clusters were determined and intercompared. The SPEI package(Beguería and Vicente-Serrano, 2017) was employed to calculate the SPI in R software.

    Further, to study the effect of climatic signals on precipitation changes in the study period, the Oceanic Ni?o Index(ONI) was extracted (https://ggweather.com/enso/oni.htm)and then El Ni?o, La Ni?a and neutral years (phases) were identified over the study period. Furthermore, the relationship between the ENSO phases and the SPI was examined and interpreted. It should be noted that events were defined as five consecutive overlapping three-month periods at or above the +0.5°C anomaly for El Ni?o events and at or below the -0.5°C anomaly for La Ni?a events. The threshold was further divided into weak for 0.5°C to 0.9°C SST anomalies, moderate (1.0°C to 1.4°C), strong (1.5°C to 1.9°C), and very strong (≥ 2.0°C) events (Rojas, 2020).

    Also, regarding the effect that ENSO may impose on the skill of different models, tercile plots (Nikulin et al.,2018) were used to examine the sensitivity of the skill to El Ni?o, La Ni?a, and neutral years. The tercile plots were drawn using the visualizeR package in R software (Frías et al., 2018).

    3. Results

    The average monthly forecasted precipitation over the 24 years of the study period with a one-month lead time is compared in Fig. 2 with the average monthly recorded precipitation at synoptic stations. One may note that the MF model yielded poor performance in most clusters and seasons. However, the ECMWF and UKMO models had acceptable performance, particularly in low-precipitation and dry months, as well as in low-precipitation clusters (such as G1). For G4 and G1, which together cover the greatest area of Iran, ECMWF and UKMO underestimated precipitation in high-precipitation months, but overestimated precipitation in G2, G3 and G5, although the UKMO performed somewhat better than the ECMWF in highland snow-covered areas. For G6 and G8, the models underestimated precipitation considerably in fall, but provided better estimates than the reference forecasts in other seasons. All in all, in G5 and G7, the UKMO and ECMWF models forecasted the precipitation closer to observations.

    In the following, first, the total precipitation forecasted during the study period will be evaluated. Then, the seasonal precipitation is examined in different clusters and terciles, and the performance of the models in different months is also interpreted. Then, the skills of multiple models are compared with individual models, and finally, based on the drought index, the quality of forecasts is evaluated.

    Fig. 2. Average monthly forecasted precipitation values of ECMWF, UKMO and MF, along with observed data, in precipitation clusters with a one-month lead time.

    3.1. Evaluation of total forecasts

    The Pearson correlation coefficient results in Fig. 3 show that ECMWF and UKMO performed best in the G5 and G7 clusters, with over 70% significant correlation with observations. But, the correlation coefficients of MF were below 0.5 in most precipitation clusters. In the G6 and G8 precipitation clusters, encompassing the north of Iran, all models had correlation coefficients of less than 0.5, which were lower than in other clusters. Moreover, the skill of the models decreased with an increase in lead time. Overall, based on this criterion, the ECMWF model was in better agreement with the observations compared with the other two models, although UKMO was slightly better than ECMWF for a three-month lead time. Moreover, the SD of the correlation coefficients of the stations in each cluster is shown to understand the degree of dispersion in the evaluation result. As seen in the figure, in most clusters the SD is low and models perform the same at most stations. Only in clusters 6 and 8 is the SD slightly higher, which also affects the evaluation results of these two clusters.

    Pearson’s product-moment correlation statistical significance test was also performed, to determine whether the degree of correlation between models and observations was reliable. The

    P

    -values of the significance test, as shown in the table within Fig. 3a, indicates that the

    P

    -value in all models and clusters is less than 0.05, implying statistical significance of the correlation between the forecast models and observations. Only for the MF model in G6 cluster is the

    P

    -value slightly greater than 0.05.

    The evaluation results of the RMSE skill score in Fig. 4 indicate that the models were more skillful in western Iran.However, in G3, encompassing northwestern Iran, the ECMWF model had a negative skill score and performed poorly in comparison with the reference forecasts. Similarly, UKMO provided poor forecasts in G3. The model also produced poor forecasts for the high-precipitation G8 cluster.

    The CRPSss results in Fig. 5 demonstrate that, in G4,G5 and G7, the models had closer CDFs to those of the observations and performed quite well in the clusters located in western Iran.

    In general, the annual evaluation results indicated that the models provided better forecasts in western and southwestern Iran. In contrast, the models performed quite poorly in precipitation clusters in northern Iran. The ECMWF model outperformed the UKMO model in most evaluation criteria.The MF forecasts were poor. The gray cells in Figs. 4 and 5 are differentiated because their values are far from the other scores.

    In Figs. 4, 5 and 7, the numbers in the squares correspond to the average score for the entire cluster (i.e., averaged over all stations within the cluster), while the small sized numbers indicate the SDs of the scores across all the stations in the cluster.

    Also, to understand why the forecasts were better in G1, G2, G4, G5 and G7, as compared with those in G3, G6 and G8, the time series of total annual precipitation is shown in Fig. 6. The results of the evaluation are fully consistent with Fig. 6, showing why the models have varying evaluation scores in different clusters.

    Fig. 3. Evaluation using the (a) mean and (b) SD of the Pearson correlation coefficient in each cluster. The table in panel (a) shows the P-values of the three models in the eight clusters (perfect score for the Pearson correlation coefficient is 1).

    3.2. Evaluation of three-month average precipitation

    In this section, the evaluation scores are reported based on three-month average observations and hindcasts at a onemonth lead time. Based on Fig. 7, both the UKMO and ECMWF models yielded positive skill scores in most clusters (with the exception of northern climates in the G3,G6 and G8 clusters), indicating that the models forecasted the probability of occurrence better than the climatology.All three models had better proficiency scores in higher and lower terciles than in the middle tercile. In the summer,according to Fig. 7, the performance of all models declined compared to those in other seasons, whereas only the UKMO model scored higher in the G7 cluster.

    For the fall season, the ECMWF model forecasted light precipitation better than heavy precipitation. In this season,the MF model yielded a more accurate detection and lower false alarm rate than in other seasons. Overall, similar to the findings of Manzanas et al. (2014), all models produced better forecasts in this season than in the other seasons. In the G5 and G7 precipitation clusters, all three models performed better compared to the other seasons and precipitation clusters. In a similar study by Shirvani and Landman(2016), the best performance was found in the fall season.In winter, the performance of models was similar to that in the spring season, but with lower skill scores.

    Fig. 4. Evaluation using the RMSESS in eight precipitation clusters (RMSESS is 1 for “perfect” forecast and 0 for“forecast equals climatology”).

    3.3. Evaluation of monthly precipitation

    In this section, the performance of the models was investigated on a monthly basis in eight precipitation clusters.According to the Pearson correlation coefficients shown in Fig. 8, it is clear that the correlation coefficient decreased with increasing lead time. There was also no specific dependence between the lead time and the month of the year. The models produced good forecasts in spring months in dry climates, and in December in wet climates at two- and threemonth lead times.

    The ECMWF model performed best in November,while it had its poorest performance in August in most clusters. For western clusters in Iran, such as G7, which receive greater precipitation compared to central clusters,the model performed better in summer months than in other seasons; however, poorer forecasts were achieved in February. As mentioned in section 3.1, ECMWF did not provide acceptable performance in northern clusters. The UKMO model was to a certain extent similar in its performance to the ECMWF model, except that it performed best in dry central clusters in the winter. The UKMO model performed quite poorly in northern clusters in May. The MF forecasts did not correlate well with the observations in all months.For most clusters, this model provided better forecasts in March and April.

    Also, to better assess the forecasts in different months,scatterplots of monthly forecasts over the study period are provided in Fig. 9. Based on the results of this study, the models had their best forecasts in G7 but their weakest ones in the G8 precipitation cluster. For the sake of brevity, scatterplots of only these two clusters are presented. According to Fig. 9a, in G7, observed precipitation was zero in the summer months, while the highest correlation occurred in January between the observations and ECMWF forecasts. From June to September, observed precipitation in this cluster was close to zero, so the highest ROCin section 3.2 was obtained in this period. This highlights that the ROCis conditioned on the observation such that the models accurately forecast dry conditions over that cluster. Also, in rainy months, the ECMWF and UKMO models performed well,and the scores shown in Fig. 9 may be attributable to the forecast skill.

    In the G8 precipitation cluster, most models underestimated the observations, with only the ECMWF model overestimating precipitation in the summer months. Although more precipitation occurred in winter than in spring, all models forecasted similar precipitation in these two seasons. All three models underestimated the fall and winter precipitation. The ECMWF model was a better performer than the other models.

    Fig. 5. Evaluation using CRPSSS in eight precipitation clusters (perfect score for CRPSSS is 1).

    3.4. Evaluation of multi-models

    According to Fig. 10a, multi-models, especially MRMM, were better correlated with observations in most seasons and clusters than individual models, although MRMM had weaker forecasts in summer. The MF model in winter,especially in rainier clusters, showed no correlation with observations. Overall, the best forecasts for the multi- and individual models occurred in fall.

    Figure 10b shows the evaluation results of multi- and individual models using the RRMSE criterion. Of note is a large error in the MF forecasts. Despite the poor performance of the MF model, the three multi-models tended to provide similar results to those of the ECWMF and UKMO models, which represent the two best-performing individual models. Note that among the three multi-models, the MRMM provided the best overall results. The MRMM results showed that the model combination based on regression had a positive effect on reducing the forecast error of individual models. In accordance with previous studies(Doblas-Reyes et al., 2009; Bundel et al., 2011; Ma et al.,2012; Manzanas et al., 2014), the overall results in this section indicate that multi-models are effective in improving the predictive skill of individual models.

    3.5. Evaluation of forecasted SPI

    Figure 11 shows the time series of observed SPI and those corresponding to the ECMWF and UKMO model forecasts in eight precipitation clusters. The MF-derived SPI values are not shown here because of its poor performance. El Ni?o, neutral and La Ni?a years were extracted during the study period and added to Fig. 11. To better understand the performance of models in forecasting dry and wet years, the correlation coefficient between the SPI of the models and that of the observation is shown in Fig. 11. In G1, both models forecasted the wet years during 1995-2002 well, albeit with a slight overestimation by the UKMO model. Furthermore, both models forecasted the observed 2000-2002 dry period, although they extended the drought until 2004. In this cluster, the models forecasted a severe wet year in 1996, but the year was neutral based on the ONI. This was also the case for the 2001 severe drought, which was a La Ni?a year based on the ONI. Overall, the two models performed quite well in this cluster. In G2, the results were similar but weaker to those of the G1 cluster, whereas the UKMO model had poorer forecasts than the ECMWF model.

    Fig. 6. Total annaul precipitation values of the ECMWF, UKMO and MF models, along with observed data, in precipitation clusters with a one-month lead time.

    In G3, the UKMO model forecasted the severity of the SPI in wet and dry periods better than the ECMWF model.It should be noted that, while 1994 and 2008 were severe wet and drought years in the G3 precipitation cluster, respectively, El Ni?o and La Ni?a prevailed in 1994 and 2008.Thus, climate indices are not a good sign of wet/dry conditions in G3 in such years. In G4, ECMWF forecasted the drought/wetness index better than the UKMO model,although both models underestimated the dry periods in 2008-2012. Considering G5, ECMWF overestimated the drought from 2000 to 2003, while both models failed to forecast the wet year in 2016. In G6, the UKMO model slightly overestimated the dry periods. Both models forecasted the wet years from 1993 to 2000. However, the forecasted duration and intensity of wet years were shorter and weaker,respectively, than those observed.

    Regarding G7, both models performed well and ECMWF forecasted the wet and dry periods quite close to the observations. Finally, in G8, both models had average performances and failed to forecast dry and wet periods. Overall, the models performed well in forecasting dry/wet periods in most clusters and, nearly in all clusters, they were able to detect the reduction in precipitation over the 2000-2001 period.

    Fig. 7. Seasonal evaluation using ROCSS in eight precipitation clusters for three (L=lower, M=middle and U=upper) tercile categories (ROCss is 1 for “perfect” forecast and 0 for “forecast equals climatology”): (a) spring (MAM, March-April); (b)summer (JJA, June-August); (c) fall (SON, September-November); (d) winter (DJF, December-February).

    Fig. 8. Monthly evaluation using the Pearson correlation coefficients in eight precipitation clusters (perfect score for Pearson correlation is 1).

    Fig. 9. Scatterplots of monthly precipitation in two precipitation clusters: (a) G7; (b) G8.

    In examining the effects of ENSO in the study period,only heavy precipitation in 1995 could be marginally attributed to this phenomenon. Based on the work of Shirvani and Landman (2016), although ENSO is the main factor in seasonal forecast skill (Manzanas et al., 2014), no clearly strong predictive capability was found in this study for Iran.This may be partly attributable to the complicated precipitation variability in relation to SST patterns.

    Finally, tercile plots were used to evaluate the impact of ENSO on the performance of models. For the sake of brevity, only the results of the G5 cluster, which had the best forecast results based on Fig. 7, are shown. According to Fig. 12a,although no clear relationship could be found between model performance and ENSO conditions, the ECMWF members of the ensemble forecast in El Ni?o years were slightly better than those in La Ni?a and neutral years. Moreover,according to Fig. 12b, the UKMO model did slightly better in the La Ni?a years. All in all, there is no clear connection between ENSO phenomena and precipitation predictability at seasonal time scales in Iran.

    Fig. 10. Individual and multi-model seasonal evaluation using the Pearson correlation coefficient and RRMSE in eight precipitation clusters (perfect score for Pearson correlation is 1, and for RRMSE it is 1). MAM, March-May(spring); JJA, June-August (summer); SON, September-November (fall); DJF, December-February (winter).

    Fig. 11. Observed and forecasted SPI time series in G1 to G8 (the small 2 × 2 tables shown in the top-right corner of each panel indicate the correlation between the model-forecasted SPI and those of observations. N,neutral; W, weak; M, medium; S, strong; VS, very strong; E, El Ni?o; L, La Ni?a.

    Fig. 12. Tercile plots for the (a) ECMWF and (b) UKMO models in the G5 cluster. The red color spectrum of each square represents the probability for each category (below, normal, above). Dots show the corresponding observed category for each year of the hindcast period. Numbers on the right show the ROCSS for each model and each tercile.(ROCSS is 1 for “perfect” forecast and 0 for “forecast equals climatology”).

    4. Conclusions

    In this study, ECMWF (SEAS5), UKMO (GloSea5)and Météo-France (System5) monthly precipitation forecasts within the C3S database over the period 1993-2017 were evaluated in eight precipitation clusters in Iran. The evaluations were performed on an annual, seasonal and monthly basis. The following conclusions may be drawn:

    (1) Based on the deterministic evaluation of the forecasts, the best correlation coefficients were achieved in the G5 and G7 clusters in western Iran, while the poorest performance was associated with G6 and G8 in northern Iran.Moreover, in forecasting the precipitation amount, all three models had their greatest error in G3 in northwestern Iran.The models performed well in the dry central and eastern clusters. The MF model had the greatest error among the models based on deterministic evaluation (Figs. 3 and 4).

    (2) All models better forecasted upper-tercile events in dry seasons and lower-tercile events in wet seasons, but gained their lowest skill score in the middle category (Fig. 7)(similar to Manzanas et al., 2014, 2018).

    (3) Based on probabilistic evaluation, all three models scored higher in the G5 and G7 precipitation clusters compared to those of other clusters (Fig. 5)

    (4) In all evaluations, the skill of the models decreased with increasing lead time.

    (5) In the monthly evaluation, the ECMWF model performed better in low-precipitation clusters in fall and in high-precipitation clusters in summer. On the contrary, the model performed poorly in northern clusters. The UKMO model yielded somewhat similar results, except the model performed better in low-precipitation clusters in the winter. In general, all three models overestimated precipitation in the summer (Fig. 8).

    (6) The MRMM multi-model had better skill than individual models (Fig. 10).

    (7) The forecast models had relatively good skill in forecasting dry and wet years, although they underestimated some dry years and overestimated some wet years (Fig. 11).

    (8) No specific relationship was found regarding the influence of the ENSO climatic signal on model performance.The models did not yield acceptable forecasts in northern clusters where higher precipitation and relatively lower temperature prevail.

    (9) In assessing the impact of climatic global signals on severe precipitation, the year 1995 may be considered to be significantly influenced by ONI phenomena, while no major effect was detected in other periods.

    All in all, the evaluation results demonstrated that both the UKMO and ECMWF models perform well in forecasting monthly precipitation in Iran, especially in western precipitation clusters. However, it is not possible to clearly indicate which of the two models performs better. In most precipitation clusters at short lead times, the ECMWF model was better correlated with the observations and its forecasts gained higher skill scores. On the contrary, at a three-month lead time, the UKMO model had higher correlation coefficients with the observations in most precipitation clusters. The MF model is not an appropriate precipitation forecast model for Iran. Furthermore, the models are generally able to forecast dry/wet occurrences in Iran.

    精品国产一区二区三区久久久樱花| 精品少妇内射三级| 亚洲国产欧美日韩在线播放 | 国产精品一区二区在线不卡| 日韩一区二区视频免费看| 国产av精品麻豆| 国内揄拍国产精品人妻在线| 精品少妇久久久久久888优播| 麻豆成人午夜福利视频| 亚洲国产精品国产精品| av福利片在线观看| 成人国产麻豆网| 日本爱情动作片www.在线观看| 亚洲精品第二区| 观看免费一级毛片| 在线看a的网站| 免费观看a级毛片全部| 日本猛色少妇xxxxx猛交久久| videossex国产| 蜜桃久久精品国产亚洲av| 欧美人与善性xxx| 免费人妻精品一区二区三区视频| 久久影院123| 亚洲专区国产一区二区| 久久精品国产亚洲av香蕉五月 | 大香蕉久久网| 日本wwww免费看| 日日夜夜操网爽| 国产无遮挡羞羞视频在线观看| 两个人免费观看高清视频| 精品人妻在线不人妻| 国产av精品麻豆| 亚洲七黄色美女视频| 久久精品久久久久久噜噜老黄| 色视频在线一区二区三区| 国产亚洲一区二区精品| 十八禁网站网址无遮挡| 亚洲色图 男人天堂 中文字幕| 久久ye,这里只有精品| 一级a爱视频在线免费观看| 在线观看免费视频网站a站| 亚洲国产精品一区二区三区在线| 国产成人啪精品午夜网站| 在线 av 中文字幕| 久热这里只有精品99| 大型av网站在线播放| 亚洲第一欧美日韩一区二区三区 | 国产亚洲av高清不卡| av视频免费观看在线观看| 欧美在线一区亚洲| av一本久久久久| 日韩 亚洲 欧美在线| 精品国产一区二区三区四区第35| 男女之事视频高清在线观看| 欧美中文综合在线视频| 九色亚洲精品在线播放| 一本一本久久a久久精品综合妖精| 这个男人来自地球电影免费观看| 精品一区二区三区av网在线观看 | 日韩欧美国产一区二区入口| 欧美日韩亚洲综合一区二区三区_| 国产一区二区在线观看av| 亚洲av成人一区二区三| 精品人妻1区二区| 免费高清在线观看日韩| 两人在一起打扑克的视频| 午夜日韩欧美国产| 新久久久久国产一级毛片| 大陆偷拍与自拍| 99久久精品国产亚洲精品| 99久久99久久久精品蜜桃| 日韩欧美国产一区二区入口| 满18在线观看网站| 777米奇影视久久| 久久热在线av| 国产成人a∨麻豆精品| 亚洲一区中文字幕在线| 建设人人有责人人尽责人人享有的| 国产野战对白在线观看| 国产精品香港三级国产av潘金莲| 久久亚洲国产成人精品v| 黄色视频,在线免费观看| a级毛片在线看网站| 国产一区二区在线观看av| 淫妇啪啪啪对白视频 | 十八禁高潮呻吟视频| 日本欧美视频一区| 国产亚洲av片在线观看秒播厂| 午夜视频精品福利| 日韩 亚洲 欧美在线| 悠悠久久av| 久久久精品免费免费高清| 黑丝袜美女国产一区| 久久精品成人免费网站| 色婷婷av一区二区三区视频| 欧美大码av| 超碰成人久久| 又紧又爽又黄一区二区| 日韩欧美国产一区二区入口| 国产欧美日韩一区二区三 | 精品欧美一区二区三区在线| 午夜两性在线视频| 性高湖久久久久久久久免费观看| 在线av久久热| 亚洲五月色婷婷综合| 99精品欧美一区二区三区四区| 欧美日韩视频精品一区| 午夜福利影视在线免费观看| 18禁国产床啪视频网站| 别揉我奶头~嗯~啊~动态视频 | 亚洲七黄色美女视频| 日本猛色少妇xxxxx猛交久久| 侵犯人妻中文字幕一二三四区| 国产av一区二区精品久久| 女性生殖器流出的白浆| 母亲3免费完整高清在线观看| 18禁黄网站禁片午夜丰满| 亚洲综合色网址| 窝窝影院91人妻| 精品少妇一区二区三区视频日本电影| 国产精品久久久av美女十八| 国产精品.久久久| 如日韩欧美国产精品一区二区三区| 国产欧美日韩综合在线一区二区| 国产亚洲精品一区二区www | 日本91视频免费播放| 亚洲av电影在线进入| 日韩一卡2卡3卡4卡2021年| 9191精品国产免费久久| 亚洲少妇的诱惑av| 国产精品自产拍在线观看55亚洲 | 丁香六月天网| 欧美日本中文国产一区发布| 交换朋友夫妻互换小说| 极品人妻少妇av视频| 色视频在线一区二区三区| 97人妻天天添夜夜摸| av在线app专区| 国产亚洲欧美在线一区二区| 青草久久国产| √禁漫天堂资源中文www| 99国产精品一区二区三区| 十八禁高潮呻吟视频| 久久久国产精品麻豆| 久久九九热精品免费| 久久人人爽av亚洲精品天堂| 999精品在线视频| 视频区欧美日本亚洲| 别揉我奶头~嗯~啊~动态视频 | 精品亚洲乱码少妇综合久久| 日韩视频在线欧美| 国产在线免费精品| 51午夜福利影视在线观看| 麻豆av在线久日| 18禁裸乳无遮挡动漫免费视频| 窝窝影院91人妻| 国产亚洲一区二区精品| av视频免费观看在线观看| 成年美女黄网站色视频大全免费| 国产国语露脸激情在线看| xxxhd国产人妻xxx| 亚洲熟女毛片儿| 欧美黑人精品巨大| 在线看a的网站| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品久久久久成人av| 国产成人欧美| 最近中文字幕2019免费版| √禁漫天堂资源中文www| 国产精品久久久久久精品古装| 亚洲av电影在线进入| 国产精品免费大片| 亚洲伊人色综图| 亚洲国产欧美一区二区综合| 精品卡一卡二卡四卡免费| 久久亚洲精品不卡| 精品一区二区三区四区五区乱码| 国产精品国产三级国产专区5o| 成人亚洲精品一区在线观看| svipshipincom国产片| 久久久久久久大尺度免费视频| 精品国产一区二区三区久久久樱花| 黑人猛操日本美女一级片| 飞空精品影院首页| 亚洲七黄色美女视频| 黄片播放在线免费| 美女中出高潮动态图| 香蕉国产在线看| 美国免费a级毛片| 国产免费av片在线观看野外av| 女警被强在线播放| 黄色视频在线播放观看不卡| 天天添夜夜摸| 久久人人爽av亚洲精品天堂| 精品视频人人做人人爽| 免费不卡黄色视频| 婷婷色av中文字幕| 狠狠狠狠99中文字幕| 两个人看的免费小视频| avwww免费| 久久香蕉激情| 97精品久久久久久久久久精品| 女性生殖器流出的白浆| 黑人欧美特级aaaaaa片| 亚洲国产毛片av蜜桃av| 久久精品久久久久久噜噜老黄| 女性被躁到高潮视频| 精品亚洲成a人片在线观看| 天堂中文最新版在线下载| 免费久久久久久久精品成人欧美视频| 777久久人妻少妇嫩草av网站| 99re6热这里在线精品视频| 18禁裸乳无遮挡动漫免费视频| 久久久精品区二区三区| 啪啪无遮挡十八禁网站| 国产成+人综合+亚洲专区| 免费在线观看完整版高清| av片东京热男人的天堂| 一区二区三区激情视频| 狠狠婷婷综合久久久久久88av| 青青草视频在线视频观看| 无遮挡黄片免费观看| 青草久久国产| 久久国产精品影院| 最近最新中文字幕大全免费视频| 中文字幕高清在线视频| 高清视频免费观看一区二区| 午夜成年电影在线免费观看| av天堂在线播放| 在线观看免费日韩欧美大片| 午夜激情av网站| 亚洲av电影在线进入| 成人影院久久| 男女无遮挡免费网站观看| 日韩有码中文字幕| 韩国高清视频一区二区三区| 肉色欧美久久久久久久蜜桃| 国产成人免费无遮挡视频| 午夜免费观看性视频| 精品国产一区二区久久| 久久久久久人人人人人| 母亲3免费完整高清在线观看| 国产成人精品久久二区二区免费| 精品少妇黑人巨大在线播放| 久久午夜综合久久蜜桃| 久久精品aⅴ一区二区三区四区| 午夜老司机福利片| 悠悠久久av| 久久免费观看电影| 搡老熟女国产l中国老女人| 国产精品二区激情视频| 老司机福利观看| 嫩草影视91久久| 一区二区三区精品91| 超碰成人久久| 亚洲九九香蕉| 黄色a级毛片大全视频| 日韩精品免费视频一区二区三区| 成人国产av品久久久| 日本91视频免费播放| 亚洲熟女精品中文字幕| 日韩中文字幕欧美一区二区| 狠狠狠狠99中文字幕| 中文字幕高清在线视频| 成年人免费黄色播放视频| 法律面前人人平等表现在哪些方面 | 黄色a级毛片大全视频| 国产成人系列免费观看| 免费在线观看黄色视频的| 久久久久国产精品人妻一区二区| 永久免费av网站大全| 色播在线永久视频| 亚洲男人天堂网一区| 自拍欧美九色日韩亚洲蝌蚪91| 欧美+亚洲+日韩+国产| 汤姆久久久久久久影院中文字幕| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲第一欧美日韩一区二区三区 | 巨乳人妻的诱惑在线观看| 69精品国产乱码久久久| 国产精品一二三区在线看| 欧美xxⅹ黑人| 国产亚洲精品一区二区www | 日韩,欧美,国产一区二区三区| av网站在线播放免费| 亚洲精品一区蜜桃| 别揉我奶头~嗯~啊~动态视频 | 一进一出抽搐动态| 热re99久久精品国产66热6| 母亲3免费完整高清在线观看| 一本综合久久免费| 国产男女超爽视频在线观看| 大片免费播放器 马上看| 99re6热这里在线精品视频| 丰满人妻熟妇乱又伦精品不卡| 国产精品久久久av美女十八| 各种免费的搞黄视频| 嫩草影视91久久| 国产一区二区三区av在线| 9热在线视频观看99| 日韩熟女老妇一区二区性免费视频| 成人18禁高潮啪啪吃奶动态图| 精品一品国产午夜福利视频| 高清av免费在线| 精品亚洲成a人片在线观看| 国产极品粉嫩免费观看在线| 12—13女人毛片做爰片一| 最近最新免费中文字幕在线| 777久久人妻少妇嫩草av网站| 性色av乱码一区二区三区2| 中文字幕高清在线视频| 一个人免费在线观看的高清视频 | 久久久欧美国产精品| 高清av免费在线| 少妇精品久久久久久久| 成人国产av品久久久| 午夜福利,免费看| 国产精品 欧美亚洲| 国产又爽黄色视频| 成人三级做爰电影| 成年动漫av网址| 久久久久国产一级毛片高清牌| 精品人妻一区二区三区麻豆| 国产97色在线日韩免费| 一个人免费看片子| 大码成人一级视频| 欧美97在线视频| 18禁裸乳无遮挡动漫免费视频| 亚洲精品一区蜜桃| 99热网站在线观看| 在线观看人妻少妇| 亚洲免费av在线视频| 亚洲av成人不卡在线观看播放网 | 午夜福利视频在线观看免费| 精品国产一区二区久久| 十八禁高潮呻吟视频| 国产成人一区二区三区免费视频网站| 99久久人妻综合| bbb黄色大片| 亚洲国产欧美日韩在线播放| 操出白浆在线播放| 丝瓜视频免费看黄片| 99热国产这里只有精品6| 亚洲国产欧美日韩在线播放| 男人操女人黄网站| 久久久久久久大尺度免费视频| 后天国语完整版免费观看| 亚洲专区国产一区二区| 亚洲人成电影观看| 秋霞在线观看毛片| 亚洲精品日韩在线中文字幕| 秋霞在线观看毛片| 欧美精品一区二区免费开放| 12—13女人毛片做爰片一| 男人操女人黄网站| 亚洲专区中文字幕在线| 国产麻豆69| 99久久人妻综合| 亚洲国产av新网站| 91成人精品电影| 十分钟在线观看高清视频www| 十八禁高潮呻吟视频| 99热网站在线观看| 欧美av亚洲av综合av国产av| 欧美日韩亚洲国产一区二区在线观看 | 法律面前人人平等表现在哪些方面 | 一进一出抽搐动态| 精品人妻1区二区| 日韩一卡2卡3卡4卡2021年| 亚洲精品一区蜜桃| 少妇的丰满在线观看| 亚洲成人免费av在线播放| 日本vs欧美在线观看视频| 黑人巨大精品欧美一区二区蜜桃| 99热网站在线观看| 国产又爽黄色视频| 免费av中文字幕在线| 午夜免费观看性视频| 久久久精品94久久精品| 欧美黄色片欧美黄色片| av国产精品久久久久影院| 老汉色av国产亚洲站长工具| 午夜福利视频精品| 亚洲九九香蕉| 久久久久网色| 又大又爽又粗| 另类精品久久| 热99re8久久精品国产| 亚洲精品一区蜜桃| 91麻豆av在线| 久久久欧美国产精品| 91大片在线观看| 久久中文看片网| 亚洲国产精品一区三区| 欧美日本中文国产一区发布| 午夜日韩欧美国产| 日本a在线网址| 岛国在线观看网站| 亚洲av国产av综合av卡| 999精品在线视频| 国产欧美日韩综合在线一区二区| 精品亚洲成a人片在线观看| 狂野欧美激情性xxxx| 搡老岳熟女国产| 日本猛色少妇xxxxx猛交久久| 国产av一区二区精品久久| 久久午夜综合久久蜜桃| 日韩三级视频一区二区三区| netflix在线观看网站| 国产亚洲精品久久久久5区| 老司机亚洲免费影院| 国产成人精品久久二区二区91| 9热在线视频观看99| 久久久久国内视频| 蜜桃国产av成人99| 国产一区二区在线观看av| 国产亚洲av片在线观看秒播厂| 精品欧美一区二区三区在线| 国产不卡av网站在线观看| 淫妇啪啪啪对白视频 | 久久久久精品国产欧美久久久 | 91麻豆av在线| 国产一卡二卡三卡精品| 日日爽夜夜爽网站| 丝袜美腿诱惑在线| 日韩 亚洲 欧美在线| 久久性视频一级片| 自线自在国产av| 一本色道久久久久久精品综合| 欧美在线黄色| 丰满迷人的少妇在线观看| 啦啦啦 在线观看视频| 精品少妇内射三级| 高清黄色对白视频在线免费看| 免费高清在线观看日韩| 亚洲av片天天在线观看| 欧美97在线视频| 丰满人妻熟妇乱又伦精品不卡| 欧美激情 高清一区二区三区| 我的亚洲天堂| 亚洲中文av在线| 一区二区三区乱码不卡18| 男女下面插进去视频免费观看| 99热全是精品| 国产极品粉嫩免费观看在线| 国产精品久久久av美女十八| 成人影院久久| 性高湖久久久久久久久免费观看| 亚洲专区字幕在线| 电影成人av| 亚洲专区国产一区二区| 久久香蕉激情| tube8黄色片| 多毛熟女@视频| 日本精品一区二区三区蜜桃| 色播在线永久视频| 啪啪无遮挡十八禁网站| 婷婷色av中文字幕| 国产成人免费观看mmmm| 国产免费av片在线观看野外av| 亚洲欧美日韩高清在线视频 | 亚洲国产av影院在线观看| 天天影视国产精品| 欧美大码av| 国产精品九九99| 无限看片的www在线观看| 久久免费观看电影| 午夜精品国产一区二区电影| 亚洲伊人久久精品综合| 精品人妻在线不人妻| 高潮久久久久久久久久久不卡| 午夜福利在线免费观看网站| 曰老女人黄片| 亚洲天堂av无毛| 精品人妻一区二区三区麻豆| 欧美另类一区| 成人免费观看视频高清| 曰老女人黄片| 别揉我奶头~嗯~啊~动态视频 | 成人18禁高潮啪啪吃奶动态图| 国产一区二区激情短视频 | 亚洲五月婷婷丁香| 精品亚洲乱码少妇综合久久| 一级毛片电影观看| 国产一区二区激情短视频 | 制服诱惑二区| 无遮挡黄片免费观看| 国产高清videossex| 久久久水蜜桃国产精品网| 搡老熟女国产l中国老女人| 国产精品 国内视频| 蜜桃国产av成人99| 爱豆传媒免费全集在线观看| 午夜福利,免费看| 亚洲欧美一区二区三区久久| 免费看十八禁软件| 亚洲 国产 在线| 亚洲综合色网址| 中文字幕人妻丝袜一区二区| 欧美人与性动交α欧美软件| 免费人妻精品一区二区三区视频| 精品人妻一区二区三区麻豆| 中文字幕最新亚洲高清| 日韩视频一区二区在线观看| 又紧又爽又黄一区二区| 日韩制服骚丝袜av| 欧美变态另类bdsm刘玥| 少妇猛男粗大的猛烈进出视频| 国产福利在线免费观看视频| 飞空精品影院首页| 一级片'在线观看视频| 亚洲精品美女久久久久99蜜臀| 女人爽到高潮嗷嗷叫在线视频| 午夜福利在线观看吧| 中文字幕制服av| 高清av免费在线| 久久久水蜜桃国产精品网| 热99国产精品久久久久久7| 91国产中文字幕| 丁香六月天网| 久久综合国产亚洲精品| 99精国产麻豆久久婷婷| 国产日韩欧美视频二区| 99久久精品国产亚洲精品| 亚洲视频免费观看视频| 精品久久久精品久久久| 欧美日韩成人在线一区二区| 亚洲成人免费av在线播放| 精品第一国产精品| 免费女性裸体啪啪无遮挡网站| 一区二区三区激情视频| 一本色道久久久久久精品综合| 亚洲男人天堂网一区| 又黄又粗又硬又大视频| 如日韩欧美国产精品一区二区三区| 少妇精品久久久久久久| 精品国产乱子伦一区二区三区 | 啦啦啦啦在线视频资源| 日韩中文字幕视频在线看片| 99国产精品免费福利视频| 国产精品熟女久久久久浪| 国产熟女午夜一区二区三区| 日韩电影二区| 99国产精品99久久久久| 免费看十八禁软件| 国产亚洲精品第一综合不卡| 91老司机精品| 亚洲国产精品一区二区三区在线| 女人被躁到高潮嗷嗷叫费观| www.精华液| 国产成人精品在线电影| 国产日韩欧美视频二区| 久久久久视频综合| 成年人免费黄色播放视频| 手机成人av网站| 18禁黄网站禁片午夜丰满| 久久久精品免费免费高清| 王馨瑶露胸无遮挡在线观看| 黄片大片在线免费观看| 欧美av亚洲av综合av国产av| 脱女人内裤的视频| 亚洲国产欧美一区二区综合| 十八禁高潮呻吟视频| 国产一区二区 视频在线| 国产亚洲精品久久久久5区| 国产一级毛片在线| 国产成人啪精品午夜网站| 国产精品一区二区在线观看99| 麻豆乱淫一区二区| 日韩中文字幕视频在线看片| 黄色视频不卡| 久久国产精品影院| 一本色道久久久久久精品综合| 久久久久精品人妻al黑| 手机成人av网站| 亚洲国产av新网站| 亚洲av美国av| 一边摸一边抽搐一进一出视频| 女警被强在线播放| 国产一卡二卡三卡精品| 不卡av一区二区三区| av福利片在线| 久久久精品免费免费高清| 精品久久久精品久久久| 国产欧美日韩一区二区精品| 国产成人系列免费观看| 久久国产亚洲av麻豆专区| 国产黄色免费在线视频| 国产亚洲av片在线观看秒播厂| 日本av手机在线免费观看| 老鸭窝网址在线观看| 国产野战对白在线观看| 国产99久久九九免费精品| 母亲3免费完整高清在线观看| 亚洲少妇的诱惑av| 久久99热这里只频精品6学生| 美女大奶头黄色视频| 97在线人人人人妻| 91老司机精品| 母亲3免费完整高清在线观看| 极品少妇高潮喷水抽搐| 精品少妇久久久久久888优播| a级毛片在线看网站| 波多野结衣av一区二区av| 丁香六月欧美| 欧美日韩成人在线一区二区| 日韩 欧美 亚洲 中文字幕| 美女福利国产在线| 成人手机av| 久热这里只有精品99| 中文字幕精品免费在线观看视频| 国产欧美日韩综合在线一区二区| 1024香蕉在线观看| 亚洲va日本ⅴa欧美va伊人久久 | 超碰97精品在线观看|