999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

Interannual Relationship between the Winter Aleutian Low and Rainfall in the Following Summer in South China

2015-11-24 09:49:56SONGLinYeandDUANWanSuo

SONG Lin-Yeand DUAN Wan-Suo

1State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics, Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100029, China

2University of Chinese Academy of Sciences, Beijing 100049, China

Interannual Relationship between the Winter Aleutian Low and Rainfall in the Following Summer in South China

SONG Lin-Ye1,2and DUAN Wan-Suo1*

1State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics, Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100029, China

2University of Chinese Academy of Sciences, Beijing 100049, China

Using observations and reanalysis data, this study investigates the interannual relationship between the winter Aleutian Low (AL) and the rainfall anomalies in the following summer in South China (SC). Results show that the winter AL is significantly positively (negatively) correlated with the SC rainfall anomalies in the following July (August). Specifically, SC rainfall anomalies have a tendency to be positive (negative) in July (August) when the preceding winter AL is stronger than normal. The winter AL-related atmospheric circulation anomalies in the following summer are also examined. When the winter AL is stronger, there is a significant anticyclonic (cyclonic) circulation anomaly over the subtropical western North Pacific in the following July (August). Southerly (northerly) wind anomalies to the west of this anomalous anticyclonic (cyclonic) circulation increase (decrease) the northward moisture transportation and contribute to the positive (negative) rainfall anomalies over SC in July (August). This study indicates that the AL in the preceding winter can be used as a potential predictor of the rainfall anomalies in the following July and August over SC. Keywords: Aleutian Low, South China, rainfall, atmospheric circulation

1 Introduction

South China (SC), with its large population and rapid socioeconomic development, is one of the most prosperous regions in China. SC receives the bulk of its annual rainfall in the boreal summer season (e.g., Wu et al., 2012; Duan et al., 2013). The interannual and intraseasonal variations of summer rainfall in SC are very large, often bringing extreme floods and droughts that can exert significant impacts on agriculture, industry, and the daily lives of people (e.g., Huang et al., 2003). Hence, to skillfully forecast the summer rainfall in SC is of great importance for both society and economic development, urging a better understanding of the summer rainfall variability in SC, as well as its influencing factors.

It has been demonstrated in previous studies that thevariability of summer rainfall in SC is influenced by many factors, such as the El Ni?o-Southern Oscillation (ENSO) cycle (e.g., Wu et al., 2003; Wang et al., 2010), Indian Ocean and Atlantic Ocean sea surface temperature (SST) (e.g., Tang et al., 2008; Wu et al., 2009; Xie et al., 2009), Eurasian or Tibetan snow cover (e.g., Wu and Kirtman, 2007), the Arctic Oscillation (e.g., Gong et al., 2011; Chen et al., 2015), the circulation anomaly over the Ural Mountains (e.g., Mao et al., 2011), the East Asian westerly jet (e.g., Liang and Wang, 1998), the western North Pacific subtropical high (e.g., Tao et al., 2001), and the East Asian summer monsoon (e.g., Yang and Lau, 2006). The Aleutian Low (AL) is the first dominant mode of North Pacific atmosphere circulation variability (e.g., Overland et al., 1999). It is characterized by a semi-permanent low pressure center located near the Aleutian Islands during boreal winter. Previous studies have demonstrated that variability of the AL is an important indicator of winter climate change around its surrounding regions (e.g., Parkinson, 1990; Trenberth and Hurrell, 1994; Niebauer et al., 1999). For instance, Trenberth and Hurrell (1994) showed that the boreal winter AL is closely connected with the Pacific-North America teleconnection pattern. Niebauer et al. (1999) reported that the winter temperature anomaly in the Bering Sea tends to be positive (negative) when the strength of the AL is above (below) normal.

*Correspongding author: DUAN Wan-Suo, duanws@lasg.iap.ac.cn

