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

Revisiting Mid-Holocene Temperature over China Using PMIP3 Simulations

2015-11-24 10:20:43TIANZhiPingandJIANGDaBang

TIAN Zhi-Pingand JIANG Da-Bang,2

1Climate Change Research Center, Chinese Academy of Sciences, Beijing 100029, China

2Nansen-Zhu International Research Centre, Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100029, China

Revisiting Mid-Holocene Temperature over China Using PMIP3 Simulations

TIAN Zhi-Ping1and JIANG Da-Bang1,2

1Climate Change Research Center, Chinese Academy of Sciences, Beijing 100029, China

2Nansen-Zhu International Research Centre, Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100029, China

Using the simulations performed by 15 climate models under the latest protocol of the Paleoclimate Modeling Intercomparison Project (PMIP) Phase 3 (PMIP3), the authors revisited the annual and seasonal temperature changes over China during the mid-Holocene. Similar to the previous results produced by PMIP Phase 1 (PMIP1) and 2 (PMIP2) models, 14 (15) of the 15 PMIP3 models reproduced colder annual (boreal winter and spring) temperature in response to mid-Holocene insolation changes, with an average cooling of 0.33 K (1.31 K and 1.58 K) over China. The mid-Holocene boreal summer (autumn) temperature increased in all (13) of the 15 PMIP3 models, with an average warming of 1.02 K (0.61 K) at the national scale. Those changes simulated by the PMIP3 models were similar to those from the PMIP2 simulations but generally weaker than those from the PMIP1 models. A considerable mismatch still existed between the simulated cooling by the PMIP3 models and the reconstructed warming for annual and winter temperatures over China during the mid-Holocene, as was also the case between the previous PMIP1/2 simulations and proxy data.

Mid-Holocene, temperature over China, PMIP3 simulations

1 Introduction

