删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

Datasets for the CMIP6 Scenario Model Intercomparison Project (ScenarioMIP) Simulations with the Cou

本站小编 Free考研考试/2022-01-02

<script type="text/javascript" src="https://cdn.bootcss.com/mathjax/2.7.2-beta.0/MathJax.js?config=TeX-AMS-MML_HTMLorMML"></script> <script type='text/x-mathjax-config'> MathJax.Hub.Config({ tex2jax: {inlineMath: [['$','$'], ['\\(','\\)']]}, "HTML-CSS": { scale: 100, minScaleAdjust: 100, mtextFontInherit: true, } }); </script> Shuwen ZHAO1,2,
Yongqiang YU1,2,3,,,
Pengfei LIN1,2,
Hailong LIU1,2,
Bian HE1,2,
Qing BAO1,2,
Yuyang GUO1,2,
Lijuan HUA4,
Kangjun CHEN1,2,
Xiaowei WANG1,2,5

Corresponding author: Yongqiang YU,yyq@lasg.iap.ac.cn;
1.State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics (LASG), Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100029, China
2.University of Chinese Academy of Sciences, Beijing 100049, China
3.Center for ocean Mega-Science, Chinese Academy of Sciences, Qingdao 266071, China
4.Chinese Academy of Meteorological Sciences, Beijing 100081, China
5.Yunnan University, Kunming 650091, China
Manuscript received: 2020-04-17
Manuscript revised: 2020-08-29
Manuscript accepted: 2020-09-08
Abstract:The datasets for the tier-1 Scenario Model Intercomparison Project (ScenarioMIP) experiments from the Chinese Academy of Sciences (CAS) Flexible Global Ocean–Atmosphere–Land System model, finite-volume version 3 (CAS FGOALS-f3-L) are described in this study. ScenarioMIP is one of the core MIP experiments in phase 6 of the Coupled Model Intercomparison Project (CMIP6). Considering future CO2, CH4, N2O and other gases’ concentrations, as well as land use, the design of ScenarioMIP involves eight pathways, including two tiers (tier-1 and tier-2) of priority. Tier-1 includes four combined Shared Socioeconomic Pathways (SSPs) with radiative forcing, i.e., SSP1-2.6, SSP2-4.5, SSP3-7.0 and SSP5-8.5, in which the globally averaged radiative forcing at the top of the atmosphere around the year 2100 is approximately 2.6, 4.5, 7.0 and 8.5 W m?2, respectively. This study provides an introduction to the ScenarioMIP datasets of this model, such as their storage location, sizes, variables, etc. Preliminary analysis indicates that surface air temperatures will increase by about 1.89°C, 3.07°C, 4.06°C and 5.17°C by around 2100 under these four scenarios, respectively. Meanwhile, some other key climate variables, such as sea-ice extension, precipitation, heat content, and sea level rise, also show significant long-term trends associated with the radiative forcing increases. These datasets will help us understand how the climate will change under different anthropogenic and radiative forcings.
Keywords: CMIP6,
ScenarioMIP,
FGOALS-f3-L,
coupled model,
data description
摘要:本文介绍了中国科学院大气物理研究所大气科学和地球流体力学数值模拟国家重点实验室研发的耦合气候系统模式CAS FGOALS-f3-L参与第六阶段耦合模式比较计划(CMIP6)情景模式比较计划(SceanrioMIP)第一级(tier-1)的实验数据集。SceanrioMIP是CMIP6的核心比较计划之一。考虑到未来的CO2、CH4、N2O和其他气体浓度,以及土地利用,SceanrioMIP的设计涉及8条路径,包括两层优先级(tier-1和tier-2)。Tier-1包括四个具有不同辐射强迫的共享社会经济路径(SSPs),即SSP1-2.6、SSP2-4.5、SSP3-7.0和SSP5-8.5,其中2100年前后全球大气顶部平均辐射强迫分别约为2.6、4.5、7.0和8.5 W m?2。本研究对CAS FGOALS-f3-L的ScenarioMIP数据集进行了介绍,如存储位置、大小、变量等。初步分析表明,这四种情景下,到2100年前后地表气温将分别升高约1.89℃、3.07℃、4.06℃和5.17℃。与此同时,其他一些关键的气候变量,如海冰面积、降水、海水热含量和海平面上升高度,也显示出与辐射强迫增加有关的显著的长期趋势。这些数据集将帮助我们了解在不同的人为和辐射强迫作用下气候将如何变化。
关键词:CMIP6,
ScenarioMIP,
FGOALS-f3-L,
耦合气候系统模式,
数据描述