Previous studies have also found that variability of the AL can influence the climate anomalies over remote regions (e.g., Zhu and Wang, 2010; Yu and Kim, 2011; Park et al., 2012; Xiao et al., 2014). For example, Xiao et al. (2014) suggested that the summer surface temperature over the Arctic is modulated by interannual variability of the preceding AL. Zhu and Wang (2010) detected a significant relationship between winter AL and the simultaneous Australian summer monsoon rainfall on interannual time scales.

As discussed above, the AL can exert significant influences on the climate anomalies in its surrounding regions and remote areas. However, it remains unexplored as to whether the AL has an influence on the summer rainfall over SC. If such an influence exists, it could potentially be a good precursory signal for seasonal rainfall prediction in SC. In addition, previous studies regarding the connection of the AL with climate anomalies focused mainly on the simultaneous boreal winter (austral sum-mer) climate. Few attempts have been made to study the delayed influence of the boreal winter AL on the climate of the following summer.

In the present study, we investigate the possible linkage between the boreal winter AL and the rainfall in the following summer over SC (hereafter referred to simply as‘following summer SC rainfall’). We find a significant connection between the intensity of the winter AL and the rainfall anomalies in July and August of the following year over SC. Hence, the objective of this study is to investigate the influence of the boreal winter AL on the rainfall in the following July and August. The remainder of the paper is organized as follows: The data used in the study are described in section 2. Section 3 examines the relationship between the AL and summer rainfall over SC, describes the atmospheric circulation, and discusses the possible mechanism by which the AL influences the rainfall. Section 4 concludes the study.

2 Data

This study uses a monthly high-resolution precipitation dataset over land from the University of Delaware (UD) for the period 1900-2010 (Legates and Willmott, 1990). This dataset has a horizontal resolution of 0.5° × 0.5°. Also employed are the 160-station rainfall data from the National Climate Center of the China Meteorological Administration (CMA) from 1951 to the present.

The monthly mean sea level pressure (SLP) dataset is from the HadSLP2 product of the UK Met Office Hadley Centre, which has a horizontal of 5° × 5° and is available for the period 1850-2004 (Allan and Ansell, 2006).

Horizontal winds and geopotential heights are extracted from the 20th Century Reanalysis, version 2 (20CR2), of the National Oceanic and Atmospheric Administration-Cooperative Institute for Research in Environmental Sciences, for the period 1871-2010 (Compo et al., 2011), with a horizontal resolution of 2° × 2°. The analysis in this study focuses on interannual variation. Hence, all data are subjected to a nine-year high-pass filter. The results shown in this analysis are not sensitive to a reasonable variation of the cutoff period from 5 to 11 years. Previous studies have found that boreal winter ENSO events could exert influences on the variability of the winter AL (e.g., Roeckner et al., 1996; Alexander et al., 2002), which may confuse the influence of the winter AL on the following summer SC rainfall. As ENSO possesses a phase-locking feature with a mature peak in boreal winter, we subtract the ENSO signal in the variables by means of linear regression with respect to the Ni?o3.4 index during December-January-February (DJF). The Ni?o3.4 index is defined as the area-averaged SST anomalies over the region (5°S-5°N, 170-120°W). The removal of signals associated with the DJF Ni?o3.4 index from the DJF AL index and analysis data can exclude the possibility that any significant correlation between the DJF AL and following summer SC rainfall is due to the influence of the DJF ENSO on both the DJF AL and following summer SC rainfall.

3 Results

3.1 Definition of the AL and SC rainfall

Following Yu and Kim (2011), an empirical orthogonal function (EOF) analysis is performed to obtain the AL index (ALI). The region chosen for the EOF analysis is (20-60°N, 120°E-80°W). Figure 1a shows the leading EOF mode of monthly SLP anomalies over the selected region for the period 1900-2004. The corresponding time series of the leading EOF mode is displayed in Fig. 1b and is defined as the ALI in this study. The first EOF mode explains about 36.5% of the total SLP variance and can be separated well from the other eigenvalues. The spatial pattern of the first EOF mode is characterized by a center of action of negative SLP anomalies near the Aleutian region (Fig. 1a), which is consistent with the result in Yu and Kim (2011).