The mid-Holocene, 6000 years before present, was a period when the climate and environment were remarkably different from today (Bartlein et al., 2011), mainly due to a change in incoming solar radiation at the top of the atmosphere (Berger, 1978). It is an ideal time period for investigating past climate change on the orbital scale, for which proxy data are abundant and relatively reliable. The mid-Holocene is also one of the benchmark periods within the framework of the Paleoclimate Modeling Intercomparison Project (PMIP), in which a hierarchy of climate models simulate the mid-Holocene climate either at the global (e.g., Braconnot et al., 2007; O'ishi and Abe-Ouchi, 2011) or regional scale over the monsoon regions of northern and southern Africa, North and South America, and northern Australia (e.g., Joussaume et al.,1999; Liu et al., 2004; Zhao and Harrison, 2012). Comparatively, however, less emphasis has been placed on the East Asian monsoon region (e.g., Wang et al., 2010; Zhou and Zhao, 2010, 2013), and thus there is an urgent need to further investigate the mid-Holocene East Asian climate.

Numerous efforts have been devoted to reconstructing the mid-Holocene climate over East Asia using a variety of proxy data, which suggest China experienced generally warmer and wetter than present conditions during that time. By contrast, some simulations focused on the mid-Holocene summer (June-August) climate over East Asia have been performed using atmospheric general circulation models (AGCMs) (e.g., Wang, 1999, 2002; Wang and Wang, 2013), regional climate models nested within AGCMs (e.g., Zheng et al., 2004), atmosphere-ocean general circulation models (AOGCMs) (e.g., Wei and Wang, 2004), and atmosphere-ocean-vegetation general circulation models (AOVGCMs) (Dallmeyer et al., 2010). In these previous studies, a warmer summer climate and stronger monsoon circulation over East Asia during the mid-Holocene have been reasonably reproduced, but there is a large degree of uncertainty among simulations. Additionally, little attention has been paid to the climate change for the annual mean and the other three seasons.

Recently, Jiang et al. (2012) examined annual and seasonal temperatures over China during the mid-Holocene using 36 climate models, including 16 AGCMs in PMIP Phase 1 (PMIP1) and 20 coupled models in Phase 2 (PMIP2). They found that 35 (36) of the 36 models reproduced colder-than-baseline annual (winter, December-February) temperature over China during that time—opposite to the warmer-than-present conditions derived from multiple proxy data sources. Similar results were obtained from the median of 10 PMIP Phase 3 (PMIP3) models (Zheng et al., 2013) and a state-of-the-art AOVGCM (Tian and Jiang, 2013). Furthermore, Jiang et al. (2012) and Liu et al. (2014) pointed out that the uncertainties in both PMIP simulations and proxy data may be responsible for the above model-data mismatch, in which quite different performances were found between PMIP1 AGCMs and PMIP2 coupled models, and different results could be interpreted from the reconstructed records with sparse spatial coverage. Considering that only the PMIP1/2 models and some PMIP3 models have been used in these earlier studies for the mid-Holocene climate over China, as well as the identified model-dependent uncertainties in the simulations, all available models of the latest PMIP3/Coupled Model Intercomparison Project Phase 5 (CMIP5) simulations need to be revisited and compared to multiple proxy data, in order to try uncover some clues about whether the inconsistency between the simulated cooling and reconstructed warming of the mid-Holocene annual and winter climates over China arises from the models, from the proxy data, or from both sides.

2 Data

2.1 Models and data

The present study was based on all 15 of the climate models under the latest PMIP3/CMIP5 protocol for mid-Holocene climate simulations, including eight AOGCMs and seven AOVGCMs, which have overall higher horizontal resolution than the former PMIP1/2 models. Basic information on the 15 models is provided in Table 1. The boundary conditions for the mid-Holocene experiment included changes in the Earth's orbital parameters (Berger, 1978), which played the most important role in the mid-Holocene climate, and atmospheric greenhouse gas concentrations. In PMIP3, the atmospheric concentrations of CO2, CH4, and N2O varied from the pre-industrial (baseline) levels of 284 ppm, 791 ppb, and 275 ppb, to 280 ppm, 650 ppb, and 270 ppb during the mid-Holocene, respectively. Sea surface temperatures (SSTs) were computed by oceanic general circulation models, and the vegetation was fixed at the present state in AOGCMs but was simulated with dynamic vegetation models in AOVGCMs for both the baseline and mid-Holocene experiments. More details about the models and experiments are provided by Taylor et al. (2012).

The temperature data used for model evaluation were taken from the CN05.2 monthly dataset for the period 1961-2000, with a half-degree horizontal resolution over China, provided by the National Climate Center of the China Meteorological Administration (Wu and Gao, 2013). Given that there are differences in the horizontal resolution among the models, all model and observational data were aggregated to a relatively mid-range resolution of 2.5° × 2° (longitude × latitude).

2.2 Evaluation of climate models

Whether or not the climate models can reasonably reproduce the present annual and seasonal temperature climatology over China, both in pattern and magnitude, will largely determine the confidence we have in their results for the mid-Holocene climate. In this sense, the spatial correlation coefficient (SCC), standard deviation, and centered root-mean-square difference (CRMSD) of each baseline experiment with respect to observation were calculated based on 179 grid points across mainland China. As shown in the Taylor diagram (Fig. 1), the SCCs ranged from 0.86 to 0.99 for annual and seasonal temperatures, all of which were statistically significant at the 99% confidence level, indicating reliable performance ofthe models in simulating the spatial pattern of the observed temperature climatology over China. Normalized standard deviations (CRMSDs) were 0.95-1.24 (0.26-0.52) for the annual mean, 0.90-1.16 (0.26-0.52) for winter, 0.97-1.44 (0.31-0.71) for spring (March-May), 0.87-1.21 (0.23-0.53) for summer, and 0.96-1.19 (0.21-0.41) for autumn (September-November) temperature. Generally, most models overestimated the spatial variability of annual and seasonal temperatures over China, since their standard deviations were larger than the observed values. They showed their best performance in autumn and poorest performance in spring. Note that most models performed better than PMIP1/2 models for modern temperature climatology, as indicated by their higher SCCs and lower CRMSDs (Jiang et al., 2012). In terms of the median of the 15 models, an SCC of 0.97-0.98, a normalized standard deviation of 1.02-1.14, and a normalized CRMSD of 0.21-0.31 were produced for annual and seasonal temperatures, showing a quantitatively higher reliability than most of the individual models (Fig. 1). Thus, the median of the 15 models is emphasized in the following analysis.

Table 1 Basic information about the 15 climate models within PMIP3 for the mid-Holocene experiments used in this study.

3 Mid-Holocene temperature over China

3.1 Annual temperature changes

With respect to the baseline period, 14 out of 15 PMIP3 models reproduced colder annual temperatures over China during the mid-Holocene (Fig. 2). Quantitatively, the regionally averaged change over China ranged from -0.97 K (FGOALS-g2; the model full name is provided in Table 1, hereafter the same) to 0.32 K (HadGEM2-ES), with two models (CNRM-CM5 and HadGEM2-CC) showing a very weak cooling of less than 0.1 K and only one model (HadGEM2-ES) indicating a warming. Though it differed somewhat from model to model, the simulated magnitude of annual temperature change was similar to that from the 20 PMIP2 coupled models but generally weaker than that from the 16 PMIP1 AGCMs (Jiang et al., 2012). Considering the median of the 15 models, annual temperature averaged over China was reduced by 0.33 K during the mid-Holocene relative to the baseline period, with a standard deviation of 0.33 K across individual models. The geographical distribution ofthe mid-Holocene annual temperature change in the 15-model median showed a generally weak cooling, with a high consistency among the models, over most of China, but a weak warming over northern Northeast China, most of Xinjiang, and part of North and central China (Fig. 3a). Annual temperature change was less than 0.6 K over most of China, except over the eastern Tibetan Plateau and South China where the cooling was 0.6-1.2 K. The spatial pattern of annual temperature change over China simulated by the 15 PMIP3 models was generally similar to the 10 PMIP3 model result in Zheng et al. (2013), but significant cooling was only shown in part of South China and warmer temperature was simulated over most of the Tibetan Plateau, excluding its central part, in the 10-model median (Zheng et al., 2013).

In response to mid-Holocene orbital forcing, SSTs changed significantly with respect to the baseline period according to the 12 models with SST data available (Fig. 4). Annual SSTs generally increased north of 30°N but decreased south of it in the Northern Hemisphere, with the latter change being consistent in most models and thus playing a major role (Fig. 4a). As a result, the consistently colder SSTs in the oceans adjacent to the East Asian continent could lead to larger losses of surface heat over East Asia during warm months and smaller gains of surface heat during cold months, inducing an overall decrease insurface air temperature over China (Fig. 3a).

3.2 Seasonal temperature change

As mentioned in Jiang et al. (2012) and Tian and Jiang (2013), the mid-Holocene change in seasonal temperature generally followed the insolation change over China, which was reduced in winter, spring, and autumn but increased in summer relative to the present day. In response to those changes, all of the 15 PMIP3 models reproduced colder winter and spring temperatures over China during the mid-Holocene, with the national-averaged change ranging from -2.13 K to -0.76 K and from -4.18 K to -0.94 K, respectively (Fig. 2). The median of the 15 models showed a winter (spring) cooling of 1.31 K (1.58 K) averaged over China, with a standard deviation of 0.40 K (0.80 K) across the models. The cooling magnitude (standard deviation) of winter and spring temperatures in the PMIP3 models was comparable to that in the PMIP2 models, but generally smaller (larger) than the PMIP1 AGCM results (Jiang et al., 2012). In summer, temperature increased by 0.50-2.52 K in all of the 15 models averaged over China during the mid-Holocene, with a warmingof 1.02 K in the median of the 15 models. The warming magnitude was quantitatively comparable to that obtained from PMIP1/2 models (Jiang et al., 2012). Autumn temperature over China was increased by 0.04-5.08 K in 13 models but decreased by 0.10 K in GISS-E2-R and by 0.37 K in FGOALS-g2 during the mid-Holocene, with a median warming of 0.61 K and a standard deviation of 1.29 K across the models. Thus, there was a large degree of scatter in the simulated autumn temperature among the PMIP3 models. Note that the overall autumn warming over China did not follow the reduced autumn insolation during the mid-Holocene, which was also the case in the PMIP2 models but opposite to the cooling reproduced by the PMIP1 models (Jiang et al., 2012). This may be because of the different baseline period used in the PMIP1, PMIP2, and PMIP3 models, which was the present day (ca. 1985) in PMIP1 but the pre-industrial period (ca. 1750) in PMIP2/3. It is also worth mentioning that there were large seasonal climate changes during the mid-Holocene in HadGEM2-ES, showing strongest cooling in winter and spring but strongest warming in summer and autumn among the 15 PMIP3 models (Fig. 2). This was because HadGEM2-ES reproduced relatively lower (higher) baseline temperature but higher (lower) mid-Holocene temperature in summer and autumn (winter and spring) compared to the other 14 models.

As shown from the geographical distribution based on the median of the 15 models, the mid-Holocene winter and spring temperatures were reduced by 0.6-2.4 K across China, both with a high consistency of cooling among models (Figs. 3b and 3c). Summer temperature was increased by 0-2.4 K over most of China, except for a cooling over part of the southern Tibetan Plateau and South China, where model consistency was less than 80% (Fig. 3d). In autumn, temperature rose over most of China, with a consistent warming of 0.3-1.8 K over the region west of 110°E but a change smaller than 0.6 K east of 110°E over China (Fig. 3e). As a whole, the large-scale spatial pattern of seasonal temperature changes over China during the mid-Holocene reproduced by the 15 PMIP3 models resembled that simulated by the 36 PMIP1/2 models (Jiang et al., 2012), except in autumn when the simulated magnitude was overall larger for the former than for the latter. The mid-Holocene winter and summer temperature changes over China in the 15-model median were also similar to the 10-model median results in Zheng et al. (2013), except for a weak warming over part of the Tibetan Plateau in winter and a warming across the whole of China in summer for the latter.

During the mid-Holocene winter and spring, SSTs increased slightly at northern high latitudes but decreased with a high consistency among the models over the oceans adjacent to the East Asian continent (Figs. 4b and 4c), leading to an overall cooling over the whole of China in those seasons (Figs. 3b and 3c). Similar to the annual mean, the mid-Holocene summer SSTs increased by 0-1.0 K north of 30°N but decreased by 0-0.6 K south of it, both with a high consistency among models (Fig. 4d), which favored a surface warming over most of China but a slight cooling south of 30°N over China in summer (Fig. 3d). During the mid-Holocene autumn, the 12-model median reproduced generally warmer SSTs in the Northern Hemisphere, contributing to the increased autumn temperature over most of China (Figs. 3e and 4e).

3.3 Model-data comparison

As derived from the 64 records of pollen, lake cores, paleosol, ice cores, peat, sediment, stalagmites, and fossil fruits detailed in Jiang et al. (2012), the mid-Holocene annual temperature over China was reconstructed to be warmer than today. In addition, a number of datasets of pollen samples from paleovegetation have all indicated warmer winter temperatures over China during the mid-Holocene (Yu et al., 2000; Ni et al., 2010), which was even stronger than the annual warming (Xu et al., 1988; Kong et al., 1990; Shi et al., 1993; Guiot et al., 2008; Jiang et al., 2010). Contrary to the reconstruction, colder-than-baseline annual (winter) climate was reproduced by 14 (15) of the 15 PMIP3 models, with a cooling of 0.33 K (1.31 K) for the median of the 15 models, averaged over China during the mid-Holocene (Fig. 2). Concerning the spatial pattern changes, the simulated mid-Holocene annual (winter) temperature was consistently colder than the baseline period over most (all) of China in the PMIP3 models (Figs. 3a and 3b). However, a slight warming was reproduced by the median of the 15 PMIP3 models over northern Northeast China, most of Xinjiang, and part of North and central China, with a high consistency among individual models, except for the former region (Fig. 3a), and thus was qualitatively in accordance with proxy data, though the spatial coverage was very sparse over those areas.

Based on the above, a considerable mismatch still exists between proxy data and PMIP3 simulations, as was also the case between the reconstruction and PMIP1/2 simulations (Jiang et al., 2012). On the other hand, as discussed in Jiang et al. (2012), the colder-than-baseline annual and winter climates reproduced by PMIP3 models were mainly attributed to the negative radiative forcing over China during the mid-Holocene, which were consistently robust between the models, albeit with some uncertainties still apparent. Additionally, uncertainty in the proxy data may be another source for the model-data inconsistency. In spite of the robust warmer annual and winter climates over China as derived from multiple proxy data, the mid-Holocene annual temperature was indicated by Guiot et al. (2008) and Bartlein et al. (2011) to be colder in parts of China. Furthermore, Liu et al. (2014) pointed out that the above reconstructions were mainly based on vegetation records, and thus a seasonal bias and precipitation effects cannot be ruled out to reconcile the contradiction between the reconstructed cooling and the simulated warming. Collectively, further study should focus on more reconstruction work as well as climate model improvement, to narrow the spread between proxy data and simulations.

4 Conclusion

The mid-Holocene annual and seasonal temperature changes over China were revisited using the 15 climate models under the latest PMIP3/CMIP5 protocol. Results showed that, with respect to the baseline period, the mid-Holocene changes in annual and seasonal temperatures over China reproduced by the PMIP3 models were similar to the PMIP1/2 simulations in Jiang et al. (2012). Following the mid-Holocene annual and seasonal insolation changes over China, 14 (15) of the 15 PMIP3 models reproduced colder-than-baseline annual (winter and spring) temperature, whereas all (13) of the 15 models simulated warmer summer (autumn) climate over China during that period. Averaged over China for the median of the 15 models, the mid-Holocene changes in annual, winter, spring, summer, and autumn temperature were -0.33 K, -1.31 K, -1.58 K, 1.02 K, and 0.61 K, respectively —similar to those from the 20 PMIP2 coupled models but generally weaker than those from the 16 PMIP1 AGCMs. A considerable mismatch still exists between the reconstructed warming and the simulated cooling in the PMIP3 models for annual and winter temperatures over China during the mid-Holocene, except over northern Northeast China, most of Xinjiang, and part of North and central China, where a slight warming was reproduced by most of the PMIP3 models.

Acknowledgments. We sincerely thank the anonymous reviewers for their insightful comments, which improved the paper. We also thank the climate modeling groups (listed in Table 1) for producing and making available their model output. This research was supported by the National Natural Science Foundation of China (Grant No. 41222034).

References

Bartlein, P., S. Harrison, S. Brewer, et al., 2011: Pollen-based continental climate reconstructions at 6 and 21 ka: A global synthesis, Climate Dyn., 37, 775-802.

Berger, A. L., 1978: Long-term variations of daily insolation and quaternary climatic changes, J. Atmos. Sci., 35, 2362-2367.

Braconnot, P., B. Otto-Bliesner, S. Harrison, et al., 2007: Results of PMIP2 coupled simulations of the mid-Holocene and last glacial maximum—Part 1: Experiments and large-scale features, Climate Past, 3, 261-277.

Dallmeyer, A., M. Claussen, and J. Otto, 2010: Contribution of oceanic and vegetation feedbacks to Holocene climate change in monsoonal Asia, Climate Past, 6, 195-218.

Guiot, J., H. B. Wu, W. Y. Jiang, et al., 2008: East Asian Monsoon and paleoclimatic data analysis: A vegetation point of view, Climate Past, 4, 137-145.

Jiang, D., X. Lang, Z. Tian, et al., 2012: Considerable model-data mismatch in temperature over China during the mid-Holocene: Results of PMIP simulations, J. Climate, 25, 4135-4153.

Jiang, W., J. Guiot, G. Chu, et al., 2010: An improved methodology of the modern analogues technique for palaeoclimate reconstruction in arid and semi-arid regions, Boreas, 39, 145-153.

Joussaume, S., K. E. Taylor, P. Braconnot, et al., 1999: Monsoon changes for 6000 years ago: Results of 18 simulations from the Paleoclimate Modeling Intercomparison Project (PMIP), Geophys. Res. Lett., 26, 859-862.

Kong, Z. C., N. Q. Du, F. S. Shan, et al., 1990: Vegetational and climatic changes in the last 11,000 years in Qinghai Lake— Numerical analysis based on palynology in core QH85-14C, Mar. Geol. Quat. Geol. (in Chinese), 10, 79-90.

Liu, Z., J. Zhu, Y. Rosenthal, et al., 2014: The Holocene temperature conundrum, Proc. Natl. Acad. Sci. USA, 111, E3501-E3505, doi:10.1073/pnas.1407229111.

Liu, Z., S. P. Harrison, J. Kutzbach, et al., 2004: Global monsoons in the mid-Holocene and oceanic feedback, Climate Dyn., 22, 157-182.

Ni, J., G. Yu, S. P. Harrison, et al., 2010: Palaeovegetation in China during the late Quaternary: Biome reconstructions based on a global scheme of plant functional types, Palaeogeogr. Palaeoclimatol. Palaeoecol., 289, 44-61.

O'ishi, R., and A. Abe-Ouchi, 2011: Polar amplification in the mid-Holocene derived from dynamical vegetation change with a GCM, Geophys. Res. Lett., 38, L14702, doi: 10.1029/ 2011GL048001.

Shi, Y., Z. Kong, S. Wang, et al., 1993: Mid-Holocene climates and environments in China, Glob. Planet. Change, 7, 219-233.

Taylor, K. E., 2001: Summarizing multiple aspects of model performance in a single diagram, J. Geophys. Res., 106, 7183-7192.

Taylor, K. E., R. J. Stouffer, and G. A. Meehl, 2012: An overview of CMIP5 and the experiment design, Bull. Amer. Meteor. Soc., 93, 485-498.

Tian, Z., and D. Jiang, 2013: Mid-Holocene ocean and vegetation feedbacks over East Asia, Climate Past, 9, 2153-2171.

Wang, H. J., 1999: Role of vegetation and soil in the Holocene megathermal climate over China, J. Geophys. Res., 104, 9361-9367.

Wang, H. J., 2002: The mid-Holocene climate simulated by a grid-point AGCM coupled with a biome model, Adv. Atmos. Sci., 19, 205-218.

Wang, T., and H. J. Wang, 2013: Mid-Holocene Asian summer climate and its responses to cold ocean surface simulated in the PMIP2 OAGCMs experiments, J. Geophys. Res., 118, 4117-4128.

Wang, T., H. J. Wang, and D. Jiang, 2010: Mid-Holocene East Asian summer climate as simulated by the PMIP2 models, Palaeogeogr. Palaeoclimatol. Palaeoecol., 288, 93-102.

Wei, J. F., and H. J. Wang, 2004: A possible role of solar radiation and ocean in the mid-Holocene East Asian monsoon climate, Adv. Atmos. Sci., 21, 1-12.

Wu, J., and X. J. Gao, 2013: A gridded daily observation dataset over China region and comparison with the other datasets, Chin. J. Geophys., 56, 1102-1111, doi:10.6038/cjg20130406.

Xu, Q. H., S. Y. Chen, Z. C. Kong, et al., 1988: Preliminary discussion of vegetation succession and climatic change since the Holocene in the Baiyangdian Lake district, Acta Phytoecol. Geobotan. Sinica (in Chinese), 12, 143-151.

Yu, G., X. Chen, J. Ni, et al., 2000: Palaeovegetation of China: A pollen data-based synthesis for the mid-Holocene and last glacial maximum, J. Biogeogr., 27, 635-664.

Zhao, Y., and S. P. Harrison, 2012: Mid-Holocene monsoons: A multi-model analysis of the inter-hemispheric differences in the responses to orbital forcing and ocean feedbacks, Climate Dyn., 39, 1457-1487.

Zheng, W., B. Wu, J. He, et al., 2013: The East Asian Summer Monsoon at mid-Holocene: Results from PMIP3 simulations, Climate Past, 9, 453-466.

Zheng, Y. Q., G. Yu, S. M. Wang, et al., 2004: Simulation of paleoclimate over East Asia at 6 ka BP and 21 ka BP by a regional climate model, Climate Dyn., 23, 513-529.

Zhou, B., and P. Zhao, 2010: Modeling variations of summer upper tropospheric temperature and associated climate over the Asian Pacific region during the mid-Holocene, J. Geophys. Res., 115, D20109, doi:10.1029/2010JD014029.

Zhou, B., and P. Zhao, 2013: Simulating changes of spring Asian-Pacific oscillation and associated atmospheric circulation in the mid-Holocene, Int. J. Climatol., 33, 529-538.

Tian, Z.-P., and D.-B. Jiang, 2015: Revisiting mid-Holocene temperature over China using PMIP3 simulations, Atmos. Oceanic Sci. Lett., 8, 358-364,

10.3878/AOSL20150040.

28 April 2015; revised 17 June 2015; accepted 26 June 2015; published 16 November 2015

TIAN Zhi-Ping, tianzhiping@mail.iap.ac.cn

主站蜘蛛池模板: 欧美亚洲国产精品第一页| 日韩精品久久久久久久电影蜜臀| 首页亚洲国产丝袜长腿综合| 久热99这里只有精品视频6| 毛片网站观看| 精品福利国产| 日韩黄色大片免费看| 欧美影院久久| 国产乱人伦AV在线A| 亚洲第一页在线观看| 亚洲无码电影| 亚洲一区二区三区香蕉| 尤物成AV人片在线观看| 成人噜噜噜视频在线观看| 精品无码日韩国产不卡av| 香蕉国产精品视频| 91九色视频网| 色婷婷丁香| 亚洲男女在线| 制服丝袜一区| 日韩av在线直播| 日本免费一级视频| 999国内精品久久免费视频| 欧美午夜在线播放| 欧美日韩精品一区二区在线线| 欧美一区日韩一区中文字幕页| 国产在线精品香蕉麻豆| 中文纯内无码H| 亚洲中文字幕无码爆乳| 午夜视频免费试看| 成人在线观看不卡| 青青操国产视频| 精品少妇三级亚洲| a毛片在线播放| 午夜精品一区二区蜜桃| 日韩精品久久久久久久电影蜜臀| 久久特级毛片| 欧美一级片在线| 久久五月天综合| 狂欢视频在线观看不卡| 一本大道东京热无码av| 国产一级视频久久| 欧美a网站| 国产一区二区在线视频观看| 高清精品美女在线播放| 人妻中文字幕无码久久一区| 91福利在线看| 亚洲美女一区二区三区| 2021国产精品自拍| 国产白浆在线| 久久综合久久鬼| 最新精品国偷自产在线| 欧美啪啪精品| 国产成人一区免费观看| 日韩精品亚洲精品第一页| 九九热精品在线视频| 高清亚洲欧美在线看| 亚洲动漫h| 真实国产乱子伦视频| 欧美国产成人在线| 啪啪国产视频| 免费一级毛片完整版在线看| 欧美日韩中文国产| 最新国语自产精品视频在| 亚洲人成在线精品| 国产欧美另类| 国产精品私拍在线爆乳| 在线看片国产| yjizz国产在线视频网| 99久久精品免费观看国产| 欧美精品一二三区| 日本精品影院| 女人18毛片久久| 国产毛片片精品天天看视频| 一本大道香蕉高清久久| 男人天堂亚洲天堂| 一级片免费网站| 久久久久青草大香线综合精品 | 一本久道久综合久久鬼色| 欧美a在线看| 国产女同自拍视频| 日韩在线欧美在线|