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

    CAS-ESM2.0 Model Datasets for the CMIP6 Flux-Anomaly-Forced Model Intercomparison Project (FAFMIP)

    2021-02-26 08:22:38JiangboJINHeZHANGXiaoDONGHailongLIUMinghuaZHANGXinGAOJuanxiongHEZhaoyangCHAIQingcunZENGGuangqingZHOUZhaohuiLINYiYUPengfeiLINRuxuLIANYongqiangYUMirongSONGandDonglingZHANG
    Advances in Atmospheric Sciences 2021年2期

    Jiangbo JIN, He ZHANG, Xiao DONG*, Hailong LIU, Minghua ZHANG, Xin GAO, Juanxiong HE,Zhaoyang CHAI, Qingcun ZENG, Guangqing ZHOU, Zhaohui LIN, Yi YU, Pengfei LIN,Ruxu LIAN, Yongqiang YU, Mirong SONG, and Dongling ZHANG

    1International Center for Climate and Environment Sciences, Institute of Atmospheric Physics,Chinese Academy of Sciences, Beijing 100029, China

    2State Key Laboratory of Satellite Ocean Environment Dynamics, Second Institute of Oceanography,Ministry of Natural Resources, Hangzhou 310012, China

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

    4College of Earth and Planetary Sciences, University of Chinese Academy of Sciences, Beijing 100049, China

    5Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100029, China

    ABSTRACT

    Key words:CAS-ESM2.0,CMIP6,FAFMIP,AMOC,ocean heat uptake,dynamic sea level change

    1.Background

    Currently, climate models exhibit large spread in their future projections of global climate. For instance, the coupled climate models of phase 5 of the Coupled Model Intercomparison Project (CMIP5) project a global mean sea level rise (SLR) of 0.32–0.63 m by 2081–2100 under the RCP4.5 scenario (Yin, 2012; IPCC, 2013), and the thermosteric contribution due to seawater expansion can be up to 30% of the global mean SLR and up to 100% in the North Atlantic region (Gregory et al., 2016). The spread of ocean climate changes may be associated with the spread of the surface flux perturbation responses to increased atmospheric COconcentrations (Bouttes and Gregory, 2014). However,the spatial patterns and magnitudes of the sea surface flux perturbations are largely model-dependent, which makes it complicated and difficult to directly compare ocean responses in coupled models. Therefore, to understand the ocean responses to increased COlevels and reduce the uncertainty of ocean state estimations, a flux-anomaly-forced model intercomparison project (FAFMIP), one of 23 Model Intercomparison Projects (MIPs) in CMIP6 (Eyring et al.,2016), has been proposed to compare the different model responses to the same prescribed surface flux perturbations(Gregory et al., 2016). According to Gregory et al. (2016),around ten modeling groups plan to participate in FAFMIP.At the time of this report, FAFMIP experiment outputs from nine models are available on the Earth System Grid Federation (ESGF) data server for CMIP6.

    In China, thirteen models from nine institutes or universities in mainland China and Taiwan are participating in CMIP6 (Zhou et al., 2019), and are using thirteen different models. Among these models, CAS-ESM2.0 is an earth system model (ESM) from the Institute of Atmospheric Physics, Chinese Academy of Sciences (IAP/CAS), which has a long history for the development of climate models. Since the 1980s, scientists at IAP/CAS have been devoted to developing climate models that include an atmospheric general circulation model (AGCM), ocean general circulation model(OGCM), and land surface model (LSM) (Zhang and Zeng,1988; Zeng et al., 1989; Dai and Zeng, 1997), which have been used for short-term climate predictions in China and have participated in each phase of CMIP since the 1990s(Zhou et al., 2020). CAS-ESM2.0 is the successor to CASESM1.0, which was released to the public in 2015. Compared to CAS-ESM1.0, CAS-ESM2.0 considerably modified codes and parameterization schemes in most components, as well as substantially improved performance in its climate simulations (Zhang et al., 2020; Zhou et al., 2020).The FAFMIP experiments were completed using CASESM2.0 by closely following the standard protocol of Gregory et al. (2016); the corresponding data have been uploaded to the ESGF data server for CMIP6 users to download and can be found at https://esgf-node.llnl.gov/projects/cmip6/.

    The purpose of this paper is to present the preliminary results of CAS-ESM2.0 for simulations of ocean changes in FAFMIP experiments for users of the CMIP6 datasets. The remainder of the paper is structured as follows: Section 2 describes the model and experiments; section 3 provides a basic technical validation of the CAS-ESM2.0 experiments;and sections 4 and 5 offer usage notes and a summary,respectively.

    2.Model and experimental design

    2.1.Model

    CAS-ESM2.0 is a fully coupled ESM that was developed at IAP/CAS. It is based on IAP AGCM version 5 and is derived from the LASG/IAP Climate System Ocean Model(LICOM2), Beijing Normal University/IAP Common Land Model (CoLM), Los Alamos Sea Ice Model (CICE version 4), and the Weather Research and Forecasting (WRF)model. The infrastructure of CESM Coupler 7 is adopted to ensure that the components are coupled together. Additional components in CAS-ESM2.0 include the IAP Vegetation Dynamics Model and IAP fire model, which are embedded within the land model; the IAP ocean biogeochemistry model, which is embedded within the ocean model; an atmospheric aerosol and chemistry model; and various emission models. However, the climate components, which include atmosphere, ocean, land and sea ice, are the main focus in the FAFMIP experiments. A schematic diagram of the CAS-ESM2.0 framework is shown in Fig.1.

    The atmospheric component of CAS-ESM2.0 is IAP AGCM5.0, which is the fifth-generation AGCM developed by IAP. The IAP AGCM is a global gridpoint model that uses the finite-difference scheme with terrain-following for the σ-coordinate. The horizontal resolution of IAP AGCM5 is approximately 1.4° latitude × 1.4° longitude and has 35 vertical levels with a model top at 2.2 hPa. The physical parameterizations in IAP AGCM5 have been largely updated and include new convection and cloud schemes and many modifications that were made to other parameterization schemes compared to its previous version, IAP AGCM4.1, in which the physical parameterizations were mostly taken from CAM5(Zhang et al., 2020).

    Fig.1. Schematic diagram of the CAS-ESM2.0 framework.

    The OGCM in CAS-ESM2.0 is based on LICOM2.0(Liu et al., 2012). The model domain is located between 78.5°S and 87.5°N with a 1° zonal resolution. The meridional resolution is refined to 0.5° between 10°S and 10°N and is increased gradually from 0.5° to 1° between 10° and 20°. There are 30 levels in the vertical direction with 10 m per layer in the upper 150 m. Based on the original version of LICOM2.0, key modifications have been made: (1) a new sea surface salinity boundary condition was introduced that is based on the physical process of air–sea flux exchange at the actual sea–air interface (Jin et al., 2017); (2)intra-daily air–sea interactions are resolved by coupling the atmospheric and oceanic model components once every 2 h;and (3) a new formulation of the turbulent air–sea fluxes (Fairall et al., 2003) was introduced. The ocean component of CAS-ESM2.0 has also participated the Ocean Model Intercomparison Project Phase 1 (OMIP1) experiments for CMIP6 (Dong et al., 2020).

    The sea-ice model is the improved Los Alamos Sea Ice Model (version 4.0) (Hunke and Lipscomb, 2008) and uses the same grid as the oceanic model. This model solves the dynamic and thermodynamic equations for five ice-thickness categories, with one snow layer and four ice layers. For the dynamic component, we use elastic–viscous–plastic rheology (Hunke and Dukowicz 1997), the mechanical redistribution scheme (Lipscomb et al., 2007), and the incremental remapping advection scheme (Lipscomb and Hunke, 2004).For the thermodynamic component, we use parameterization with a more realistic sea-ice salinity budget (Liu, 2010)and an improved CCSM3 radiation scheme for albedo and radiative fluxes at the sea-ice surface by incorporating an explicit melt pond and its impact on the albedo.

    The land component of CAS-ESM2.0 is CoLM (Dai et al., 2003). The initial version of CoLM was adopted as the Community Land Model (CLM) for use with the Community Climate System Model (CCSM; Collins et al.,2006), which was later adopted as the land component for the Beijing Normal University Earth System Model (BNUESM; Ji et al., 2014). The improved version of CoLM has been adopted in CAS-ESM2, e.g., an improved two-stream approximation model of radiation transfer of the canopy, a photosynthesis–stomatal conductance model for sunlit and shaded leaves and for simultaneous transfer of COand water vapor into and out of the leaves, and the embedded IAP Dynamic Global Vegetation Model (IAP-DGVM) (e.g.,Dai et al., 2004; Zeng et al., 2014).

    2.2.Experimental design

    The FAFMIP experiments are flux-anomaly-forced experiments and have been conducted using CAS-ESM2.0.A detailed description of the experimental design of the FAFMIP experiments was supplied by Gregory et al.(2016). FAFMIP presents the prescribed surface flux perturbations of wind stress, heat and fresh water for the doubled COconcentrations (Figs. 2a–c), which are derived from the ensemble mean differences between years 61–80 of 13 CMIP5 CGCM 1pctCOexperiments and all years of piControl control simulations and held seasonal cycles. Seven experiments (namely, faf-stress, faf-heat, faf-water, faf-heat-NA50pct, faf-all, faf-passiveheat, and faf-heat-NA0pct)were conducted using CAS-ESM2.0. Among the experiments, the former four experiments are Tier1, while the latter three are Tier2. The corresponding experimental designs are briefly described as follows:

    Fig.2. Spatial patterns of the surface flux perturbations of (a) momentum (units: 10?3 Pa; colors indicate magnitude of the vector, arrows indicate direction), (b) heat (units: W m?2), and (c) water (units: 10?6 kg m?2 s?1) as added atmospheric forcing in FAFMIP experiments. (d) Spatial pattern of change in the surface heat flux into seawater in the time mean of the final decade of the faf-heat experiment relative to the control (units: W m?2).

    (1) In the faf-stress experiment, a perturbation occurs in the surface eastward (zonal) and northward (meridional)momentum fluxes, i.e., wind stress is added to the momentum balance of the seawater surface. Figure 2a shows the annual mean surface momentum flux perturbations supplied by FAFMIP. The prominent feature is the westerly wind stress in the Southern Ocean, for which the largest positive change exceeds 0.02 Pa.

    (3) In the faf-water experiment, a perturbation of the freshwater flux is added to the sea water surface. Figure 2c shows the annual mean surface water flux perturbation. The pattern is dominated by precipitation changes and the positive center is located in the western equatorial Pacific and in the midlatitudes of the two hemispheres. In contrast, negative water flux occurs in the subtropics.

    (4) The design of the faf-heat-NA50pct experiment is exactly the same as that of the faf-heat experiment except that within a portion of the North Atlantic Ocean, the heat flux perturbation is multiplied by 0.5, as was proposed at the FAFMIP meeting in April 2019. The purpose of the experiment is to produce an Atlantic Meridional Overturning Circulation (AMOC) weakening that is similar to that in 1pctCO2, whereas the original faf-heat experiment provides greater weakening because of the redistribution feedback.

    (5) In the faf-all experiment, the surface flux perturbations of momentum, heat and freshwater are simultaneously imposed on seawater. Method B is also applied to compute the surface flux, which is the same method as in the faf-heat experiment. The experiment is used to quantify the nonlinearities of the effects of the three perturbations.

    (6) In the faf-passiveheat experiment, the heat flux perturbation is applied instead to the passive tracer Tinitialized to zero, which can be used to diagnose the effect of added heat on ocean temperatures due to ocean circulation changes through a comparison of faf-passiveheat with fafheat. The experiment is equivalent to the piControl experiment. The results of this experiment are used as a control experiment in this paper.

    (7) The faf-heat-NA0pct experiment is similar to fafheat-NA50pct but with a zero perturbation in the North Atlantic region, as was also proposed at the FAFMIP meeting in April 2019.

    It is noteworthy that all perturbations are imposed on the seawater surface and do not directly affect the sea ice, to maintain the preindustrial state of the sea ice. The initial conditions of all experiments are the same, and are derived from the end of a 300-year spin up of the piControl experiment of CAS-ESM2.0. Further details of the implementation of each of the experiments can be found at the following websites: http://www.fafmip.org and http://www.fafmip.org/design_apr19.html.

    3.Technical validation

    Surface heat flux perturbation plays a dominant role in the weakening of the AMOC and increases the OHU and SLR among the surface flux perturbations (Rahmstorf and Ganopolski, 1999; Gregory et al., 2005; Bouttes and Gregory, 2014; Gregory et al., 2016). Therefore, minimization of the effect of the redistributed heat flux change ?Q due to circulation changes is very important in the faf-heat and faf-all experiments. Figure 2d shows the change in the surface ?Q of the faf-heat experiment relative to that of the control. The spatial pattern is consistent with the multimodel mean results shown in Gregory et al. (2016). The most significant positive change is observed in the North Atlantic and a negative flux change can be observed in the equatorial Pacific and Atlantic, western boundary in the North Pacific and around the Aleutian Islands, and along the western and eastern coasts of the South American continent. The global ocean area mean ?Q in the final decade is 0.15 W m, which is slightly greater than the multimodel mean result (0.07 W m) presented in Gregory et al.(2016). The global ocean area mean heat flux perturbation is 1.86 W m, and thus 108% of the heat flux perturbation is added to the ocean, which results in a greater increase in the ocean temperature than that in faf-passiveheat (Fig.3c).

    3.1.Time series of changes

    Fig.3. Annual mean time series of the (a) surface air temperature change (units: °C), (b) net heat flux change (units:W m?2), (c) volume mean ocean temperature change (units: °C), and (d) maximum transport change of the AMOC(units: Sv), in six experiments: faf-heat (red line), faf-stress (black line), faf-water (blue line), faf-all (dark green line), faf-heat-NA50pct (orange line) and faf-heat-NA0pct (green line). The dashed black line indicates the fafpassiveheat experiment.

    Figure 3 shows the time series of several key metrics in the experiments following Gregory et al. (2016), which include the global mean surface air temperature change(Fig.3a), net surface heat flux change (Fig.3b), ocean volume mean temperature change (Fig.3c), and maximum AMOC transport change (Fig.3d). From the results of the surface air temperature change, faf-stress varies up and down slightly above the zero line, whereas the faf-heat, faf-water and faf-all experiments all show surface cooling with a magnitude within ?0.5°C. Our results are similar to the model performance in Gregory et al. (2016, top row of their Fig.5),especially regarding the consistency of the results simulated by GFDL-ESM2M. The magnitudes of air temperature change in faf-heat-NA50pct and faf-heat-NA0pct are notably smaller compared to those in faf-heat due to the reduced heat flux perturbation in the North Atlantic region.These results show a decrease in the surface air temperature change in the faf-heat and faf-all experiments, which is mainly associated with the surface heat flux perturbation,especially for the North Atlantic region. The imposed heat flux perturbation in the North Atlantic region could reduce heat transport by AMOC in the heat flux perturbation experiments, which would further induce cooling of air temperatures in the North Atlantic region and for similar-latitude regions of Eurasia and North America (Stouffer et al.,2006). In faf-water, cooling of the surface air temperature is mainly due to the suppression of upward heat transport caused by fresher surface salinity (Gregory et al. 2016).

    As shown in Fig.3b, the net surface heat flux change is positive (~2 W m) in faf-heat and faf-all, whereas in fafstress and faf-water the net surface heat flux change oscillates near the zero line, which indicates that the net heat flux in faf-all is mainly attributable to the surface heat flux change. In our experiments, the global mean heat flux over the last 10 years is 2.01 (2.15) W min faf-heat (faf-all),and the global mean prescribed heat flux perturbations in faf-heat and faf-all are 1.86 W m. Compared with faf-heat and faf-all, the positive heat flux changes in faf-heat-NA50pct and faf-heat-NA0pct are suppressed due to the experimental designs, and the global mean surface heat flux changes in the final decade of the two experiments are 1.55 W mand 1.28 W m, respectively.

    In Fig.3c, the ocean volume mean temperatures increase continuously in the faf-heat and faf-all experiments; however, the faf-water and faf-stress experiments both show minor changes and the faf-water experiment has a slightly larger magnitude of increase than the faf-stress experiment. This feature is highly consistent with the results of Gregory et al. (2016, second row of their Fig.5). At the end of the 70th year of integration, the temperature change in faf-heat (faf-all) is approximately 0.30°C(0.32°C), with the same magnitude as MPI-ESM-LR in Gregory et al. (2016), which is slightly larger than the T(0.28°C) in faf-passiveheat forced only by the prescribed heat flux perturbation. The larger increases in temperature changes in faf-heat and faf-all compared to Tin faf-passiveheat are due to the warming effects of ?Q (Fig.3a). Similar to the air temperature and surface heat flux changes, the variations of ocean volume-mean temperatures in the faf-heat-NA50pct and faf-heat-NA0pct experiments are consistent with those in the faf-heat and faf-all experiments but with reduced magnitudes. The positive ocean volume-mean temperature anomalies at the end of the integrated 70th year reach 0.24°C and 0.20°C in faf-heat-NA50pct and faf-heat-NA0pct, respectively.

    The strength of AMOC is decreased in the faf-heat and faf-all experiments with similar magnitudes (approximately?10 Sv) (Fig.3d). Our results falls in the uncertainty range(e.g., 6–12 Sv) of the models in Gregory et al. (2016, third row of their Fig.5), with a magnitude similar to CanESM2.In the first 30 years in the faf-heat-NA50pct experiment, the reduction of AMOC is nearly consistent with those in the faf-all and faf-heat experiments. After 30 years, the difference between faf-heat-NA50pct and faf-heat/faf-all increases with an ultimate reduction in AMOC to approximately half of that in the faf-all and faf-heat experiments at the end of the 70th year. In the faf-heat-NA0pct experiment,the AMOC change is largely suppressed (within ?4 Sv) due to the absence of the surface heat flux perturbation imposed in the North Atlantic region. Compared with the faf-all and faf-heat experiments, no significant AMOC changes occur in the faf-stress and faf-water experiments (Figs. 3c and d).These time series of the global mean variables further demonstrate that the surface heat flux perturbation plays a dominant role in ocean climate change in response to COforcing among the surface flux perturbations.

    3.2.Spatial patterns of changes

    Changes in ocean heat content (OHC) and dynamic sea level (DSL) are also important metrics for evaluating ocean climate change, which is an essential scientific issue that the community aims to resolve by conducting the FAFMIP experiments. We show the spatial patterns of the DSL and OHC changes in the time means of the final decade of the FAFMIP experiments in Fig.5 and Fig.6, respectively. The greatest local changes in DSL and in OHC occur in faf-heat among the surface flux perturbation experiments (Fig.5 and Fig.6). The most striking feature in faf-all is the three dipole patterns of DSL change: one is located in the North Atlantic (a positive change at approximately 40°N and a negative change southward); the second is located across the Antarctic Circumpolar Circulation (ACC, a positive change to the north and negative change to the south); and the last is over the North Pacific with an opposite dipole pattern of passive change in subtropical areas and a negative change to the north, which are highly similar to the model-mean results in Gregory et al. (2016, their Fig.8). The North Atlantic dipole pattern is primarily due to the heat flux perturbation(Figs. 5a and d), which is consistent with a larger increase in OHC north of 40°N in the Atlantic Ocean in faf-heat and faf-all (Figs. 6a and d). Additionally, the water flux perturbation can partially reinforce the North Atlantic dipole due to the halosteric contribution because an increase in OHC occurs in most areas of the North Atlantic in faf-water(Fig.6b), although the AMOC does not change substantially (Fig.3b). The dipole pattern across the ACC results from the combined effects of the heat and wind stress perturbations and the water flux perturbations provide a counteracting contribution (Figs. 5b–d). The contributions of the heat,water and wind stress perturbations to the DSL changes across the ACC are mainly thermosteric instead of halosteric, as is inferred from the distributions of OHC change (Figs. 6a–d). The opposite dipole pattern over the North Pacific is due to the combined effects of surface heat flux and freshwater flux perturbations (Figs. 6a, b and d).

    Fig.4. Latitude–depth plots of the zonal mean ocean temperature change (units: °C) in the time mean of the final decade in FAFMIP minus that in the control experiment, in six experiments: (a) faf-heat, (b) faf-water, (c) faf-stress,(d) faf-all, (e) faf-heat-NA50pct, and (f) faf-heat-NA0pct. In addition, (g) and (h) represent added and redistributed portions in the faf-heat experiment, respectively.

    For the two new experiments (i.e., faf-heat-NA50pct and faf-heat-NA0pct), although the dipole pattern across the ACC still exists in faf-heat-NA50pct and faf-heat-NA0pct(Figs. 5e and f), the magnitude of the DSL change gradient across the ACC is significantly lower when compared with that of faf-heat, which is associated with the smaller increase in OHC (Figs. 6g and h). For the North Atlantic,the increase in OHC north of 40°N in faf-heat-NA50pct is greatly reduced compared with that in faf-heat, and faf-heat-NA0pct shows almost no increase in this region, which mainly results from the suppressed heat flux perturbation imposed in the North Atlantic in both experiments. The change in the North Atlantic dipole pattern in DSL in fafheat-NA50pct is not obvious, and faf-heat-NA0pct exhibits an almost opposite dipole pattern.

    Fig.5. Spatial patterns of the DSL change (units: m) averaged in the final decade relative to the control in six experiments: (a) faf-heat, (b) faf-water, (c) faf-stress, (d) faf-all, (e) faf-heat-NA50pct, and (f) faf-heat-NA0pct.

    Using the FAFMIP multimodel output datasets, further investigations can be performed to examine future ocean circulations and energy changes with increasing COlevels.

    4.Usage notes

    The FAFMIP datasets, including those from the faf-all,faf-heat, faf-water, faf-stress, faf-passiveheat, faf-heat-NA50pct and faf-heat-NA0pct experiments, have been uploaded to the ESGF data server for users to download and can be found at https://esgf-nodes.llnl.gov/projects/cmip6/.The dataset format is the NetCDF (Network Common Data Form), version 4. The horizontal grid numbers of the model outputs are 360 and 196 in the zonal and meridional directions, respectively. The data have 30 vertical levels with 10 m per layer in the upper 150 m and the original horizontal resolutions and vertical levels are not changed on the ESGF nodes. The variables of Priority 1 for FAFMIP are shown in Table 1. The corresponding piControl data are also available from the same website.

    5.Summary

    The purpose of this paper is to introduce the FAFMIP experimental datasets in CMIP6 as simulated by CASESM2.0 and to provide a brief evaluation of the model performances in this regard. Seven FAFMIP experiments (including faf-heat, faf-stress, faf-water, faf-all, faf-passiveheat,faf-heat-NA50pct and faf-heat-NA0pct) have been conducted using CAS-ESM2.0, with each experiment integrated over 70 years. As expected, the temporal evolutions of the global mean variables (including the changes in ocean temperature, surface air temperature, surface heat flux, and AMOC) and the spatial patterns of ocean climate changes(including DSL and OHC changes) in CAS-ESM2.0 are highly consistent with the preliminary results of other models reported in Gregory et al. (2016).

    Fig.6. Spatial patterns of the OHC change (units: GJ m?2; vertical integral of the change in the tracer multiplied by the volumetric heat capacity) averaged in the final decade relative to the control in eight experiments. Experiments(a) to (f) is the same as in Fig.5, while (g) and (h) are experiments that separate the added and redistributive components in the heat flux perturbation experiments. Detailed information of these experiments can refer to Section 2.2 (Experimental design).

    As in many previous studies, surface heat flux perturbations play a dominant role in ocean climate changes (including DSL, AMOC and OHC changes) in response to COforcing among the surface flux perturbations. Heat flux and water flux perturbations cause the DSL dipole pattern in the North Atlantic, while the Southern Ocean features of DSL change are mainly caused by heat and wind stress perturbations, and water flux perturbations have a counteracting contribution. The AMOC decreases in response to surface heat flux perturbations and the wind stress and water flux perturbations have almost no effect on AMOC change. The two newly added experiments (i.e., faf-heat-NA50pct and fafheat-NA0pct) greatly reduce the weakening of AMOC, and the dipole pattern of DSL in the North Atlantic is markedly different from those in faf-heat and faf-all. Therefore, the rationality of the two new experiments needs to be further investigated. For the OHC changes, the increased OHC infaf-heat is primarily determined by the added OHC change,especially for the Atlantic Ocean and Southern Ocean, and the redistributed OHC changes due to ocean circulation changes are also important for influencing the geographical patterns of OHC changes.

    Table 1. Information on the CAS-ESM2.0 uploaded variables for FAFMIP, including name, description and frequency.

    Based on the reasonable behaviors of the model results,further research can be performed to explore the changes in OHC, ocean circulation and sea level response to COforcing by using our model outputs along with other FAFMIP datasets. The FAFMIP experimental outputs can be compared with the 1pctCO2 experiments in CMIP6 to investigate the uncertainty of ocean responses to increasing COlevels. Our experimental datasets have been uploaded to the CMIP6 website and all scientists and researchers are welcome to download and use the data. The results in this paper also serve as a basic reference of CAS-ESM2.0 performance for all users who are interested in the FAFMIP experiments and who might use our datasets.

    Acknowledgements.

    This research was supported by the National Major Research High PerformanceComputing Program of China (Grant No. 2016YFB0200804), the National Natural Science Foundation of China (Grant Nos. 41706036, 41706028,41975129 and 41630530), the open fund of State Key Laboratory of Satellite Ocean Environment Dynamics, Second Institute of Oceanography (Grant No. QNHX2017), and the National Key Scientific and Technological Infrastructure project entitled “Earth System Science Numerical Simulator Facility” (EarthLab) and key operation construction projects of Chongqing Meteorological Bureau-"Construction of chongqing short-term climate numerical prediction platform".

    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网站| 国产99久久九九免费精品| 大码成人一级视频| av又黄又爽大尺度在线免费看| 国产一区二区在线观看av| 老司机影院毛片| 国产av精品麻豆| 欧美av亚洲av综合av国产av| 亚洲专区国产一区二区| 久久人人爽av亚洲精品天堂| 亚洲人成电影免费在线| 国产精品.久久久| av线在线观看网站| 婷婷色综合www| 黑人猛操日本美女一级片| 亚洲欧美精品自产自拍| 亚洲久久久国产精品| 日本黄色日本黄色录像| 免费看不卡的av| 日本五十路高清| 国产老妇伦熟女老妇高清| 黑人欧美特级aaaaaa片| 国产精品秋霞免费鲁丝片| 国产免费福利视频在线观看| 一区福利在线观看| 国产真人三级小视频在线观看| 精品少妇一区二区三区视频日本电影| 国产亚洲午夜精品一区二区久久| 久久人妻熟女aⅴ| 国产亚洲精品第一综合不卡| xxxhd国产人妻xxx| 成人黄色视频免费在线看| 99国产综合亚洲精品| 日本91视频免费播放| 超碰成人久久| 国产一区二区 视频在线| 最新在线观看一区二区三区 | 老司机影院毛片| 亚洲欧美色中文字幕在线| 涩涩av久久男人的天堂| 男人舔女人的私密视频| 老汉色∧v一级毛片| 99热全是精品| 午夜激情av网站| 亚洲精品日本国产第一区| 最近手机中文字幕大全| 精品卡一卡二卡四卡免费| 脱女人内裤的视频| av一本久久久久| 一区二区av电影网| 色网站视频免费| avwww免费| www.自偷自拍.com| 最近手机中文字幕大全| 国产亚洲午夜精品一区二区久久| 国产精品国产三级专区第一集| 男的添女的下面高潮视频| 国产男女超爽视频在线观看| 91精品三级在线观看| 午夜福利,免费看| 1024香蕉在线观看| 欧美av亚洲av综合av国产av| 午夜日韩欧美国产| 一级片'在线观看视频| 国产麻豆69| bbb黄色大片| 欧美精品一区二区免费开放| 日韩av在线免费看完整版不卡| 精品国产乱码久久久久久男人| 看免费成人av毛片| 久久精品国产亚洲av涩爱| 两个人免费观看高清视频| 欧美在线黄色| 91成人精品电影| 晚上一个人看的免费电影| 久热爱精品视频在线9| 国产女主播在线喷水免费视频网站| 中文字幕另类日韩欧美亚洲嫩草| 欧美黑人欧美精品刺激| 欧美日本中文国产一区发布| 国产不卡av网站在线观看| 女人被躁到高潮嗷嗷叫费观| 久久久久久久大尺度免费视频| 久久精品久久久久久久性| 一区二区三区乱码不卡18| 午夜日韩欧美国产| 手机成人av网站| 好男人视频免费观看在线| 丁香六月欧美| 日日摸夜夜添夜夜爱| 中文字幕高清在线视频| 好男人视频免费观看在线| 午夜免费鲁丝| 丰满人妻熟妇乱又伦精品不卡| 爱豆传媒免费全集在线观看| 国产成人免费观看mmmm| 午夜老司机福利片| 亚洲av国产av综合av卡| 90打野战视频偷拍视频| 熟女av电影| 搡老岳熟女国产| 丝袜喷水一区| 欧美日韩黄片免| 午夜日韩欧美国产| 久久久久久久久免费视频了| 成人18禁高潮啪啪吃奶动态图| bbb黄色大片| 18禁黄网站禁片午夜丰满| 中国美女看黄片| 久久性视频一级片| 电影成人av| 永久免费av网站大全| 日本a在线网址| 亚洲九九香蕉| 18在线观看网站| 亚洲激情五月婷婷啪啪| 1024视频免费在线观看| 亚洲黑人精品在线| 观看av在线不卡| 大片电影免费在线观看免费| 免费在线观看日本一区| av天堂久久9| 国产一区亚洲一区在线观看| 可以免费在线观看a视频的电影网站| 大片电影免费在线观看免费| 国产麻豆69| 自拍欧美九色日韩亚洲蝌蚪91| 日韩中文字幕视频在线看片| 亚洲国产中文字幕在线视频| 50天的宝宝边吃奶边哭怎么回事| 校园人妻丝袜中文字幕| 一区二区三区四区激情视频| 欧美乱码精品一区二区三区| 亚洲欧美清纯卡通| 中文字幕精品免费在线观看视频| 热99国产精品久久久久久7| 久久精品国产综合久久久| 少妇猛男粗大的猛烈进出视频| 99国产精品99久久久久| 咕卡用的链子| 美女视频免费永久观看网站| 丝瓜视频免费看黄片| 女人高潮潮喷娇喘18禁视频| 国产亚洲一区二区精品| 老汉色∧v一级毛片| 秋霞在线观看毛片| 麻豆乱淫一区二区| 热re99久久精品国产66热6| 国产精品一二三区在线看| 亚洲五月色婷婷综合| 国产无遮挡羞羞视频在线观看| 午夜两性在线视频| 日韩中文字幕欧美一区二区 | 午夜影院在线不卡| 国产主播在线观看一区二区 | 人人妻人人澡人人看| av网站在线播放免费| 啦啦啦在线观看免费高清www| 多毛熟女@视频| 欧美中文综合在线视频| 丰满迷人的少妇在线观看| 王馨瑶露胸无遮挡在线观看| 久久久精品区二区三区| 精品亚洲乱码少妇综合久久| 在线精品无人区一区二区三| 精品福利永久在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 9色porny在线观看| 日韩 亚洲 欧美在线| 国产真人三级小视频在线观看| 成人免费观看视频高清| 久久国产精品大桥未久av| 狠狠婷婷综合久久久久久88av| 亚洲av欧美aⅴ国产| 欧美精品一区二区免费开放| 午夜激情久久久久久久| 好男人视频免费观看在线| 国产精品久久久久成人av| 国产视频一区二区在线看| av线在线观看网站| 久久av网站| 波野结衣二区三区在线| 久久99一区二区三区| 亚洲精品自拍成人| 伊人亚洲综合成人网| 欧美人与性动交α欧美精品济南到| www.精华液| 黄色毛片三级朝国网站| 国产片特级美女逼逼视频| 男人舔女人的私密视频| 国产亚洲精品第一综合不卡| 一区二区日韩欧美中文字幕| 国产av一区二区精品久久| 国产精品九九99| 久久精品亚洲av国产电影网| 日韩一卡2卡3卡4卡2021年| 熟女少妇亚洲综合色aaa.| 热re99久久国产66热| 国产精品免费视频内射| 亚洲精品中文字幕在线视频| 天天影视国产精品| 80岁老熟妇乱子伦牲交| 欧美日韩综合久久久久久| 色精品久久人妻99蜜桃| 免费在线观看黄色视频的| 一级片'在线观看视频| 高潮久久久久久久久久久不卡| 国产在视频线精品| 国产精品久久久久久精品古装| 欧美在线黄色| 日韩,欧美,国产一区二区三区| 男女下面插进去视频免费观看| 国产极品粉嫩免费观看在线| 精品欧美一区二区三区在线| 色视频在线一区二区三区| 成年人午夜在线观看视频| 天天躁夜夜躁狠狠躁躁| avwww免费| 亚洲人成网站在线观看播放| 女人久久www免费人成看片| 男女免费视频国产| 中国美女看黄片| 一本—道久久a久久精品蜜桃钙片| 国产97色在线日韩免费| 久久热在线av| 天天躁日日躁夜夜躁夜夜| 免费看av在线观看网站| 久久狼人影院| 19禁男女啪啪无遮挡网站| 亚洲黑人精品在线| 国产精品秋霞免费鲁丝片| 欧美日韩视频高清一区二区三区二| 天天躁日日躁夜夜躁夜夜| 欧美 日韩 精品 国产| 好男人电影高清在线观看| 丁香六月天网| 久久亚洲国产成人精品v| 你懂的网址亚洲精品在线观看| 欧美亚洲 丝袜 人妻 在线| e午夜精品久久久久久久| 两个人看的免费小视频| 午夜福利影视在线免费观看| 久久精品久久久久久久性| 国产亚洲精品第一综合不卡| 国产又色又爽无遮挡免| 黄色视频不卡| 欧美亚洲 丝袜 人妻 在线| 男女边吃奶边做爰视频| 亚洲av片天天在线观看| 精品少妇一区二区三区视频日本电影| 亚洲男人天堂网一区| 如日韩欧美国产精品一区二区三区| 多毛熟女@视频| 国产在视频线精品| 在线观看www视频免费| 国产精品香港三级国产av潘金莲 | 这个男人来自地球电影免费观看| 99精品久久久久人妻精品| 日本猛色少妇xxxxx猛交久久| 手机成人av网站| 国产亚洲精品第一综合不卡| 成人亚洲欧美一区二区av| 青草久久国产| 久久久国产一区二区| 午夜91福利影院| www.自偷自拍.com| 国产成人精品在线电影| 女人精品久久久久毛片| 一本久久精品| 亚洲第一av免费看| 国产伦人伦偷精品视频| 精品国产乱码久久久久久小说| 久久综合国产亚洲精品| 久久久久网色| 新久久久久国产一级毛片| 免费看不卡的av| 久久免费观看电影| 又粗又硬又长又爽又黄的视频| 在现免费观看毛片| 国产淫语在线视频| 交换朋友夫妻互换小说| 久久久久国产一级毛片高清牌| 91精品伊人久久大香线蕉| 日韩制服丝袜自拍偷拍| 亚洲久久久国产精品| 男人操女人黄网站| 亚洲精品国产一区二区精华液| av不卡在线播放| 亚洲伊人色综图| 女人久久www免费人成看片| 午夜老司机福利片| 国产成人a∨麻豆精品| 黑人巨大精品欧美一区二区蜜桃| 久久鲁丝午夜福利片| avwww免费| 香蕉丝袜av| 一本—道久久a久久精品蜜桃钙片| 色播在线永久视频| 丝袜喷水一区| 欧美xxⅹ黑人| 久久久精品国产亚洲av高清涩受| 91老司机精品| 韩国高清视频一区二区三区| av欧美777| 男的添女的下面高潮视频| 日本vs欧美在线观看视频| 欧美另类一区| www日本在线高清视频| 亚洲av成人不卡在线观看播放网 | 超碰成人久久| 久久久久久人人人人人| 成人亚洲精品一区在线观看| 极品人妻少妇av视频| 王馨瑶露胸无遮挡在线观看| 欧美成人精品欧美一级黄| 国产在线一区二区三区精| 欧美变态另类bdsm刘玥| 大陆偷拍与自拍| 只有这里有精品99| 国产精品一区二区在线不卡| 国产精品一区二区精品视频观看| 国产真人三级小视频在线观看| 精品久久久久久久毛片微露脸 | 久久久精品区二区三区| 国产精品久久久av美女十八| 亚洲成色77777| 亚洲美女黄色视频免费看| 国产成人啪精品午夜网站| 我要看黄色一级片免费的| videos熟女内射| 免费观看av网站的网址| 亚洲精品乱久久久久久| 欧美性长视频在线观看| 国产成人系列免费观看| 男女之事视频高清在线观看 | 脱女人内裤的视频| 91字幕亚洲| 日韩,欧美,国产一区二区三区| 国产成人影院久久av| 伦理电影免费视频| 韩国高清视频一区二区三区| 亚洲精品美女久久久久99蜜臀 | 一边摸一边抽搐一进一出视频| 一个人免费看片子| 欧美日本中文国产一区发布| 久久久精品国产亚洲av高清涩受| 欧美 亚洲 国产 日韩一| av线在线观看网站| 99国产精品免费福利视频| cao死你这个sao货| 黑人猛操日本美女一级片| 久久久久网色| 人人妻,人人澡人人爽秒播 | 在线观看免费视频网站a站| 欧美变态另类bdsm刘玥| 欧美精品一区二区大全| e午夜精品久久久久久久| av国产久精品久网站免费入址| 999久久久国产精品视频| 人人妻人人爽人人添夜夜欢视频| 久久久久久久久久久久大奶| 久久久精品94久久精品| 麻豆乱淫一区二区| 亚洲,欧美,日韩| 麻豆国产av国片精品| 亚洲色图综合在线观看| 亚洲 国产 在线| 国产精品av久久久久免费| 狂野欧美激情性xxxx| 亚洲精品日本国产第一区| 国产午夜精品一二区理论片| 亚洲精品自拍成人| av一本久久久久| 秋霞在线观看毛片| 亚洲久久久国产精品| 久久这里只有精品19| 91九色精品人成在线观看| 男女无遮挡免费网站观看| 久久人妻熟女aⅴ| 中文字幕色久视频| 亚洲国产精品成人久久小说| kizo精华| 欧美日韩福利视频一区二区| 99国产精品99久久久久| 亚洲伊人色综图| 成人亚洲精品一区在线观看| 性少妇av在线| 欧美 亚洲 国产 日韩一| 男人添女人高潮全过程视频| 性色av乱码一区二区三区2| 婷婷丁香在线五月| 少妇 在线观看| 久久亚洲国产成人精品v| 国产黄色免费在线视频| 国产人伦9x9x在线观看| 777久久人妻少妇嫩草av网站| 18禁国产床啪视频网站| 亚洲精品久久久久久婷婷小说| av国产精品久久久久影院| 女人高潮潮喷娇喘18禁视频| 久久久久久久久久久久大奶| 色视频在线一区二区三区| 伊人久久大香线蕉亚洲五| 精品国产乱码久久久久久男人| www.精华液| 老汉色∧v一级毛片| 国产成人精品在线电影| 久久久国产一区二区| 国产亚洲一区二区精品| 叶爱在线成人免费视频播放| 日本猛色少妇xxxxx猛交久久| 国产片特级美女逼逼视频| 亚洲第一av免费看| cao死你这个sao货| 久久青草综合色| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲欧美一区二区三区黑人| 校园人妻丝袜中文字幕| 亚洲国产欧美日韩在线播放| 亚洲中文av在线| 国产精品秋霞免费鲁丝片| 国产成人精品久久二区二区免费| 午夜免费观看性视频| 女警被强在线播放| 亚洲欧美中文字幕日韩二区| 日本黄色日本黄色录像| 国产成人欧美在线观看 | 色播在线永久视频| 国产成人系列免费观看| 亚洲色图 男人天堂 中文字幕| 老司机影院成人| 亚洲成色77777| 男女高潮啪啪啪动态图| 亚洲三区欧美一区| 女性生殖器流出的白浆| 精品一区在线观看国产| 日本五十路高清| 日韩人妻精品一区2区三区| 久久99热这里只频精品6学生| av欧美777| 欧美亚洲 丝袜 人妻 在线| 欧美久久黑人一区二区| av天堂在线播放| 最新的欧美精品一区二区| 男女高潮啪啪啪动态图| 亚洲精品在线美女| 国产高清不卡午夜福利| 国产高清国产精品国产三级| 天天躁夜夜躁狠狠躁躁| 免费女性裸体啪啪无遮挡网站| 精品第一国产精品| av天堂在线播放| 日韩一本色道免费dvd| 我的亚洲天堂| 亚洲国产av影院在线观看| 69精品国产乱码久久久| 久久精品亚洲熟妇少妇任你| 亚洲精品美女久久av网站| 亚洲精品国产一区二区精华液| 国产精品久久久久久精品古装| 亚洲国产精品一区二区三区在线| 真人做人爱边吃奶动态| 国产午夜精品一二区理论片| 精品久久久久久电影网| 亚洲精品第二区| 免费观看av网站的网址| 91精品国产国语对白视频| kizo精华| xxx大片免费视频| 亚洲精品国产av蜜桃| 日韩av不卡免费在线播放| 黄色毛片三级朝国网站| 久久久精品免费免费高清| 视频区图区小说| 亚洲国产精品成人久久小说| 日本黄色日本黄色录像| 18禁裸乳无遮挡动漫免费视频| 久久精品久久久久久久性| 99国产综合亚洲精品| 国产成人av激情在线播放| 97在线人人人人妻| 亚洲人成电影观看| 中文字幕精品免费在线观看视频| 老汉色av国产亚洲站长工具| 成年女人毛片免费观看观看9 | 免费不卡黄色视频| 国产亚洲av高清不卡| 一级毛片黄色毛片免费观看视频| 我的亚洲天堂| 午夜免费观看性视频| 亚洲免费av在线视频| 国产无遮挡羞羞视频在线观看| 久久精品国产亚洲av高清一级| 美女福利国产在线| 欧美激情 高清一区二区三区| 国产精品人妻久久久影院| 天堂俺去俺来也www色官网| 久久人妻熟女aⅴ| 啦啦啦在线免费观看视频4| 久久亚洲国产成人精品v| 国产精品免费大片| 女人高潮潮喷娇喘18禁视频| 最黄视频免费看| 天堂俺去俺来也www色官网| 纵有疾风起免费观看全集完整版| 美女福利国产在线| 黄片小视频在线播放| 一边亲一边摸免费视频| av在线app专区| bbb黄色大片| 韩国高清视频一区二区三区| 男女下面插进去视频免费观看| 在线观看免费高清a一片| 日韩 欧美 亚洲 中文字幕| 婷婷色综合大香蕉| 国产视频一区二区在线看| cao死你这个sao货| 久久精品国产综合久久久| 久久av网站| 少妇的丰满在线观看| 久久人人97超碰香蕉20202| 久久99一区二区三区| 热re99久久国产66热| 中文字幕人妻丝袜一区二区| 日韩中文字幕视频在线看片| 久久人人97超碰香蕉20202| 这个男人来自地球电影免费观看| 精品少妇内射三级| 五月开心婷婷网| 日本av手机在线免费观看| 欧美老熟妇乱子伦牲交| svipshipincom国产片| 精品一区二区三卡| 亚洲伊人久久精品综合| 久久久久久久国产电影| 亚洲色图 男人天堂 中文字幕| 亚洲欧美一区二区三区国产| 狠狠婷婷综合久久久久久88av| 中文字幕高清在线视频| 欧美激情极品国产一区二区三区| 天堂8中文在线网| 蜜桃在线观看..| 久久久精品94久久精品| 日韩精品免费视频一区二区三区| 亚洲成国产人片在线观看| 成年av动漫网址| 69精品国产乱码久久久| 国产成人欧美| 午夜av观看不卡| 热re99久久精品国产66热6| 电影成人av| 亚洲欧美中文字幕日韩二区| 成年美女黄网站色视频大全免费| 国语对白做爰xxxⅹ性视频网站| 亚洲中文av在线| 亚洲色图综合在线观看| 黑人欧美特级aaaaaa片| 久久精品熟女亚洲av麻豆精品| 亚洲欧洲日产国产| 男人添女人高潮全过程视频| 国产一区亚洲一区在线观看| 欧美av亚洲av综合av国产av| 99热网站在线观看| 成人国产一区最新在线观看 | 国产成人精品久久二区二区免费| 看免费成人av毛片| 我的亚洲天堂| 一个人免费看片子| 亚洲伊人久久精品综合| 99久久人妻综合| 久久久精品国产亚洲av高清涩受| 嫩草影视91久久| 99热国产这里只有精品6| 亚洲精品国产区一区二| 十八禁人妻一区二区| 欧美日韩av久久| 欧美xxⅹ黑人| 最近手机中文字幕大全| 国产av精品麻豆| 日日夜夜操网爽| 国产在线免费精品| 丁香六月欧美| 亚洲av电影在线进入| 最近手机中文字幕大全| 午夜福利在线免费观看网站| 男男h啪啪无遮挡| 亚洲精品美女久久av网站| 亚洲色图综合在线观看| 在线天堂中文资源库| 亚洲 欧美一区二区三区| 性高湖久久久久久久久免费观看| 亚洲av片天天在线观看| 大话2 男鬼变身卡| 国产老妇伦熟女老妇高清| 日韩av不卡免费在线播放| 日韩一区二区三区影片| 视频在线观看一区二区三区| 亚洲欧美精品综合一区二区三区| 亚洲精品自拍成人| 一边摸一边做爽爽视频免费| 免费高清在线观看视频在线观看| 久久久精品区二区三区| 亚洲人成电影免费在线| 十八禁高潮呻吟视频| av又黄又爽大尺度在线免费看| 在线观看免费视频网站a站| 国产人伦9x9x在线观看| 亚洲久久久国产精品| 国产1区2区3区精品| 免费高清在线观看日韩| 99国产精品一区二区三区| 国产福利在线免费观看视频|