In this study, the region of SC is selected as (22-30°N, 105-118°E). There are 29 China Meteorological Administration (CMA) stations over this region, all distributed to the south of the Yangtze River. The time series of SC rainfall is defined as the area-mean rainfall over the selected SC region for the gridded rainfall dataset. For the CMA station rainfall dataset, the time series of SC rainfall is defined as the 29-station-averaged rainfall.

3.2 AL-rainfall connection

Figure 1 (a) The first EOF mode of monthly sea level pressure (SLP) anomalies over the North Pacific during 1900-2004. (b) The corresponding principle component (PC) time series.

Table 1 Correlation coefficients of South China summer rainfall (JJA, June, July, and August, separately) with the preceding winter (December-January-February, DJF) Aleutian Low index. The rainfall dataset are from the University of Delaware (UD) grid dataset and China Meteorological Administration (CMA) station dataset. All indices have been subjected to a nine-year high-pass filter and the DJF ENSO signal has been linearly removed in constructing this table.

In this section we examine the possible connection between the AL and SC rainfall by computing their correlation coefficients. Table 1 shows the correlation coefficients between the winter (DJF-averaged) ALI and the following summer (June-July-August, JJA-averaged and three individual months) rainfall over SC, after removing the influence of ENSO signal. Interestingly, the winter ALI is significantly and positively correlated with July rainfall but significantly and negatively correlated with August rainfall over SC. The standardized time series of the ALI and SC rainfall in July and August are displayed in Fig. 2. The correlation coefficient of the ALI with July (August) SC rainfall from the UD rain dataset is 0.32 (-0.33) for the period 1901-2004, significant at the 95% confidence level based on the Student’s t-test (Table 1). This indicates that an anomalously strong AL in the preceding winter tends to be followed by flood (drought) conditions in July (August) over SC. Note that no significant relationship is found between the ALI and JJA-averaged or June rainfall (Table 1). The insignificant correlation in the JJA-averaged rainfall might be attributable to the remarkably different behavior in these three individual months. In addition, the above correlation is computed using the station rain dataset for the period 1951-2004. The interannual variations of the station rainfall and UD rainfall are generally consistent with each other (Fig. 2). The correlation coefficient of the winter ALI is 0.37 (-0.37) with July (August) rainfall for the period 1951-2004 using the station rain dataset, significant at the 95% confidence level—consistent with that using the UD rain dataset (Table 1).

Figure 2 Normalized time series of DJF Aleutian low index (ALI) (bars) and (a) July rainfall indices from the University of Delaware (UD) rainfall dataset (red lines) and from China Meteorological Administration (CMA) 160-station dataset (blue lines). (b) As in (a) but for August rainfall indices. DJF ALI is also shown in (b) for comparison. Note that all indices have been subjected to a nine-year high-pass filter and the DJF ENSO signal has been linearly removed in constructing this figure.

Previous studies regarding the summer SC rainfall have focused mainly on the JJA-averaged rainfall (e.g., Wu et al., 2003, 2012). Our present study finds that the correlation between the boreal winter ALI and the following JJA-averaged SC rainfall is statistically insignificant. This may be due to the fact that the correlation of the ALI with the following July (August) SC rainfall is significantly positive (negative). Therefore, the subseasonal features of summer rainfall should be taken into account when investigating the SC summer rainfall. For example, Wang et al. (2009) suggested that a separate prediction of the bi-monthly anomalies for May-June and July-August may help improve the seasonal prediction of summer rainfall over East Asia.

Figures 3a and 3b show the correlation patterns of the boreal winter ALI with SC rainfall in July and August, respectively. Very different features appear between the two periods. In July, significant and positive correlations are seen over the SC region, with the maximum value centered at about (25°N, 110°E) (Fig. 3a). In contrast, significant and negative correlations are observed over the SC region in August, with the center of maximum correlation located more northeastward compared to that in July (Fig. 3b). We also examine the correlations of the ALI in other preceding months with the July and August SC rainfall. Figures 3c and 3d display the lag-lead correlations of the three-month means of the ALI with the SC rainfall index in July and August, respectively. The ALI is from preceding autumn (September-October-November, SON) to simultaneous summer (June-July-August, JJA). The results clearly show that the correlations are significant at the 90% confidence level both in July and August from the preceding early to late winter, with the maximum value appearing in DJF (Figs. 3c and 3d). The above results indicate that SC rainfall anomalies tend to be positive (negative) in July (August) when the AL is stronger than normal in the preceding winter.