--> --> -->
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
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.2. Experimental design
--> 3
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 range of 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.
Scenario nameForcing category2100 forcing (W m?2)SSP
SSP1-2.6Low2.61
SSP2-4.5Medium4.52
SSP3-7.0High7.03
SSP5-8.5High8.55


Table1. ScenarioMIP experiment designs.



3
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 ScenarioMIP-provided concentrations for all long-lived GHGs, such as CO2, CH4, N2O and CFCs. The concentrations of CO2, CH4 and N2O in our model are shown in Fig. 1, in which it can be seen that, for concentrations of CH4 and N2O, 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). CO2 still 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-to-date 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.
Figure1. 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.



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.

2
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).
Figure2. 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.


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 m2 around 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 m?2 above the aggregate mean of CMIP5 (Collins et al., 2013).
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 106 km2 (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 106 km2 for SSP1-2.6 after 2055. The September sea-ice extent in the Arctic of the CMIP5 multi-model 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 106 km2 around 2060 and continues to decline to 0 km2. 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)?1, respectively. Moreover, comparing surface air temperature, heat content and sea level rise in SSP1-2.6, we find that, after 2050, the concentrations of CO2 stop 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% °C?1 under global warming, which is roughly similar to previous estimates by Held and Soden (2006).

2
3.2. Changes in surface air temperature
--> 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.
Figure3. 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).



2
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).
Figure4. 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).



2
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).
Figure5. 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).



2
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.
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 CO2 experiments, 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/most-recent/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 atmospheric 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.
Variable abbreviationLong nameunits
cltTotal cloud cover percentage%
evspsblsoiWater evaporation from soilkg m?2 s?1
evspsblvegEvaporation from canopykg m?2 s?1
hfbasinNorthward ocean heat transportW
hfdsDownward heat flux at seawater surfaceW m?2
hurRelative humidity%
hursNear-surface relative humidity%
husSpecific humidity1
hussNear-surface specific humidity1
mlotstOcean mixed layer thickness defined by sigma Tm
mrfsoSoil frozen water contentkg m?2
mrroTotal runoffkg m?2 s?1
mrrosSurface runoffkg m?2 s?1
mrsoTotal soil moisture contentkg m?2
mrsosMoisture in upper portion of soil columnkg m?2
msftbarotOcean barotropic mass streamfunctionkg s?1
msftmzOcean meridional overturning mass streamfunctionkg s?1
msftmzmpaOcean meridional overturning mass streamfunction due to parameterized mesoscale advectionkg s?1
prPrecipitationkg m?2 s?1
prsnSnowfall fluxkg m?2 s?1
prvegPrecipitation onto canopykg m?2 s?1
pslSea level pressurePa
psSurface air pressurePa
sfcWindNear-surface wind speedm s?1
siconcSea-ice area percentage (ocean grid)%
siconcaSea-ice area percentage (atmospheric grid)%
siprRainfall rate over sea icekg m?2 s?1
sisnconcSnow area percentage%
sisnthickSnow thicknessm
siuX-component of sea-ice velocitym s?1
sivY-component of sea-ice velocitym s?1
sosSea surface salinity0.001
rlutTOA outgoing longwave radiationW m?2
rsdtTOA incident shortwave radiationW m?2
rsutTOA outgoing shortwave radiationW m?2
taAir temperatureK
tasNear-surface air temperatureK
tauuSurface downward eastward wind stressPa
tauvSurface downward northward wind stressPa
tsSurface temperatureK
tslTemperature of soilK
thetaoSeawater potential temperature°C
thetaogaGlobal average seawater potential temperature°C
tobSeawater potential temperature at sea floor°C
tosSea surface temperature°C
tosgaGlobal average sea surface temperature°C
uasEastward near-surface windm s?1
uoSeawater X velocitym s?1
vasNorthward near-surface windm s?1
voSeawater Y velocitym s?1
wapOmega (= dp / dt)Pa s?1
wfoWater flux into seawaterkg m?2 s?1
zgGeopotential heightm
zosSea surface height above geoidm


Table2. 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).


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.

相关话题/Datasets CMIP6 Scenario