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

    Datasets for the CMIP6 Scenario Model Intercomparison Project(ScenarioMIP) Simulations with the Coupled Model CAS FGOALS-f3-L

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

    Shuwen ZHAO, Yongqiang YU*,3, Pengfei LIN, Hailong LIU, Bian HE, Qing BAO,Yuyang GUO, Lijuan HUA, Kangjun CHEN, and Xiaowei WANG,5

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

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

    3Center for ocean Mega-Science, Chinese Academy of Sciences, Qingdao 266071, China

    4Chinese Academy of Meteorological Sciences, Beijing 100081, China

    5Yunnan University, Kunming 650091, China

    ABSTRACT

    Key words:CMIP6,ScenarioMIP,FGOALS-f3-L,coupled model,data description

    1.Introduction

    The state-of-the-art earth/climate system model is an important tool to investigate multi-sphere interaction of the climate system as well as reproduce the past climate change and project future climate change. A modeling group from the State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics (LASG),Institute of Atmospheric Physics (IAP), has been developing the Flexible Global Ocean–Atmosphere–Land System model (FGOALS) since the early 2000s (Yu et al., 2002,2004). Its second spectral version, FGOALS-s2 (Bao et al.,2013; Zhou et al., 2014b), was released in 2012 and contributed to phase 5 of the Coupled Model Intercomparison Project (CMIP5) (Taylor et al., 2012). This model shows good simulation performance in global and regional climate characteristics, such as sea surface temperature (SST),precipitation, and atmospheric circulation. Specifically,FGOALS-s2 reproduces the annual and semi-annual SST cycles in the equatorial Pacific Ocean, the main characteristics of the Asian summer monsoon, the extratropical teleconnection patterns of large-scale climate, the irregularity of ENSO periods, etc. (Bao et al., 2013). The third version is FGOALS-f3-L, and has been employed to carry out many numerical experiments following the CMIP6 design (Eyring et al., 2016).

    The long-term trend of deep ocean in the Pi-Control experiment of FGOALS-f3-L is smaller than that of the previous version (Guo et al, 1996), so responses of the coupled model to the external forcings in the Scenario Model Intercomparison Project (ScenarioMIP) experiment are less affected as before by long-term climate drift. This suggests that FGOALS-f3-L is more reliable in its future climate projections. In addition, model biases in the tropical region are also significantly reduced, which further narrows the uncertainty of the projected climate change.

    As one of primary activities in CMIP6, ScenarioMIP will improve understanding of climate changes under different future scenarios and their impacts on society. According to the model participation information provided by the CMIP6 website (https://rawgit.com/WCRP-CMIP/CMIP6_CVs/master/src/CMIP6_source_id.html), more than 50 models around the world have participated in the model comparison program so far. The ScenarioMIP experiments are driven by new emissions and land-use scenarios (Riahi et al., 2017), which are produced with integrated assessment models (IAMs) considering the Shared Socioeconomic Pathways (SSPs) and future pathways of societal development.The ScenarioMIP aims at: (1) facilitating integrated research in order to better understand the impact of different scenarios on the physical processes of the climate system and the impact of climate change on society; (2) providing a data basis for ScenarioMIP specific scientific issues and other CMIP6 projects, aimed at the climate impact of particular forcing relative to scenario projections; and (3) providing a basis for the development of a new method of quantitative assessment and prediction of uncertainty based on multi-model ensembles, considering model performance,model dependence and observational uncertainty (O’Neill et al., 2016; Zhang et al., 2019).

    The ScenarioMIP simulations of FGOALS-f3-L were completed in October 2019 and the data have been released on the Earth System Grid Federation (ESGF) data server(https://esgf-node.llnl.gov/projects/cmip6/). To provide a description of the ScenarioMIP model outputs, we present some fundamental analysis on some key output variables,such as temperature, precipitation etc. Descriptions of the model and experimental design are documented in section 2. The projections of long-term climate change are described in section 3. And finally, section 4 provides data records and usage notes.

    2.Model and experiments

    2.1.Model description

    The climate model CAS FGOALS-f3-L includes four individual components: an atmospheric component, version 2.2 of the Finite-volume Atmospheric Model (FAMIL)(Zhou et al., 2015; Bao et al., 2019; Li et al., 2019), which has a horizontal resolution of 1° (lat) by 1.25° (lon); an oceanic component, version 3 of the LASG/ IAP Climate Ocean Model (LICOM3) (Lin et al., 2020), in which the horizontal resolution is also proximately equal to 1°; a land component, version 4.0 of the Community Land Model (CLM4)(Oleson et al., 2010), which has a horizontal resolution of 0.9° (lat) by 1.25° (lon); and a sea-ice component, version 4 of the Los Alamos Sea Ice Model (CICE4) (Hunke and Lipscomb, 2010), in which the horizontal resolution is also proximately equal to 1°. A flux coupler module, version 7 of NCAR’s coupler (Craig, 2014), is used to connect these four components.

    2.2.Experimental design

    2.2.1.Introduction to the scenarios

    CMIP, now in its sixth phase (CMIP6), organized by the World Climate Research Programme’s Working Group on Coupled Modelling, has become one of the cornerstones of climate science. ScenarioMIP, one of the major activities of CMIP6, can solve the problem raised in CMIP6 about assessing future climate changes induced by human activity.

    For the first time, future scenario experiments are to be run as a MIP. The new scenarios in ScenarioMIP are based on the SSPs and forcing levels. The SSPs describe different scenarios of social development in the future without the impact of climate change or climate policies. SSP1, SSP2,SSP3, SSP4 and SSP5 represent pathways of sustainable development, middle of the road, regional rivalry, inequality, and fossil-fueled development, respectively (O’Neill et al., 2016; Zhang et al., 2019).

    Considering future emissions, greenhouse gas (GHG)concentrations and particular land use, the design of ScenarioMIP involves eight pathways, including two tiers (tier-1 and tier-2) of priority. The scenarios in tier-1 are important for climate science, IAM, and IAV (impacts, adaptation and vulnerability) research, and there are four experiments in tier-1 under different forcing scenarios, named SSP1-2.6,SSP2-4.5, SSP3-7.0 and SSP5-8.5. There are other scenarios of interest, including additional 21st century scenarios,initial condition ensemble members for the SSP3-7.0 scenario, and long-term extension scenarios. Due to the limited computer resource, only tier-1 experiments are conducted with CASFGOALS-f3-L, and thus we only describe some details of tier-1 as shown in Table 1. SSP1-2.6, SSP2-4.5 and SSP5-8.5 are updates of the previous RCP2.6, RCP4.5 and RCP8.5 in CMIP5, respectively. SSP3-7.0 is a new scenario, which represents the medium to high end of the rangeof future forcing pathways (O’Neill et al., 2016; Zhang et al., 2019). In general, ScenarioMIP’s adoption of SSP pathways can provide more reliable information about climate change than before to policymakers.

    Table 1. ScenarioMIP experiment designs.

    2.2.2.Experiment forcings

    Anthropogenic forcings in ScenarioMIP experiments consist of land use, emissions of long-lived GHGs, concentrations of long-lived GHGs, and emissions of air pollutants.In addition, there are some requirements for these SSP simulations as follows: The model version participating in ScenarioMIP must be identical to that participating in the CMIP6 historical run, thus ensuring continuity of climate simulation.All ScenarioMIP experiments must use the ScenarioMIPprovided concentrations for all long-lived GHGs, such as CO, CH, NO and CFCs. The concentrations of CO, CHand NO in our model are shown in Fig.1, in which it can be seen that, for concentrations of CHand NO, the concentrations in SSP3-7.0 are higher than those in SSP5-8.5.However, this feature is not reflected in the time series of anthropogenic radiative forcing (O’Neill et al., 2016). COstill dominates the GHGs in the future of the 21st century.In addition, ScenarioMIP provides other anthropogenic forcings (i.e., short-lived GHGs, aerosol, and ozone), up-todate solar radiation time series (Matthes et al., 2017), and volcanic forcing prescribed in terms of historical volcanic eruptions. Moreover, as a climate system model, the coupled model CASFGOALS-f3-L is driven by concentrations of GHGs in ScenarioMIP experiments. All forcing data from CMIP6 are bilinearly interpolated into the grid points of the atmospheric component model.

    3.Projections of long-term climate change

    A portion of the ScenarioMIP datasets of the CAS FGOALS-f3-L have been submitted online. In this section,the responses of the climate system to the radiative forcings from ScenarioMIP experiments are presented.

    3.1.Temporal evolution in SSP experiments

    Figure 2 shows the annual global mean net radiation at the top of the atmosphere (TOA), surface air temperature at 2 m, sea-ice extent in the Northern Hemisphere in September, the upper 700 m and 700–2000 m ocean heat content anomalies, sea level rise, and precipitation in 2015–2100,for SSP1-2.6, SSP2-4.5, SSP3-7.0 and SSP5-8.5. The upper 700 m and 700–2000 m ocean heat content anomalies are calculated relative to the 20-year mean from 1995 to 2014 in historical runs of our model. The sea level rise is calculated based on the density change relative to the 1995–2014 climatological mean state that resulted from heat expansion (Landerer et al., 2007).

    The net radiation at the TOA (Fig.2a) in the four different scenarios are positive, which means the energy of the earth system is imbalanced and increased. With gradual stabilization of GHG concentrations for SSP1-2.6 (black curves in Fig.1), TOA net radiation presents a definite decreasing trend for SSP1-2.6 and the value is getting closer to 0.5 W maround 2100, revealing a quasi-stabilization characteristic.This shows that the net radiation at the TOA will gradually become balanced and a future equilibrium state will be achieved after a few centuries. Also, there is a slight upward trend of TOA net radiation during 2015–60 and a downward trend after about 2060, for SSP2-4.5. However, TOA net radiation continues to increase for SSP3-7.0 and SSP5-8.5 till 2100. In fact, these changes in the net TOA radiation budget basically correspond to the temporal evolution of GHG concentrations in the 21st century in Fig.1.However, it should be noted that the experiments are still beyond stabilization at the end of the 21st century, except for the SSP1-26 experiment, because the climate system might take a very long time (perhaps centuries) to reach a new steady state as a response to sufficiently strong external radiation forcings. The variation trend of TOA net radiation in ScenarioMIP is similar to that in CMIP5 in climate projection experiments; it is a little more than about 0.2 W mabove the aggregate mean of CMIP5 (Collins et al., 2013).

    Fig.1. Time series of (a) CO2 (units: ppm), (b) CH4 (units: ppm), and (c) N2O (units: 10?3 ppm) concentrations from 2015 to 2100 in ScenarioMIP experiments.

    Fig.2. The annual global mean (a) TOA net radiation (units: W m?2, 5-year running mean), (b) surface air temperature at 2 m (SAT, units: °C, 5-year running mean), (c) Arctic sea-ice extent (units: 106 km2, 5-year running mean, dashed line represents near ice-free conditions), (d) upper 700 m ocean heat content anomalies relative to the 20-year mean from 1995 to 2014 in historical runs (0–700 m OHCA, units: 1023 J, annual mean), (e) 700–2000 m ocean heat content anomalies relative to the 20-year mean from 1995 to 2014 in historical runs (700–2000 m OHCA,units: 1023 J, annual mean), (f) sea level rise due to heat expansion relative to the 20-year mean from 1995 to 2014 in historical runs (sea level rise, units: cm, annual mean), and (g) precipitation (units: mm d?1, 5-year running mean)from 2015 to 2100, for SSP1-2.6, SSP2-4.5, SSP3-7.0 and SSP5-8.5.

    A significant feature of climate change is the changes of surface air temperature. It can be seen from the time series of annual global mean surface air temperature (Fig.2b) that all four scenarios show persistent warming throughout the whole of the 21st century. Temperature increases are almost the same for all the scenarios during the first 10 years and then differ from each other after 2026.The warming rate is highest in SSP5-8.5, and remarkably lower in SSP1-2.6, in which the global mean surface air temperature reaches a quasi-steady state. Specifically, for SSP1-2.6, SSP2-4.5, SSP3-7.0 and SSP5-8.5, projected global mean surface air temperatures between 2080 and 2100 relative to the climatological mean between 1850 and 1900 from historical runs increase by 1.97°C, 2.84°C,3.91°C and 4.66°C, respectively. Compared with the CMIP5 multi-model ensemble mean temperature changes,the coupled model simulates higher temperature changes(Collins et al., 2013), which should be attributable to the larger climate sensitivity of FGOALS-f3-L. The climate sensitivity of our model is about 3.0°C and slightly larger than the result of the CMIP5 multi-model ensemble mean (Collins et al., 2013; He et al., 2020).

    In addition to temperature, the change in sea-ice extent,particularly in the Arctic, is a crucial feature of climate change. Also, September is the end of the summer melting season. The annual sea-ice extents in the Arctic in September (Fig.2c) in SSP2-4.5, SSP3-7.0 and SSP5-8.5 exhibit downward trends and appear to have steepened with time.An ice-free summer condition is defined by an ice extent of less than 10km(Stroeve et al., 2012). The Arctic will be nearly ice-free in September for SSP2-4.5 after 2049, for SSP3-7.0 after 2044, and for SSP5-8.5 after 2041, which results from the strong warming. The September sea-ice extent in the Arctic shows a rapid reduction during 2015–55 and fluctuates around 10kmfor SSP1-2.6 after 2055. The September sea-ice extent in the Arctic of the CMIP5 multimodel ensemble mean decreases faster than the simulated sea-ice extent in FGOALS-f3-L. Meanwhile, in the CMIP5 RCP2.6, the September sea-ice extent in the Arctic drops to below 10kmaround 2060 and continues to decline to 0 km. The main reason for this difference is the different pathway settings.

    Not only is the increased TOA net radiation used to heat the atmosphere, but some of it is stored in the oceans.Insight into the ocean response to climate change can be gained by examining the changes of ocean heat content. The annual global mean upper 700 m (Fig.2d) and 700–2000 m(Fig.2e) ocean heat content anomalies increase for all four scenarios in 2015–2100, implying that the ocean will become warmer in response to positive forcings. When the SST increases as a result of positive forcing, the interior water masses respond to the integrated signal at the surface and then propagate down to greater depth. The differences in the upper 700 m ocean heat content anomalies between the four scenarios are not significant before 2040, but become more significant after 2040. The differences in the 700–2000 m ocean heat content anomalies between the four scenarios are not significant before 2050, but become more significant after 2050. Compared with the time series of surface air temperature, the results show that the deep ocean changes slowly with climate change. Based on the related calculation and analysis, the heat absorbed by the ocean accounts for about 76.4%, 78.9%, 83.6% and 84.6% of the added heat of the earth system from 2015 to 2100 for SSP1-2.6, SSP2-4.5, SSP3-7.0 and SSP5-8.5, respectively, which means most of the increased earth energy is stored in the oceans and the proportion of the increased energy stored in the ocean increases as the positive forcing increases.

    The observed sea level rise is mainly caused by thermal expansion of seawater and melting of land glaciers. Thermal expansion of seawater contributes about one third to the observed rise, while melting of land glaciers contributes the other two thirds (Cazenave and Remy, 2011). In a climate system model without a glacier component, only the contribution of thermal expansion to sea level, i.e., thermosteric sea level, can be estimated. The global mean thermosteric sea level rise between 2015 and 2100 is shown in Fig.2f. The features of the time series of thermosteric sea level for all four scenarios are similar to the upper 700 m ocean heat content,which illustrates that the upper 700 m ocean warming is crucial to thermosteric sea level rise. For SSP1-2.6, SSP2-4.5,SSP3-7.0 and SSP5-8.5, the average thermosteric sea levels rise is 1.29, 1.88, 2.59 and 3.06 cm (10 yr), respectively.Moreover, comparing surface air temperature, heat content and sea level rise in SSP1-2.6, we find that, after 2050, the concentrations of COstop increasing or decreasing, surface air temperature has warmed weakly or ground to a halt,but the deep sea is still warming and the thermosteric sea level is still rising. That means that responses of temperature in the deep ocean and thermosteric sea level to the external radiation forcing have long-term hysteretic effects—committed climate warming (Mauritsen and Pincus, 2017). Although the variation trend of thermosteric sea level of the CMIP5 ensemble mean is consistent with that of our model, the result of our model is smaller than that in CMIP5 (Collins et al., 2013).

    In the context of global warming, the water vapor in the air increases and the global precipitation increases. The annual global mean precipitation (Fig.2g) for all four scenarios is projected to increase gradually during 2015–2100.For SSP1-2.6, SSP2-4.5, SSP3-7.0 and SSP5-8.5, the precipitation in 2015 accounts for about 2.3%, 4.4%, 5.1% and 7.2%, respectively, of that in 2100. The precipitation change of the CMIP5 ensemble average by 2100 is less than that simulated by the model (Collins et al., 2013). Based on calculations and analyses, the correlation relationship between global precipitation and global temperature is positive and approximately linear. The changes in atmospheric energy balance lead to the changes in global precipitation mainly through the response of the surface temperature(Andrews et al., 2009; Bala et al., 2010). In addition, the precipitation sensitivity, i.e., the changes in the percentage of global precipitation relative to global temperature, is about 2% °Cunder global warming, which is roughly similar to previous estimates by Held and Soden (2006).

    3.2.Changes in surface air temperature

    Fig.3. Surface air temperature changes in 2080–2100 for SSP1-2.6, SSP2-4.5, SSP3-7.0 and SSP5-8.5 relative to the 20-year mean from 1995 to 2014 in historical runs (units: °C).

    The geographical distributions of surface 2 m air temperature change in 2081–2100 for SSP1-2.6, SSP2-4.5, SSP3-7.0 and SSP5-8.5 relative to the 20-year mean from 1995 to 2014 of the corresponding historical simulation are shown in Fig.3. The geographical patterns of future surface air temperature warming are similar in all four scenarios, and similar to those obtained in previous models of CMIP5, such as the polar amplification of temperature change and stronger warming over the continents (Hansen et al., 2007). Polar amplification refers to amplified surface warming in Arctic latitudes, resulting primarily from positive feedbacks from the retreat of ice and snow. The lack of polar amplification in the Antarctic is relevant to deep ocean mixing, strong heat uptake and the persistence of the vast Antarctic ice sheet.The key feature of warming is that it is larger over land areas than over oceans. The warming over all continents in the SSP1-2.6 scenario (Fig.3a) is four to five times smaller than that in the SSP5-8.5 case (Fig.3d), and two to three times smaller than that in the SSP2-4.5 (Fig.3b) and SSP3-7.0 (Fig.3c) experiments. The North Atlantic is slightly warmer (Figs. 3c and d) and even cooler (Figs. 3a and b),which is called the North Atlantic warming hole (NAWH)(Drijfhout et al., 2012) and is also reflected in the result of the CMIP5 ensemble mean (Collins et al., 2013). Based on previous research, it is likely that the NAWH is also related to the weakening of the Atlantic Meridional Overturning Circulation (AMOC) (Knutti et al., 2004; Timmermann et al.,2007). In addition, the NAWH is also related to warming of the tropical Indian Ocean (TIO). The rainfall related to the warming in the TIO induces a positive NAO-like response with enhanced surface westerly winds over the subpolar North Atlantic, cooling the underlying ocean via stronger turbulent heat flux and southward Ekman transport, resulting in an NAWH pattern (Hu and Fedorov, 2020, personal communication). In summary, large-scale warming patterns are very similar in all four scenarios.

    3.3.Changes in precipitation

    The geographical distributions of projected precipitation changes in 2081–2100 for SSP1-2.6, SSP2-4.5, SSP3-7.0 and SSP5-8.5 relative to the 20-year mean from 1995 to 2014 of the corresponding historical simulation are shown in Fig.4. Similar to the surface 2 m air temperature, the changes in precipitation are smallest in SSP1-2.6 (Fig.4a)and largest in SSP5-8.5 (Fig.4d). Although global water vapor increases with global warming, the changes in precipitation in a warming world are not uniform across regions,which is mainly related to the regional change of atmospheric circulation and the regional increase of water vapor.Precipitation increases along the intertropical convergence zone (ITCZ), which is due to strong moisture convergence mainly associated with atmospheric circulation changes(Seager et al., 2010). The increases of water vapor coupled to the pre-existing circulation will increase the moisture convergence in the ITCZ, where there is already moisture convergence and rainfall (Stevens et al., 2013). Furthermore,through our analysis of meridional winds (not shown), the convergence intensity of meridional winds over the ITCZ region also increases; plus, the ITCZ is situated slightly southward, which should be related to the strength of the simulated meridional wind. By comparing the average precipitation of the four experiments from 2080 to 2100 with the average precipitation of the Global Precipitation Climatology Project monthly precipitation dataset from 1980 to 2014, it can be seen that the precipitation in the central and eastern equatorial Pacific Ocean increases by more than 100% in the four experiments. In addition, precipitation increases in middle and high latitudes on account of enhanced poleward transport of moisture, which is due to increased water vapor in the lower troposphere and enhanced circulation.However, precipitation decreases over much of the subtropics on account of weaker moisture transport, which is because there is a robust decrease in the difference between precipitation-minus-evaporation in these dry regions (Held and Soden, 2006; Seager et al., 2010). Specifically, precipitation increases over most of the Eurasian continent, North America, and much of South America, and precipitation decreases in the oceanic subtropics and from there into the Amazon Basin, Mediterranean, and southern regions of Africa. These features are generally similar to the results of CMIP5 (Collins et al., 2013; Nazarenko et al., 2015).

    Fig.4. Projected precipitation changes averaged between 2080 and 2100 from SSP1-2.6, SSP2-4.5, SSP3-7.0 and SSP5-8.5 simulations relative to the 20-year mean from 1995 to 2014 from historical runs (units: mm d?1).

    3.4.Changes in zonally averaged ocean temperature

    The changes in zonal-mean ocean temperature for 2080–2100 under the SSP1-2.6, SSP2-4.5, SSP3-7.0 and SSP5-8.5 scenarios relative to the 20-year mean from 1995 to 2014 in historical runs are shown in Fig.5. At the end of the 21st century, there is a general warming above 2000 m,while there is little change below 2000 m in most of the ocean. That means most of the extra energy in the ocean is stored above 2000 m. The extra heat is gradually transferred to a greater depth of 2000 m through the processes of mixing and advection. In the deep sea, especially from 2000 m to 3000 m, due to deep ocean mixing, the warming of the Southern Ocean is more obvious. Also, there is a slight cooling in parts of the ocean north of 20°N below 2000 m,which is likely related to the projected decrease in the strength of the AMOC. Due to vertical advection directly from the surface, the descending branch of the AMOC results in a reduction in downward heat transfer and a slight cooling north of 20°N below 2000 m (Banks and Gregory,2006).

    3.5.Summary

    If GHG emissions continue unabated, both the heat content in the ocean and the globally averaged surface air temperature will continue to rise over the 21st century, as expected. For the four scenarios, the projected change in global mean surface air temperature relative to the climatological mean between 1850 and 1900 from historical runs is about 1.97°C, 2.84°C, 3.91°C and 4.66°C, respectively, around the end of the 21st century. Temperature changes will not be regionally uniform, and much larger over land than ocean areas. In particular, there is polar amplification in the Arctic region. Due to this polar amplification, the sea-ice extent in the Arctic, especially in September, will significantly decrease in the coming decades. The global ocean shows continued warming under all four scenarios. The deep ocean responds slowly to the climate change. Most of the extra energy released from human activity is stored in the oceans. At the end of the 21st century, there is projected to be a general warming above 2000 m, while there is little change below 2000 m in most of the ocean. However,warming of the Southern Ocean is also evident below 2000 m. Due to global ocean warming, the global mean thermosteric sea level rises under all four scenarios. Global mean precipitation will also increase with increased global mean surface temperature in the future 21st century. In addition to changes in long-term trends and climatological mean states,the response of internal climate variability, such as ENSO,to global warming is another hot topic. In fact, preliminary analysis shows that ENSO amplitudes in ScenarioMIP are lightly larger than those in the piControl and historical runs,but there are no significant changes in ENSO frequency among these experiments. In future work, detailed dynamical feedback analyses of the ENSO cycle are required to explain the changes in ENSO behavior as a result of global warming.

    Fig.5. Changes in zonal mean temperature in the ocean for 2080–2100 under the SSP1-2.6, SSP2-4.5, SSP3-7.0 and SSP5-8.5 scenarios relative to the 20-year mean from 1995 to 2014 in historical runs (units: °C).

    In general, the projected climate changes from CMIP6 simulations by FGOALS-f3-L are qualitatively consistent with CMIP5 simulations by FGOALS-s2, with only some quantitative results that are slightly different. Based on the abrupt four times COexperiments, the estimated climate sensitivity has been reduced from 4.5°C in FGOALS-s2 to 3.0°C in FGOALS-f3-L due to changes in some subgrid parametrization schemes (He et al., 2020). As a result, the projected global mean surface air temperature by FGOALS-s2 RCP8.5 experiments, about 5.3°C, is much larger than that of about 4.66°C by FGOALS-f3-L SSP5-8.5 experiments.Compared with the spatial distributions of surface air temperature from these two models under the high GHG emissions scenario, the simulated cold anomalies in the North Atlantic and the Southern Ocean, and the simulated warm anomalies of the Southern Ocean in FGOASL-s2, are larger than those in FGOASL-f3-L, which might be partly related to changes in thermohaline circulations in the deep ocean. In addition, under different GHG emission pathways,the maxima of different models appear at different locations (Zhou et al., 2014a). The difference in quantitative results may be resulted from differences in the pathways between CMIP5 and CMIP6, climate sensitivities, and the climatological mean climatology simulated by these two models, but the detailed mechanisms still need to be further analyzed.

    4.Data records and usage notes

    Four datasets of four scenarios (SSP1-2.6, SSP2-4.5,SSP3-7.0 and SSP5-8.5) have been submitted to the ESGF node and are available at https://esgf-node.llnl.gov/search/cmip6/. These datasets have been post-processed by version 3 of the Climate Model Output Rewriter software(CMOR3, https://cmor.llnl.gov/). The dataset format is version 4 of the Network Common Data Form (NetCDF),which is easy to read and written by professional common software such as Climate Data Operators (CDO, https://www.unidata.ucar.edu/software/netcdf/workshops/mostrecent/third_party/CDO.html), NCAR Command Language(NCL, http://www.ncl.ucar.edu), NetCDF Operator (NCO,http://nco.sourceforge.net) and Python (https://www.python.org).

    All variables are uploaded to the ESGF node with their original grid, except for atmospheric variables, in which the original grid in cube–sphere form has been interpolated into a global latitude–longitude grid. The oceanic model outputs are on a tripolar grid with 30 vertical levels. The atmo-spheric model outputs are on a global latitude–longitude grid with 32 vertical levels. Please refer to section 2.1 for horizontal resolution. The 125 variables of the four scenarios,such as SST, sea water salinity, sea water potential temperature, sea surface height above geoid, and so on, are available. In addition, in the ScenarioMIP experiments, a large number of high-frequency outputs are added, such as 3 h, 6 h, etc., which can be directly used to drive regional climate models for dynamic or statistical downscaling. The most frequently used monthly variables are listed in Table 2. Scalar data in the dataset can be interpolated directly, while vector data need to be rotated according to the angles between the original grid and the latitude–longitude grid before interpolation. If users need other related data variables, please feel free to contact us.

    Table 2. Descriptions of common monthly variables in ScenarioMIP datasets of CAS FGOALS-f3-L. Horizontal resolution: 1° (lat) ×1.25° (lon) (atmosphere); 1° × 1° (ocean, sea ice); 0.9° (lat) × 1.25° (lon)(land).

    Table 2. (Continued.)

    Acknowledgements.

    This study was jointly supported by the Strategic Priority Research Program of the Chinese Academy of Sciences (Grant Nos. XDA19060102 and XDB42000000) and the National Natural Science Foundation of China (Grants Nos.41530426 and 91958201).The authors acknowledge the technical support from the National Key Scientific and Technological Infrastructure project “Earth System Science Numerical Simulator Facility ” (EarthLab).

    Disclosure statement

    No potential conflict of interest was reported by the authors.

    Open Access

    This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source,provide a link to the Creative Commons license, and indicate if changes were made.

    亚洲国产精品合色在线| 午夜福利在线在线| 日韩中字成人| 超碰av人人做人人爽久久| 美女cb高潮喷水在线观看| 国产亚洲最大av| 久久精品国产亚洲网站| 日产精品乱码卡一卡2卡三| 一个人看视频在线观看www免费| 校园人妻丝袜中文字幕| 一边亲一边摸免费视频| 99久久中文字幕三级久久日本| 黄片无遮挡物在线观看| 神马国产精品三级电影在线观看| 欧美三级亚洲精品| 大香蕉久久网| 国产免费男女视频| 亚洲成av人片在线播放无| 人人妻人人看人人澡| 高清午夜精品一区二区三区| 国产私拍福利视频在线观看| 日本一二三区视频观看| 精品99又大又爽又粗少妇毛片| 18+在线观看网站| 日本-黄色视频高清免费观看| 久久精品国产自在天天线| 我要看日韩黄色一级片| 国产亚洲一区二区精品| 国产伦一二天堂av在线观看| 亚洲乱码一区二区免费版| 97在线视频观看| 一本久久精品| 观看免费一级毛片| 国产高清不卡午夜福利| 色综合亚洲欧美另类图片| 两个人视频免费观看高清| 免费在线观看成人毛片| 国产精品不卡视频一区二区| 国产亚洲5aaaaa淫片| 九色成人免费人妻av| 插阴视频在线观看视频| 久久久a久久爽久久v久久| 伦理电影大哥的女人| 久久精品国产亚洲av涩爱| 麻豆精品久久久久久蜜桃| h日本视频在线播放| 最近最新中文字幕大全电影3| 99热这里只有精品一区| 麻豆精品久久久久久蜜桃| 色视频www国产| 男人和女人高潮做爰伦理| 午夜老司机福利剧场| 久久6这里有精品| 欧美一区二区亚洲| 久久综合国产亚洲精品| 欧美一区二区国产精品久久精品| 午夜a级毛片| 网址你懂的国产日韩在线| 秋霞在线观看毛片| 男插女下体视频免费在线播放| 欧美成人一区二区免费高清观看| 一级毛片我不卡| 国产淫片久久久久久久久| 亚洲国产色片| 在线免费十八禁| 国产黄色视频一区二区在线观看 | av天堂中文字幕网| 国产精品综合久久久久久久免费| 自拍偷自拍亚洲精品老妇| 亚洲18禁久久av| 一个人免费在线观看电影| 卡戴珊不雅视频在线播放| 99热6这里只有精品| 国产激情偷乱视频一区二区| av播播在线观看一区| 欧美性猛交╳xxx乱大交人| 99热全是精品| 少妇熟女aⅴ在线视频| 久久精品国产自在天天线| 免费搜索国产男女视频| 国产高清三级在线| 亚洲激情五月婷婷啪啪| 日本爱情动作片www.在线观看| 超碰97精品在线观看| 欧美人与善性xxx| 婷婷色综合大香蕉| 午夜福利在线观看吧| 级片在线观看| 国产免费一级a男人的天堂| 久久人人爽人人爽人人片va| 国产视频内射| 国产精品.久久久| 热99在线观看视频| 久久久色成人| av在线亚洲专区| 国产午夜精品久久久久久一区二区三区| 偷拍熟女少妇极品色| 成人鲁丝片一二三区免费| 欧美人与善性xxx| 黄片wwwwww| 91久久精品国产一区二区成人| 午夜精品国产一区二区电影 | 国产精品蜜桃在线观看| 女人久久www免费人成看片 | 国产精品不卡视频一区二区| 久久亚洲精品不卡| 好男人在线观看高清免费视频| 免费人成在线观看视频色| 久久热精品热| 两个人视频免费观看高清| 色尼玛亚洲综合影院| 日本黄色视频三级网站网址| 狂野欧美激情性xxxx在线观看| 国产视频内射| 亚洲色图av天堂| 菩萨蛮人人尽说江南好唐韦庄 | 日韩,欧美,国产一区二区三区 | 中文字幕久久专区| 一个人看的www免费观看视频| 国产精品综合久久久久久久免费| 国产亚洲最大av| 51国产日韩欧美| 人妻系列 视频| 男人舔女人下体高潮全视频| 中文欧美无线码| 国产成人精品婷婷| 国产精品一区二区三区四区久久| 亚洲av男天堂| 国产男人的电影天堂91| 综合色av麻豆| 一边摸一边抽搐一进一小说| 建设人人有责人人尽责人人享有的 | 免费av观看视频| 黑人高潮一二区| 一区二区三区免费毛片| 日本三级黄在线观看| 亚洲人成网站在线观看播放| 成人毛片a级毛片在线播放| 亚洲自偷自拍三级| av卡一久久| 亚洲精品456在线播放app| 精品久久久久久久人妻蜜臀av| 日韩av在线免费看完整版不卡| 亚洲欧洲日产国产| 免费在线观看成人毛片| 国产伦理片在线播放av一区| 婷婷色综合大香蕉| 亚洲高清免费不卡视频| 亚洲美女搞黄在线观看| 91狼人影院| 国产三级中文精品| 七月丁香在线播放| 久久久久久国产a免费观看| 国产片特级美女逼逼视频| 亚洲久久久久久中文字幕| 中文字幕制服av| 免费电影在线观看免费观看| 成人国产麻豆网| 有码 亚洲区| 亚洲综合精品二区| 中文字幕免费在线视频6| 尾随美女入室| ponron亚洲| 在线免费观看的www视频| 91精品伊人久久大香线蕉| 99热这里只有是精品在线观看| 亚洲三级黄色毛片| 日本免费在线观看一区| 久久亚洲精品不卡| 久久精品久久精品一区二区三区| 尤物成人国产欧美一区二区三区| 汤姆久久久久久久影院中文字幕 | 久久精品综合一区二区三区| 亚洲美女搞黄在线观看| 日产精品乱码卡一卡2卡三| 男女视频在线观看网站免费| 国产伦精品一区二区三区四那| 国产一区二区三区av在线| 亚洲五月天丁香| 亚洲三级黄色毛片| 亚洲美女视频黄频| 国产黄色视频一区二区在线观看 | 日本欧美国产在线视频| 国产成人freesex在线| 亚洲av男天堂| 久热久热在线精品观看| 日韩欧美 国产精品| 超碰av人人做人人爽久久| 亚洲天堂国产精品一区在线| 日韩av不卡免费在线播放| 岛国在线免费视频观看| 大又大粗又爽又黄少妇毛片口| 深夜a级毛片| 性色avwww在线观看| 日韩欧美在线乱码| 黑人高潮一二区| 成人鲁丝片一二三区免费| 国产精品久久久久久久电影| 天美传媒精品一区二区| 看片在线看免费视频| 日本午夜av视频| 午夜免费激情av| 国产成人午夜福利电影在线观看| 亚洲精品色激情综合| 久久久国产成人精品二区| 一个人看的www免费观看视频| 国产成人福利小说| 乱码一卡2卡4卡精品| 国产三级中文精品| 日本欧美国产在线视频| 免费av不卡在线播放| 91久久精品国产一区二区三区| 午夜精品国产一区二区电影 | 久久精品熟女亚洲av麻豆精品 | 97超视频在线观看视频| 欧美极品一区二区三区四区| 日韩av不卡免费在线播放| 嫩草影院入口| 久久久久久久久久成人| 国产黄a三级三级三级人| 久久久久久久国产电影| 在线观看av片永久免费下载| 最后的刺客免费高清国语| 亚洲无线观看免费| 色哟哟·www| 中文字幕亚洲精品专区| 人人妻人人澡欧美一区二区| 成人性生交大片免费视频hd| 黄色日韩在线| av在线播放精品| 啦啦啦观看免费观看视频高清| 精品久久久久久久久亚洲| 欧美色视频一区免费| 一个人看视频在线观看www免费| 五月玫瑰六月丁香| 亚洲美女搞黄在线观看| 欧美又色又爽又黄视频| 麻豆成人午夜福利视频| 精品不卡国产一区二区三区| 久久99热这里只频精品6学生 | 久久久成人免费电影| 日韩一区二区三区影片| 国产精品综合久久久久久久免费| 男人的好看免费观看在线视频| 亚洲无线观看免费| 欧美一区二区国产精品久久精品| 久久久精品大字幕| 免费电影在线观看免费观看| 国产精品蜜桃在线观看| 国产三级在线视频| 亚洲人成网站在线观看播放| 国产精品,欧美在线| 麻豆精品久久久久久蜜桃| 美女脱内裤让男人舔精品视频| 一级二级三级毛片免费看| 久久人妻av系列| 免费观看a级毛片全部| 国产老妇伦熟女老妇高清| 黄色一级大片看看| 超碰av人人做人人爽久久| 精品国内亚洲2022精品成人| 一区二区三区高清视频在线| .国产精品久久| 大香蕉久久网| 菩萨蛮人人尽说江南好唐韦庄 | 国产亚洲一区二区精品| 久久久a久久爽久久v久久| 国产午夜精品论理片| 精品一区二区三区人妻视频| 国内揄拍国产精品人妻在线| 久久久久性生活片| 成年女人看的毛片在线观看| 亚洲av熟女| 亚洲av中文字字幕乱码综合| 亚州av有码| 亚洲欧美精品综合久久99| 日韩人妻高清精品专区| 男人和女人高潮做爰伦理| 一个人看的www免费观看视频| 青春草亚洲视频在线观看| av在线老鸭窝| av免费观看日本| 在线免费十八禁| 国产老妇女一区| 欧美性猛交╳xxx乱大交人| 大又大粗又爽又黄少妇毛片口| 国产精品av视频在线免费观看| 久久久久性生活片| 日日啪夜夜撸| 一夜夜www| 少妇的逼好多水| 亚洲成人精品中文字幕电影| 国产亚洲最大av| 在线免费十八禁| 丝袜喷水一区| 日韩av不卡免费在线播放| 亚洲欧美清纯卡通| 国产精品乱码一区二三区的特点| 国产黄a三级三级三级人| 中文字幕精品亚洲无线码一区| 国产单亲对白刺激| 卡戴珊不雅视频在线播放| 三级毛片av免费| 一级爰片在线观看| 日韩欧美 国产精品| 美女黄网站色视频| 亚洲va在线va天堂va国产| 国产一区有黄有色的免费视频 | 日韩成人伦理影院| 国产国拍精品亚洲av在线观看| 国产久久久一区二区三区| 黄片无遮挡物在线观看| 建设人人有责人人尽责人人享有的 | 亚洲自拍偷在线| 在现免费观看毛片| 白带黄色成豆腐渣| 两个人视频免费观看高清| 91精品一卡2卡3卡4卡| 国产成人福利小说| 啦啦啦韩国在线观看视频| 国产成人91sexporn| 国产视频内射| 亚洲精品自拍成人| 乱系列少妇在线播放| 成年女人看的毛片在线观看| 欧美成人一区二区免费高清观看| av天堂中文字幕网| 有码 亚洲区| 最近最新中文字幕大全电影3| 国产精品av视频在线免费观看| 日韩一区二区三区影片| 欧美日本视频| 国产私拍福利视频在线观看| 亚洲国产精品国产精品| 久久这里只有精品中国| 成人亚洲精品av一区二区| 蜜臀久久99精品久久宅男| 亚洲av免费高清在线观看| 日韩欧美精品v在线| 亚洲高清免费不卡视频| 尾随美女入室| 听说在线观看完整版免费高清| 国产精品电影一区二区三区| 天堂√8在线中文| 国产精品麻豆人妻色哟哟久久 | 少妇猛男粗大的猛烈进出视频 | 国产成人福利小说| 我的老师免费观看完整版| 亚洲av福利一区| 欧美xxxx黑人xx丫x性爽| 国产免费又黄又爽又色| 免费看a级黄色片| 日韩欧美精品免费久久| 国产一区有黄有色的免费视频 | 国产成人精品久久久久久| 精品久久久久久电影网 | 我的老师免费观看完整版| 日日啪夜夜撸| 变态另类丝袜制服| 建设人人有责人人尽责人人享有的 | 成人av在线播放网站| 特级一级黄色大片| 五月伊人婷婷丁香| 久久久国产成人免费| 国产精品国产高清国产av| 亚洲内射少妇av| 精品无人区乱码1区二区| 午夜激情欧美在线| 搡老妇女老女人老熟妇| 国产人妻一区二区三区在| 成人亚洲精品av一区二区| 中文字幕av成人在线电影| 久久6这里有精品| av专区在线播放| 亚洲欧美成人综合另类久久久 | 国产精品电影一区二区三区| 少妇人妻一区二区三区视频| 97热精品久久久久久| 中国美白少妇内射xxxbb| 高清日韩中文字幕在线| 简卡轻食公司| 国产白丝娇喘喷水9色精品| av免费在线看不卡| 亚洲欧美清纯卡通| av在线亚洲专区| 久久久久久久久久黄片| 欧美丝袜亚洲另类| 高清午夜精品一区二区三区| 淫秽高清视频在线观看| 日韩一区二区三区影片| 黄色一级大片看看| 成人毛片a级毛片在线播放| 亚洲欧洲国产日韩| 最近2019中文字幕mv第一页| 久久6这里有精品| 美女高潮的动态| 99热精品在线国产| 最近中文字幕2019免费版| 国产私拍福利视频在线观看| av国产久精品久网站免费入址| 精品人妻熟女av久视频| 亚洲人成网站高清观看| 国产精品一区二区在线观看99 | 国产精品一区二区三区四区免费观看| 寂寞人妻少妇视频99o| 边亲边吃奶的免费视频| 3wmmmm亚洲av在线观看| 美女xxoo啪啪120秒动态图| 久久久国产成人精品二区| 国产色爽女视频免费观看| 女人久久www免费人成看片 | 精品一区二区免费观看| 国产一区二区亚洲精品在线观看| 青春草国产在线视频| 在线天堂最新版资源| 丰满人妻一区二区三区视频av| 女人被狂操c到高潮| 欧美日本视频| 国产国拍精品亚洲av在线观看| 少妇人妻精品综合一区二区| 久久国产乱子免费精品| 久久久久网色| 日日啪夜夜撸| 中文字幕制服av| 欧美又色又爽又黄视频| 女的被弄到高潮叫床怎么办| 黄色配什么色好看| 国产中年淑女户外野战色| 久久久久久久久久成人| 久久精品久久精品一区二区三区| 亚洲成av人片在线播放无| 日本免费a在线| 天堂影院成人在线观看| 国产亚洲91精品色在线| 久久精品国产鲁丝片午夜精品| 爱豆传媒免费全集在线观看| 免费看av在线观看网站| 又粗又爽又猛毛片免费看| 三级国产精品片| 亚洲精品aⅴ在线观看| 精品无人区乱码1区二区| 在线播放无遮挡| 久久久久久伊人网av| 亚洲美女搞黄在线观看| 自拍偷自拍亚洲精品老妇| 午夜福利视频1000在线观看| av在线亚洲专区| 国产成人精品婷婷| 国产一级毛片在线| 日韩制服骚丝袜av| 91精品国产九色| 亚洲国产精品专区欧美| 午夜激情福利司机影院| 亚洲国产精品sss在线观看| 国产一级毛片在线| 成人三级黄色视频| 久久精品国产鲁丝片午夜精品| 少妇高潮的动态图| 两性午夜刺激爽爽歪歪视频在线观看| 一区二区三区乱码不卡18| 欧美高清性xxxxhd video| 国语自产精品视频在线第100页| 丰满人妻一区二区三区视频av| 精品欧美国产一区二区三| 午夜视频国产福利| 日本免费一区二区三区高清不卡| 狂野欧美激情性xxxx在线观看| 亚洲av不卡在线观看| 黄色日韩在线| 少妇的逼好多水| 欧美成人免费av一区二区三区| 国产免费福利视频在线观看| 亚洲精品影视一区二区三区av| 青春草视频在线免费观看| 中文字幕av成人在线电影| 亚洲国产精品成人综合色| 特级一级黄色大片| 99热网站在线观看| 亚洲aⅴ乱码一区二区在线播放| 日韩成人av中文字幕在线观看| 嫩草影院入口| 如何舔出高潮| 亚洲中文字幕日韩| a级毛片免费高清观看在线播放| 美女高潮的动态| 午夜福利在线观看免费完整高清在| 欧美性猛交╳xxx乱大交人| 啦啦啦观看免费观看视频高清| 精品国产三级普通话版| 亚洲在线观看片| 亚洲国产欧美人成| 成人欧美大片| 99国产精品一区二区蜜桃av| 久久精品国产亚洲网站| av在线观看视频网站免费| 老司机影院成人| 国产精品无大码| 老师上课跳d突然被开到最大视频| 九九久久精品国产亚洲av麻豆| 日韩视频在线欧美| 日韩欧美三级三区| 色5月婷婷丁香| 久久韩国三级中文字幕| 看十八女毛片水多多多| 久久久色成人| 亚洲国产高清在线一区二区三| 成年女人永久免费观看视频| 久久鲁丝午夜福利片| 嫩草影院新地址| 亚洲第一区二区三区不卡| 岛国在线免费视频观看| 国产精品一区二区性色av| 成人二区视频| 亚洲成人精品中文字幕电影| 国产 一区 欧美 日韩| 亚洲精品色激情综合| 亚洲无线观看免费| 综合色av麻豆| 国产精品1区2区在线观看.| 精品免费久久久久久久清纯| 国产成人精品久久久久久| 国产精品熟女久久久久浪| 只有这里有精品99| 成人鲁丝片一二三区免费| 久久久久久伊人网av| 久久久国产成人精品二区| 爱豆传媒免费全集在线观看| 欧美人与善性xxx| 国产黄色小视频在线观看| 欧美不卡视频在线免费观看| 国产精品久久视频播放| 黄色一级大片看看| 能在线免费观看的黄片| 日韩av在线免费看完整版不卡| 亚洲av不卡在线观看| 久久精品国产自在天天线| 99久国产av精品国产电影| 国产 一区精品| 人妻少妇偷人精品九色| 少妇丰满av| 中文字幕亚洲精品专区| 男女边吃奶边做爰视频| 欧美高清性xxxxhd video| 欧美一区二区国产精品久久精品| 日韩欧美在线乱码| 一本久久精品| 亚洲av日韩在线播放| 晚上一个人看的免费电影| 亚洲国产精品sss在线观看| 午夜精品一区二区三区免费看| 超碰av人人做人人爽久久| 国内精品一区二区在线观看| 国产人妻一区二区三区在| 午夜福利网站1000一区二区三区| 中文字幕免费在线视频6| 成年女人永久免费观看视频| kizo精华| 日本一二三区视频观看| 1000部很黄的大片| 亚洲国产成人一精品久久久| 国产精品av视频在线免费观看| 久久这里有精品视频免费| 日本猛色少妇xxxxx猛交久久| 成人午夜精彩视频在线观看| 成人特级av手机在线观看| 欧美日韩综合久久久久久| 久久国内精品自在自线图片| 国产午夜精品久久久久久一区二区三区| 久久精品国产99精品国产亚洲性色| 亚洲精品乱久久久久久| 精品国内亚洲2022精品成人| 欧美zozozo另类| 22中文网久久字幕| 非洲黑人性xxxx精品又粗又长| 国产精品久久久久久精品电影| 国产在线男女| 搡老妇女老女人老熟妇| 国产爱豆传媒在线观看| 久久精品夜色国产| av免费在线看不卡| 大又大粗又爽又黄少妇毛片口| 一个人观看的视频www高清免费观看| 国产高清不卡午夜福利| 国产69精品久久久久777片| 免费观看在线日韩| 亚洲国产精品成人久久小说| ponron亚洲| 免费看av在线观看网站| 少妇的逼好多水| 在线观看66精品国产| 日韩亚洲欧美综合| 亚洲精品乱码久久久久久按摩| 国产大屁股一区二区在线视频| 久久99精品国语久久久| 午夜精品国产一区二区电影 | 免费大片18禁| 久久人人爽人人片av| 久久亚洲国产成人精品v| 精品久久国产蜜桃| 亚洲国产成人一精品久久久| 国内精品美女久久久久久| 日本欧美国产在线视频| 在线a可以看的网站| 国产亚洲av嫩草精品影院| 国产白丝娇喘喷水9色精品| 最后的刺客免费高清国语| 在线观看66精品国产| 狠狠狠狠99中文字幕| 日韩精品青青久久久久久| 成人性生交大片免费视频hd| 免费av毛片视频| 午夜久久久久精精品| 久久综合国产亚洲精品| 国产精品一区二区三区四区免费观看| 老司机影院成人|