The relationship between the winter AL and following summer SC rainfall is unstable. Figure 3e (3d) shows the sliding correlation between the DJF AL index and July (August) SC rainfall index with a window of 21 years. The correlation between the DJF AL and July SC rainfall is positive and statistically significant at the 90% confidence level around the 1910s, 1920s, late 1930s, and after the late 1970s; while the correlation is weak during the 1940s to the late 1970s (Fig. 3e). In addition, significant negative correlation between the DJF AL and August SC rainfall can be observed during the 1940s to 1950s and after the early 1970s; while before the 1940s and during the late 1950s to 1960s, the connection between the DJF AL and August SC rainfall is weak (Fig. 3f). The above results suggest that the connection of boreal winter AL with July and August SC rainfall experienced significant interdecadal changes in the past. The possible reasons leading to interdecadal change in the connection between the DJF AL and following summer SC rainfall will bepursued in a future study. As the connection between the DJF AL and following summer SC rainfall is better after the late 1970s (Figs. 3e and 3f), in the following analysis, we focus on the period 1979-2004 to investigate the possible physical process underlying the influence of the DJF AL on the following summer SC rainfall. It should be pointed out that atmospheric circulation anomalies in the following summer over the subtropical western North Pacific associated with the DJF AL during 1979-2004 are in close agreement with the results during 1900-2004 (not shown).

Figure 3 Correlation of DJF ALI with (a) July and (b) August rainfall in the following year. Heavy (light) shading in (a, b) indicates correlation coefficients significant at the 95% (90%) confidence level. Contour intervals are 0.1 in (a, b), and the green box is the region of South China (22-30°N, 105-118°E) in this study. Lag-lead correlation of ALI from preceding SON to JJA with (c) July and (d) August rainfall indices obtained from the UD rainfall dataset during 1901-2004. The dashed lines indicate significance at the 95% and 90% confidence level. (e) Sliding correlations between the winter (DJF) AL index and July South China rainfall index displayed at the central year of a 21-yr window. (f) As in (e) but for the August rainfall index. Horizontal dashed lines in (e, f) indicate the correlation coefficient is significant at the 90% confidence level. Note that all variables have been subjected to a nine-year high-pass filter and the DJF ENSO signal has been linearly removed in constructing this figure.

3.3 Associated atmospheric circulation anomalies in July and August

Changes consistent with rainfall anomalies are identified in July and August atmospheric circulation anomalies in association with the DJF ALI over East Asia and the western North Pacific. Figures 4a-c display the anomalies of geopotential height at 1000 hPa, winds at 850 hPa, and geopotential height at 500 hPa, in July, as obtained by regression on the DJF ALI during 1979-2004. At 1000 hPa, significant and positive geopotential height anomalies are seen over the subtropical western North Pacific and significant and negative anomalies are observed over the midlatitudes of the North Pacific and northwest of China (Fig. 4a). At 850 hPa, a significant anticyclonic circulation appears in the subtropical western North Pacific, together with easterly wind anomalies over the tropical western North Pacific and southwesterly wind anomalies over SC (Fig. 4b). These southwesterly wind anomalies, in association with the DJF ALI over SC, increase the northward transportation of moisture from the South China Sea and then concur with anomalous ascending motion there (not shown), finally resulting in positive rainfall anomalies over SC (Fig. 3a). The geopotential height anomalies over the North Pacific at 500 hPa are similar to those at 1000 hPa (Figs. 4a and 4c), with significant positive anomalies over the subtropical westernNorth Pacific and negative anomalies over the midlatitudes of the North Pacific. Note that significant negative geopotential height anomalies at 1000 hPa over Northwest China are replaced by positive anomalies at 500 hPa (Figs. 4a and 4c).

Figure 4 Anomalies of (a, d) 1000 hPa geopotential height (units: m), (b, e) 850 hPa horizontal winds (units: m s-1) and (c, f) 500 hPa geopotential height (units: m) in (a-c) July and (d-f) August regressed on the normalized DJF ALI during 1979-2004. Heavy (light) shading indicates geopotential anomalies that are significantly different from zero at the 95% (90%) confidence level. Note that all variables have been subjected to a nine-year high-pass filter and the DJF ENSO signal has been linearly removed in constructing this figure.

Correspondingly, Figs. 4d-f show atmospheric circulation anomalies in August, as obtained by regression on the DJF ALI. The anomalies in August display remarkable differences to those in July. There is an anomalous cyclone over the Philippine Sea in the lower troposphere (Fig. 4e), accompanied by negative geopotential height in the lower-to-middle troposphere (Figs. 4d and 4f), but with a slightly northeastward location in contrast with July. Significant northeasterly wind anomalies are observed over SC (Fig. 4e) associated with anomalous descending motion (not shown), which can explain the negative rainfall anomalies over SC in August (Fig. 3b).

At present, we find that the boreal winter AL index has a significant positive (negative) connection with the following July (August) SC rainfall. Hence, the preceding winter AL may be used as an effective predictor of July and August rainfall anomalies. The key system connecting the preceding winter AL and July (August) SC rainfall is the anticyclonic (cyclonic) circulation anomaly over the subtropical western Northern Pacific. However, the possible physical processes responsible for the formation of the winter AL-related anticyclonic (cyclonic) circulation anomalies in July (August) remain unclear. We speculate that winter AL-related boundary forcing (e.g., land/ocean) may play important roles in the formation of atmospheric circulation anomalies over the subtropical western North Pacific in the following summer.

4 Conclusion

This study investigates the influence of the winter AL on SC rainfall in the following summer. Results show that the wintertime AL is significantly positively (negatively) correlated with the SC rainfall in the following July (August). Specifically, SC rainfall anomalies tend to be positive (negative) when the preceding winter AL is stronger. The results of this study indicate that the preceding winter ALI can be used as a potential predictor of the following July and August rainfall anomalies over SC.

A significant anticyclonic (cyclonic) circulation anomaly can be observed over the subtropical western North Pacific in the following July (August), which can explain the pronounced positive (negative) correlation of SC rainfall in July (August) with the preceding winter AL. Southerly (northerly) wind anomalies to the west of the anomalous anticyclonic (cyclonic) circulation increase (decrease) the northward moisture transportation and thus contribute to the positive (negative) rainfall anomalies over SC in July (August).

Acknowledgments. We thank the two anonymous reviewers for their constructive suggestions and comments, which helped to improve the paper. This work was sponsored by the National Basic Research Program of China (Grant No. 2012CB955202).

References

Alexander, M. A., I. Blade, M. Newman, et al., 2002: The atmospheric bridge: The influence of ENSO teleconnections on air-sea interaction over the global oceans, J. Climate, 15, 2205-2231.

Allan, R., and T. Ansell, 2006: A new globally-complete historical gridded mean sea level pressure data set (HadSLP2), 1850-2004,J. Climate, 19, 5816-5842.

Compo, G. P., J. S. Whitaker, P. D. Sardeshmukh, et al., 2011: The twentieth century reanalysis project, Quart. J. Roy. Meteor. Soc., 137, 1-28.

Chen, S., W. Chen, and R. Wu, 2015: An interdecadal change in the relationship between boreal spring Arctic Oscillation and the East Asian summer monsoon around the early 1970s, J. Climate, 28, 1527-1542.

Duan, W., L. Y. Song, Y. Li, et al., 2013: Modulation of PDO on the predictability of interannual variability of early summer rainfall over south China, J. Geophys. Res., 118, 13008-13021.

Gong, D. Y., J. Yang, S. J. Kim, et al., 2011: Spring Arctic Oscillation-East Asian summer monsoon connection through circulation changes over the western North Pacific, Climate Dyn., 37, 2199-2216.

Huang, R., L. Zhou, and W. Chen, 2003: The progresses of recent studies on the variabilities of the East Asian monsoon and their causes, Adv. Atmos. Sci, 20(1), 55-69.

Legates, D. R., and C. J. Willmott, 1990: Mean seasonal and spatial variability in gauge corrected, global precipitation, Int. J. Climatol., 10, 111-127.

Liang, X. Z., and W. C. Wang, 1998: Association between China monsoon rainfall and tropospheric jets, Quart. J. Roy. Meteor. Soc., 124, 2597-2623.

Mao, J., J. C. Chan, and G. Wu, 2011: Interannual variations of early summer monsoon rainfall over South China under different PDO backgrounds, Int. J. Climatol., 31(6), 847-862.

Niebauer, H. J., N. A. Bond, L. P. Yakunin, et al., 1999: An update on the climatology and sea ice of the Bering Sea, in: Dynamics of the Bering Sea: A Summary of Physical, Chemical, and Biological Characteristics, and A Synopsis of Research on the Bering Sea, T. R. Loughlin (Eds.), University of Alaska Sea Grant, Fairbanks, 29-59.

Overland, J. E., J. M. Adams, and N. A. Bond, 1999: Decadal variability of the Aleutian Low and its relation to high-latitude circulation, J. Climate, 12, 1542-1548.

Park, Y. H., J. H. Yoon, Y.-H. Youn, et al., 2012: Recent warming in the western North Pacific in relation to rapid changes in the atmospheric circulation of the Siberian High and Aleutian Low systems, J. Climate, 25, 3476-3493.

Parkinson, C. L., 1990: The impact of the Siberian high and Aleutian low on the sea-ice cover of the Sea of Okhotsk, Ann. Glaciol., 14, 226-229.

Roeckner, E., J. M. Oberhuber, A. Bacher, et al., 1996: ENSO variability and atmospheric response in a global coupled atmosphere-ocean GCM, Climate Dyn., 12, 737-754.

Tao, S. Y., Q. Y. Zhang, and S. L. Zhang, 2001: An observational study on the behavior of the subtropical high over the west Pacific in summer, Acta Meteor. Sinica (in Chinese), 59, 747-758.

Trenberth, K. E., and J. W. Hurrell, 1994: Decadal atmospheric-ocean variations in the Pacific, Climate Dyn., 9, 303-319.

Tang, W., Z. Sun, and G. Tan, 2008: Effect of Indian Ocean Dipole on rainfall anomaly in south of China, J. Nanjing Inst. Meteor. (in Chinese), 31(6), 836-843.

Wang, B., J. Liu, J. Yang, et al., 2009: Distinct principal modes of early and late summer rainfall anomalies in East Asia, J. Climate, 22, 3864-3875.

Wang, B., R. Wu, and X. Fu, 2010: Pacific-East Asian teleconnection: How does ENSO affect East Asian climate? J. Climate, 13(9), 1517-1536.

Wu, R., Z. Z, Hu, and B. P. Kirtman, 2003: Evolution of ENSO-related rainfall anomalies in East Asia, J. Climate, 16, 3742-3758.

Wu, R., and B. P. Kirtman, 2007: Observed relationship of spring and summer East Asian rainfall with winter and spring Eurasian snow, J. Climate, 20(7), 1285-1304.

Wu, R., S. Yang, Z. Wen, et al., 2012: Interdecadal change in the relationship of southern China summer rainfall with tropical Indo-Pacific SST, Theor. Appl. Climatol., 108(1), 119-133.

Wu, Z. W, B. Wang, J. P. Li, et al., 2009: An empirical seasonal prediction model of the east Asian summer monsoon using ENSO and NAO, J. Geophys. Res., 114, doi:10.1029/2009JD011733.

Xiao, F., H. Jiang, B. Wu, et al., 2014: The leading correlation of the winter Aleutian Low with surface air temperature during the subsequent summer over the Arctic and its possible mechanism, Chin. Sci. Bull., 59(9), 904-912.

Xie, S. P., K. Hu, J. Hafner, et al., 2009: Indian Ocean capacitor effect on Indo-Western Pacific climate during the summer following El Ni?o, J. Climate, 22, 730-747.

Yang, S., and W. Lau, 2006: Interannual variability of the Asian monsoon, in: The Asian Monsoon, B. Wang (Ed.), Springer Praxis, 259-293.

Yu, J. Y., and S. T. Kim, 2011: Relationships between extratropical sea level pressure variations and the central Pacific and eastern Pacific types of ENSO, J. Climate, 24, 708-720.

Zhu, Y., and H. Wang, 2010: The Relationship between the Aleutian Low and the Australian summer monsoon at interannual time scales, Adv. Atmos. Sci., 27(1), 177-184.

Song, L.-Y., and W.-S. Duan, 2015: Interannual relationship between the winter Aleutian Low and rainfall in the following summer in South China,

Atmos. Oceanic Sci. Lett., 8, 271-276, 10.3878/ AOSL20150021.

3 March 2015; revised 3 May 2015; accepted 4 May 2015; published 16 September 2015

主站蜘蛛池模板: 国产性生交xxxxx免费| 特级精品毛片免费观看| 国产麻豆福利av在线播放| 最新精品国偷自产在线| 日本人妻丰满熟妇区| 亚洲专区一区二区在线观看| 国产视频 第一页| 亚洲国内精品自在自线官| 免费不卡视频| 日本一区二区三区精品视频| 欧美a级在线| 欧美黄色网站在线看| 久操线在视频在线观看| 在线免费a视频| 亚洲v日韩v欧美在线观看| 一本大道无码日韩精品影视| 九色国产在线| 久久精品这里只有国产中文精品| 国产精品久久久久久搜索| 欧美日韩亚洲国产| 国产一区二区人大臿蕉香蕉| 国产精品七七在线播放| 国产区91| 亚洲69视频| 白丝美女办公室高潮喷水视频 | 国产一级毛片在线| 精品久久久久无码| AV无码无在线观看免费| 在线观看免费黄色网址| 精品视频一区二区三区在线播 | 亚洲高清免费在线观看| 亚洲成人77777| 99一级毛片| 亚洲最大在线观看| 中文无码精品a∨在线观看| 免费一看一级毛片| 国产精品网拍在线| 激情综合网激情综合| 91午夜福利在线观看| 在线亚洲小视频| 99偷拍视频精品一区二区| 亚洲最大情网站在线观看| 国产手机在线小视频免费观看| 久久精品视频一| 97久久人人超碰国产精品| 亚洲人在线| 亚洲国产欧洲精品路线久久| 777午夜精品电影免费看| 最新精品久久精品| 成人av专区精品无码国产| 女人18一级毛片免费观看| 国产成人AV男人的天堂| 免费观看三级毛片| 欧美综合在线观看| 亚洲国产成人在线| 91久久国产成人免费观看| 99精品在线视频观看| 99国产在线视频| 免费国产好深啊好涨好硬视频| 熟妇人妻无乱码中文字幕真矢织江 | 欧美精品三级在线| 欧洲在线免费视频| www.91中文字幕| 67194成是人免费无码| 国产高清无码麻豆精品| 国产福利小视频高清在线观看| 2021最新国产精品网站| 欧洲精品视频在线观看| 日本一区二区三区精品AⅤ| 婷婷99视频精品全部在线观看 | 日韩人妻无码制服丝袜视频| 55夜色66夜色国产精品视频| 中文字幕永久在线观看| 啪啪免费视频一区二区| 国产福利小视频在线播放观看| 欧美日韩精品在线播放| 在线一级毛片| 特级aaaaaaaaa毛片免费视频| 91视频区| 伊人天堂网| 亚洲国产亚综合在线区| 狠狠v日韩v欧美v|