Advanced Search
Article Contents

Surface Temperature Changes Projected by FGOALS Models under Low Warming Scenarios in CMIP5 and CMIP6

Fund Project:

This work was supported by the Natural Science Foundation of China (Grant Nos. 41706026 and 41831175), the National Key Research and Development Program of China (Grant Nos. 2017YFA0604600, 2016YFA0601804 and 2018YFA0605702), the Fundamental Research Funds for the Central Universities (Grant Nos. 2018B04814 and 2018B03114), the open fund of the State Key Laboratory of Satellite Ocean Environment Dynamics, Second Institute of Oceanography (Grant No. QNHX1808), and the National Basic Research Program of China (Grant No. 2012CB955600)


doi: 10.1007/s00376-020-0177-5

  • To meet the low warming targets proposed in the 2015 Paris Agreement, substantial reduction in carbon emissions is needed in the future. It is important to know how surface climates respond under low warming targets. The present study investigates the surface temperature changes under the low-forcing scenario of Representative Concentration Pathways (RCP2.6) and its updated version (Shared Socioeconomic Pathways, SSP1-2.6) by the Flexible Global Ocean–Atmosphere–Land System (FGOALS) models participating in phases 5 and 6 of the Coupled Model Intercomparison Project (CMIP5 and CMIP6, respectively). In both scenarios, radiative forcing (RF) first increases to a peak of 3 W m−2 around 2045 and then decreases to 2.6 W m−2 by 2100. Global mean surface air temperature rises in all FGOALS models when RF increases (RF increasing stage) and declines or holds nearly constant when RF decreases (RF decreasing stage). The surface temperature change is distinct in its sign and magnitude between the RF increasing and decreasing stages over the land, Arctic, North Atlantic subpolar region, and Southern Ocean. Besides, the regional surface temperature change pattern displays pronounced model-to-model spread during both the RF increasing and decreasing stages, mainly due to large intermodel differences in climatological surface temperature, ice-albedo feedback, natural variability, and Atlantic Meridional Overturning Circulation change. The pattern of tropical precipitation change is generally anchored by the spatial variations of relative surface temperature change (deviations from the tropical mean value) in the FGOALS models. Moreover, the projected changes in the updated FGOALS models are closer to the multi-model ensemble mean results than their predecessors, suggesting that there are noticeable improvements in the future projections of FGOALS models from CMIP5 to CMIP6.
    摘要: 为了实现2015年巴黎协定提出的低增温目标,人类社会需要在未来大量地减少碳排放量,相应的未来情景称之为低增温情景。全球表面温度在低增温背景下如何演变是气候变化中的重要问题。本文研究了在CMIP5中低强迫情景(RCP2.6)和CMIP6中相应升级的低强迫情景(SSP1-2.6)下,4个FGOALS模式预估的未来全球表面温度的变化趋势及相应机制。在RCP2.6和SSP1-2.6两种情景下,外辐射强迫(RF)都是先不断增长,在2045年左右达到3W m-2的峰值,然后在2100年下降到2.6W m-2。在4个FGOALS模式中,全球平均表面气温都是在RF增长的时候(RF上升阶段)不断上升,而在RF下降的时候(RF下降阶段)保持稳定不变或者下降。在RF上升和下降两个不同阶段,表面温度的变化符号和变化强度在陆地、北冰洋、北大西洋副极地地区和南极都有非常显著的差异。此外,表面温度变化的区域分布特征在这两个阶段则都表现出显著的模式间差异性,这主要是因为不同模式在模拟气候态表面温度、冰-反照率反馈机制、自然变率和北大西洋经向翻转环流的变化方面存在非常大的区别。在FGOALS模式中,热带降水的空间分布特征大都受表面温度相对变化(即与热带平均增温强度的偏差)的空间结构所调控。FGOALS模式在低增温情景下得到的预估结果,在CMIP6中相对于在CMIP5中更接近于多模式集合平均的结果。这表明从CMIP5到CMIP6,FGOALS模式对低增温情景下的未来预估模拟有显著提高。
  • 加载中
  • Figure 1.  Spatial distribution of annual mean biases in SST (units: °C) and 2 m air temperature (TAS; units: °C) in (a, e) FGOALS-g2, (b, f) FGOALS-g3, (c, g) FGOALS-s2, and (d, h) FGOALS-f3-L, during 1979–2005. The ERSST.v5 and CRU TS4.04 data are referenced as the observations. The maximum value (Max), minimum value (Min), RMSE, and mean value (Mean) of global biases are labeled at the top of each panel.

    Figure 2.  Global mean annual surface air temperature change (units: °C) relative to pre-industrial level (1850–99 mean) in each FGOALS model (colored solid lines) and the CMIP5 and CMIP6 MMEs (black lines), with an 11-year running mean applied. The vertical dashed lines indicate the separation point for trend calculation. Note that the gray lines show the pathway of RF [right-hand axis in (b, d)], and the light green shadows indicate the range of the CMIP models.

    Figure 3.  Global mean annual surface air temperature change (units: °C) over ocean (green lines) and land (red lines) and their difference (land–sea warming contrast, brown lines) relative to the pre-industrial level (1850–99 mean), with an 11-year running mean applied.

    Figure 4.  Linear trend of annual surface skin temperature (units: °C), which is equivalent to SST over ice-free ocean, during (a–d) 1850–2050 and (e–h) 2050–2100 in the four FGOALS models. The white (black) dots indicate the trend is insignificant (significant) at the 95% level. The global mean value is labeled in the title of each plot.

    Figure 5.  MME linear trends of annual surface temperature (units: °C) during (a, b) 1850–2050 and (c, d) 2050–2100 in (a, c) CMIP5 and (b, d) CMIP6, along with (e–h) their SNRs, defined as the absolute value of MME change divided by the intermodel standard deviation. The magenta and black contours indicate SNRs of 1 and 3, respectively.

    Figure 6.  Relative surface temperature change (∆TS*, contours, interval = 0.15°C) and precipitation change (shading) in the tropical oceans during (a–d) 1850–2050 and (e–h) 2050–2100.

    Figure 7.  Area-weighted mean annual surface temperature (units: °C) in the (a, b) tropical oceans (20°S–20°N), (c, d) Arctic (60°–90°N), and (e, f) Southern Ocean (65°–50°S), in FGOALS-g2 (blue lines), FGOALS-g3 (red lines), FGOALS-s2 (green lines), and FGOALS-f3-L (magenta lines), with an 11-year running mean applied. The colored dashed lines are the least-squares linear fitting lines for 1850–2050 and 2050–2100.

    Figure 8.  Linear trend of annual 850 hPa zonal wind (units: m s−1) during (a–d) 1850–2050 and (e–h) 2050–2100 in the four FGOALS models. The white (black) dots indicate the trend is insignificant (significant) at the 95% level.

    Figure 9.  Area-weighted mean annual surface temperature change (unit: °C) in the NA warming hole region (WH index), with an 11-year running mean applied. The regions are marked as blue boxes in Fig. 3, which are (45°–60°N, 65°–30°W) in FGOALS-g2 and FGOALS-g3 and (40°–60°N, 40°–10°W) in FGOALS-s2 and FGOALS-f3-L. The brown lines indicate the AMO index, defined as the detrended area-weighted mean SST over (0°–60°N, 80°W–0°). The change is calculated as the difference between pre-industrial control simulations and simulations from historical and future scenarios. The colored dashed lines are the least-squares linear fitting lines for 1850–2050 and 2050–2100.

    Figure 10.  Area-weighted mean annual surface temperature (unit: °C) in the (a, b) NA warming hole region and (c, d) AMOC index (unit: Sv), with an 11-year running mean applied. Note that the curves display the annual mean AMOC index in each year, but not the AMOC change. The colored dashed lines are the least-squares linear fitting lines for 1850–2050 and 2050–2100.

    Table 1.  Model components and corresponding horizontal resolutions of FGOALS models in CMIP5 and CMIP6. The components of all the models are the Finite-volume Atmospheric model (FAMIL), the Spectral Atmospheric Model of IAP/LASG (SAMIL), the LASG/IAP Climate system Ocean Model (LICOM), the Community Land Model (CLM), the Community Sea Ice Model (CSIM), and the Los Alamos sea ice model (CICE). The version number of the component models are labeled after the acronyms.

    CMIPModelComponent
    OceanSea iceAtmosphereLand
    CMIP5 FGOALS-g2 LICOM2 CICE4 GAMIL2 CLM3
    360 × 196 128 × 60
    30 levels 26 levels
    FGOALS-s2 LICOM2 CSIM5 SAMIL2 CLM3
    360 × 196 128 × 108
    30 levels 26 levels
    CMIP6 FGOALS-g3 LICOM3 CICE4 GAMIL3 CLM4
    360 × 218 188 × 80
    30 levels 26 levels
    FGOALS-f3-L LICOM3 CICE4 FAMIL CLM4
    360 × 218 288 × 180
    30 levels 32 levels
    DownLoad: CSV

    Table 2.  CMIP models used in the present study.

    CMIP5CMIP6
    1BCC_CSM1.1(m)BCC-CSM2-MR
    2CanESM2CanESM5
    3CESM1-CAM5CESM2
    4CNRM-CM5CNRM-CM6-1
    5GFDL-CM3GFDL-ESM4
    6GISS-E2-RGISS-E2-1-G
    7HadGEM2-ESHadGEM3-GC31-LL
    8IPSL-CM5A-LRIPSL-CM6A-LR
    9MIROC5MIROC6
    10MIROC-ESMMIROC-ES2L
    11MPI-ESM-LRMPI-ESM1-2-LR
    12MPI-ESM-MRMPI-ESM1-2-HR
    13MRI-CGCM3MRI-ESM2-0
    14NorESM1-MNorESM2-LM
    15NorESM1-MENorESM2-MM
    DownLoad: CSV

    Table 3.  Pattern correlations between FGOALS models and their corresponding MME results. The bold values indicate that the correlations are significant at 95% confidence level.

    CMIP MME and modelsGlobalTropical oceans
    1850–20502050–21001850–20502050–2100
    CMIP5 MMEFGOALS-g20.86−0.110.73−0.21
    FGOALS-s20.930.380.370.48
    CMIP6 MMEFGOALS-g30.790.500.80−0.04
    FGOALS-f3-L0.93−0.020.59−0.10
    DownLoad: CSV
  • Allen, M. R., and W. J. Ingram, 2002: Constraints on future changes in climate and the hydrologic cycle. Nature, 419, 224−232, https://doi.org/10.1038/nature01092.
    Andrews, T., P. M. Forster, O. Boucher, N. Bellouin, and A. Jones, 2010: Precipitation, radiative forcing and global temperature change. Geophys. Res. Lett., 37, L14701, https://doi.org/10.1029/2010GL043991.
    Armour, K. C., J. Marshall, J. R. Scott, A. Donohoe, and E. R. Newsom, 2016: Southern ocean warming delayed by circumpolar upwelling and equatorward transport. Nature Geoscience, 9, 549−554, https://doi.org/10.1038/ngeo2731.
    Bao, Q., and Coauthors, 2013: The flexible global ocean-atmosphere-land system model, spectral version 2: FGOALS-s2. Adv. Atmos. Sci., 30, 561−576, https://doi.org/10.1007/s00376-012-2113-9.
    Bindoff, N. L., and Coauthors, 2019: Changing ocean, marine ecosystems, and dependent communities. IPCC Special Report on the Ocean and Cryosphere in a Changing Climate, H.-O. Pörtner et al., Eds., In press.
    Chadwick, R., I. Boutle, and G. Martin, 2013: Spatial patterns of precipitation change in CMIP5: Why the rich do not get richer in the tropics. J. Climate, 26, 3803−3822, https://doi.org/10.1175/JCLI-D-12-00543.1.
    Chen, C. L., G. H. Wang, S. P. Xie, and W. Liu, 2019a: Why does global warming weaken the gulf stream but intensify the kuroshio? J. Climate, 32, 7437−7451, https://doi.org/10.1175/JCLI-D-18-0895.1.
    Chen, L., X. Qu, G. Huang, and Y. F. Gong, 2019b: Projections of East Asian summer monsoon under 1.5 C and 2 C warming goals. Theor. Appl. Climatol., 137, 2187−2201, https://doi.org/10.1007/s00704-018-2720-1.
    Cheng, L. J., J. Abraham, Z. Hausfather, and K. E. Trenberth, 2019: How fast are the oceans warming? Science, 363, 128−129, https://doi.org/10.1126/science.aav7619.
    Chou, C., and J. D. Neelin, 2004: Mechanisms of global warming impacts on regional tropical precipitation. J. Climate, 17, 2688−2701, https://doi.org/10.1175/1520-0442(2004)017<2688:MOGWIO>2.0.CO;2.
    Chou, C., J. D. Neelin, C. A. Chen, and J. Y. Tu, 2009: Evaluating the “rich-get-richer” mechanism in tropical precipitation change under global warming. J. Climate, 22, 1982−2005, https://doi.org/10.1175/2008JCLI2471.1.
    Church, J. A., and N. J. White, 2011: Sea-level rise from the late 19th to the early 21st century. Surveys in Geophysics, 32, 585−602, https://doi.org/10.1007/s10712-011-9119-1.
    Chylek, P., C. K. Folland, G. Lesins, M. K. Dubey, and M. Y. Wang, 2009: Arctic air temperature change amplification and the Atlantic Multidecadal Oscillation. Geophys. Res. Lett., 36, L14801, https://doi.org/10.1029/2009GL038777.
    Cohen, J., and Coauthors, 2014: Recent Arctic amplification and extreme mid-latitude weather. Nature Geoscience, 7, 627−637, https://doi.org/10.1038/ngeo2234.
    Drijfhout, S., G. J. van Oldenborgh, and A. Cimatoribus, 2012: Is a decline of AMOC causing the warming hole above the North Atlantic in observed and modeled warming patterns? J. Climate, 25, 8373−8379, https://doi.org/10.1175/JCLI-D-12-00490.1.
    Eyring, V., S. Bony, G. A. Meehl, C. A. Senior, B. Stevens, R. J. Stouffer, and K. E. Taylor, 2016: Overview of the Coupled Model Intercomparison Project Phase 6 (CMIP6) experimental design and organization. Geoscientific Model Development, 9, 1937−1958, https://doi.org/10.5194/gmd-9-1937-2016.
    Francis, J. A., and E. Hunter, 2007: Changes in the fabric of the Arctic’s greenhouse blanket. Environmental Research Letters, 2, 045011, https://doi.org/10.1088/1748-9326/2/4/045011.
    Graversen, R. G., T. Mauritsen, M. Tjernström, E. Källén, and G. Svensson, 2008: Vertical structure of recent Arctic warming. Nature, 451, 53−56, https://doi.org/10.1038/nature06502.
    Guo, Y. Y., Y. Q. Yu, P. F. Lin, H. L. Liu, B. He, Q. Bao, S. W. Zhao, and X. W. Wang, 2020: Overview of the CMIP6 historical experiment datasets with the climate system model CAS FGOALS-f3-L. Adv. Atmos. Sci., 1057−1066, https://doi.org/10.1007/s00376-020-2004-4.
    Harris, I., T. J. Osborn, P. Jones, and D. Lister, 2020: Version 4 of the CRU TS monthly high-resolution gridded multivariate climate dataset. Scientific Data, 7, 109, https://doi.org/10.1038/s41597-020-0453-3.
    He, B., and Coauthors, 2019: CAS FGOALS-f3-L model datasets for CMIP6 historical atmospheric model intercomparison project simulation. Adv. Atmos. Sci., 36, 771−778, https://doi.org/10.1007/s00376-019-9027-8.
    He, B., and Coauthors, 2020: CAS FGOALS-f3-L model datasets for CMIP6 GMMIP tier-1 and tier-3 experiments. Adv. Atmos. Sci., 37, 18−28, https://doi.org/10.1007/s00376-019-9085-y.
    Held, I. M., and B. J. Soden, 2006: Robust responses of the hydrological cycle to global warming. J. Climate, 19, 5686−5699, https://doi.org/10.1175/JCLI3990.1.
    Held, I. M., M. Winton, K. Takahashi, T. Delworth, F. R. Zeng, and G. K. Vallis, 2010: Probing the fast and slow components of global warming by returning abruptly to preindustrial forcing. J. Climate, 23, 2418−2427, https://doi.org/10.1175/2009JCLI3466.1.
    Holland, M. M., and C. M. Bitz, 2003: Polar amplification of climate change in coupled models. Climate Dyn., 21, 221−232, https://doi.org/10.1007/s00382-003-0332-6.
    Huang, B. Y., and Coauthors, 2017: Extended reconstructed Sea surface temperature, Version 5 (ERSSTv5): Upgrades, validations, and intercomparisons. J. Climate, 30, 8179−8205, https://doi.org/10.1175/JCLI-D-16-0836.1.
    Huang, P., S. P. Xie, K. M. Hu, G. Huang, and R. Huang, 2013: Patterns of the seasonal response of tropical rainfall to global warming. Nature Geoscience, 6, 357−361, https://doi.org/10.1038/ngeo1792.
    IPCC, 2013: Climate Change 2013: The Physical Science Basis. Cambridge University Press, 1535 pp., https://doi.org/10.1017/CBO9781107415324.
    IPCC, 2018: Global Warming of 1.5°C. An IPCC Special Report on the impacts of global warming of 1.5°C above pre-industrial levels and related global greenhouse gas emission pathways, in the context of strengthening the global response to the threat of climate change, sustainable development, and efforts to eradicate poverty. V. Masson-Delmotte,, et al., Eds., In Press.
    Johnson, N. C., and S. P. Xie, 2010: Changes in the sea surface temperature threshold for tropical convection. Nature Geoscience, 3, 842−845, https://doi.org/10.1038/ngeo1008.
    Kamae, Y., and M. Watanabe, 2013: Tropospheric adjustment to increasing CO2: Its timescale and the role of land-sea contrast. Climate Dyn., 41, 3007−3024, https://doi.org/10.1007/s00382-012-1555-1.
    Keil, P., T. Mauritsen, J. Jungclaus, C. Hedemann, D. Olonscheck, and R. Ghosh, 2020: Multiple drivers of the North Atlantic warming hole. Nature Climate Change, 10, 667−671, https://doi.org/10.1038/s41558-020-0819-8.
    Kim, H., and S. I. An, 2013: On the subarctic North Atlantic cooling due to global warming. Theor. Appl. Climatol., 114, 9−19, https://doi.org/10.1007/s00704-012-0805-9.
    Kumar, A., and Coauthors, 2010: Contribution of sea ice loss to Arctic amplification. Geophys. Res. Lett., 37, L21701, https://doi.org/10.1029/2010GL045022.
    Lambert, F. H., and M. J. Webb, 2008: Dependency of global mean precipitation on surface temperature. Geophys. Res. Lett., 35, L16706, https://doi.org/10.1029/2008GL034838.
    Li, D. H., T. J. Zhou, L. W. Zou, W. X. Zhang, and L. X. Zhang, 2018: Extreme high-temperature events over east Asia in 1.5°C and 2°C warmer futures: Analysis of NCAR CESM low-warming experiments. Geophys. Res. Lett., 45, 1541−1550, https://doi.org/10.1002/2017GL076753.
    Li, L. J., and Coauthors, 2013: The flexible global ocean-atmosphere-land system model, grid-point version 2: FGOALS-g2. Adv. Atmos. Sci., 30, 543−560, https://doi.org/10.1007/s00376-012-2140-6.
    Li, L. J., and Coauthors, 2020: The flexible global ocean-atmosphere-land system model grid-point version 3 (FGOALS-g3): Description and evaluation. Journal of Advances in Modeling Earth Systems, https://doi.org/10.1029/2019ms002012.
    Lin, P. F., H. L. Liu, Y. Q. Yu, and T. J. Zhou, 2013: Long-term behaviors of two versions of FGOALS2 in preindustrial control simulations with implications for 20th century simulations. Adv. Atmos. Sci., 30, 577−592, https://doi.org/10.1007/s00376-013-2186-0.
    Long, S. M., S. P. Xie, X. T. Zheng, and Q. Y. Liu, 2014: Fast and slow responses to global warming: Sea surface temperature and precipitation patterns. J. Climate, 27, 285−299, https://doi.org/10.1175/JCLI-D-13-00297.1.
    Long, S. M., S. P. Xie, and W. Liu, 2016: Uncertainty in tropical rainfall projections: Atmospheric circulation effect and the Ocean Coupling. J. Climate, 29, 2671−2687, https://doi.org/10.1175/jcli-d-15-0601.1.
    Long, S. M., S. P. Xie, Q. Y. Liu, X. T. Zheng, G. Huang, K. M. Hu, and Y. Du, 2018: Slow ocean response and the 1.5 and 2℃ warming targets. Chinese Science Bulletin, 63, 558−570, https://doi.org/10.1360/N972017-01115. (in Chinese with English abstract)
    Long, S.-M., S.-P. Xie, Y. Du, Q. Y. Liu, X.-T. Zheng, G. Huang, K.-M. Hu, and J. Ying, 2020: Effects of ocean slow response under low warming targets. J. Climate, 33, 477−496, https://doi.org/10.1175/jcli-d-19-0213.1.
    Ma, J., S. P. Xie, and Y. Kosaka, 2012: Mechanisms for tropical tropospheric circulation change in response to global warming. J. Climate, 25, 2979−2994, https://doi.org/10.1175/JCLI-D-11-00048.1.
    Mann, M. E., 2009: Defining dangerous anthropogenic interference. Proceedings of the National Academy of Sciences of the United States of America, 106, 4065−4066, https://doi.org/10.1073/pnas.0901303106.
    Menary, M. B., and R. A. Wood, 2018: An anatomy of the projected North Atlantic warming hole in CMIP5 models. Climate Dyn., 50, 3063−3080, https://doi.org/10.1007/s00382-017-3793-8.
    Mitchell, J. F. B., C. A. Wilson, and W. M. Cunnington, 1987: On CO2 climate sensitivity and model dependence of results. Quart. J. Roy. Meteor. Soc., 113, 293−322, https://doi.org/10.1002/qj.49711347517.
    Nangombe, S., T. J. Zhou, W. X. Zhang, B. Wu, S. Hu, L. W. Zou, and D. H. Li, 2018: Record-breaking climate extremes in Africa under stabilized 1.5°C and 2°C global warming scenarios. Nature Climate Change, 8, 375−380, https://doi.org/10.1038/s41558-018-0145-6.
    Nie, Y., L. J. Li, Y. L. Tang, and B. Wang, 2019: Impacts of changes of external forcings from CMIP5 to CMIP6 on surface temperature in FGOALS-g2. SOLA, 15, 211−215, https://doi.org/10.2151/SOLA.2019-038.
    Palter, J. B., T. L. Frölicher, D. Paynter, and J. G. John, 2018: Climate, ocean circulation, and sea level changes under stabilization and overshoot pathways to 1.5 K warming. Earth System Dynamics, 9, 817−828, https://doi.org/10.5194/esd-9-817-2018.
    Pithan, F., and T. Mauritsen, 2014: Arctic amplification dominated by temperature feedbacks in contemporary climate models. Nature Geoscience, 7, 181−184, https://doi.org/10.1038/ngeo2071.
    Purkey, S. G., and G. C. Johnson, 2010: Warming of global abyssal and deep Southern Ocean waters between the 1990s and 2000s: Contributions to global heat and sea level rise budgets. J. Climate, 23, 6336−6351, https://doi.org/10.1175/2010JCLI3682.1.
    Qu, X., and G. Huang, 2018: Different multi-year mean temperature in mid-summer of South China under different 1.5°C warming scenarios. Scientific Reports, 8, 13794, https://doi.org/10.1038/s41598-018-32277-6.
    Sanderson, B. M., B. C. O’Neill, and C. Tebaldi, 2016: What would it take to achieve the Paris temperature targets? Geophys Res. Lett., 43, 7133−7142, https://doi.org/10.1002/2016GL069563.
    Schaeffer, M., W. Hare, S. Rahmstorf, and M. Vermeer, 2012: Long-term sea-level rise implied by 1.5°C And 2°C warming levels. Nature Climate Change, 2, 867−870, https://doi.org/10.1038/nclimate1584.
    Schleussner, C. F., P. Pfleiderer, and E. M. Fischer, 2017: In the observational record half a degree matters. Nature Climate Change, 7, 460−462, https://doi.org/10.1038/nclimate3320.
    Schweiger, A. J., R. W. Lindsay, S. Vavrus, and J. A. Francis, 2008: Relationships between Arctic sea ice and clouds during autumn. J. Climate, 21, 4799−4810, https://doi.org/10.1175/2008JCLI2156.1.
    Screen, J. A., and I. Simmonds, 2010: The central role of diminishing sea ice in recent Arctic temperature amplification. Nature, 464, 1334−1337, https://doi.org/10.1038/nature09051.
    Seager, R., N. Naik, and G. A. Vecchi, 2010: Thermodynamic and dynamic mechanisms for large-scale changes in the hydrological cycle in response to global warming. J. Climate, 23, 4651−4668, https://doi.org/10.1175/2010JCLI3655.1.
    Serreze, M. C., A. P. Barrett, J. C. Stroeve, D. N. Kindig, and M. M. Holland, 2009: The emergence of surface-based Arctic amplification. Cryosphere, 3, 11−19, https://doi.org/10.5194/tc-3-11-2009.
    Serreze, M. C., and R. G. Barry, 2011: Processes and impacts of Arctic amplification: A research synthesis. Global and Planetary Change, 77, 85−96, https://doi.org/10.1016/j.gloplacha.2011.03.004.
    Sgubin, G., D. Swingedouw, S. Drijfhout, Y. Mary, and A. Bennabi, 2017: Abrupt cooling over the North Atlantic in modern climate models. Nature Communications, 8, 14375, https://doi.org/10.1038/ncomms14375.
    Simmonds, I., and K. Keay, 2009: Extraordinary september arctic sea ice reductions and their relationships with storm behavior over 1979-2008. Geophys. Res. Lett., 36, L19715, https://doi.org/10.1029/2009GL039810.
    Song, M. R., L. J. Li, and J. P. Liu, 2014: Sea ice simulations of FGOALS. Flexible Global Ocean-Atmosphere-Land System Model: A Modeling Tool for the Climate Change Research Community, T. J. Zhou et al., Eds., Springer, 303−309, https://doi.org/10.1007/978-3-642-41801-3_36.
    Taylor, K. E., R. J. Stouffer, and G. A. Meehl, 2012: An overview of CMIP5 and the experiment design. Bull. Amer. Meteor. Soc., 93, 485−498, https://doi.org/10.1175/BAMS-D-11-00094.1.
    van Vuuren, D. P., and Coauthors, 2011: The representative concentration pathways: An overview. Climatic Change, 109, 5−31, https://doi.org/10.1007/s10584-011-0148-z.
    Vecchi, G. A., and B. J. Soden, 2007: Global warming and the weakening of the tropical circulation. J. Climate, 20, 4316−4340, https://doi.org/10.1175/JCLI4258.1.
    Vecchi, G. A., B. J. Soden, A. T. Wittenberg, I. M. Held, A. Leetmaa, and M. J. Harrison, 2006: Weakening of tropical Pacific atmospheric circulation due to anthropogenic forcing. Nature, 441, 73−76, https://doi.org/10.1038/nature04744.
    Wang, G. H., S. P. Xie, R. X. Huang, and C. J. Chen, 2015: Robust warming pattern of global subtropical oceans and its mechanism. J. Climate, 28, 8574−8584, https://doi.org/10.1175/JCLI-D-14-00809.1.
    Winton, M., 2006: Amplified Arctic climate change: What does surface albedo feedback have to do with it? Geophys. Res. Lett., 33, L03701, https://doi.org/10.1029/2005GL025244.
    Wu, P. L., R. Wood, J. Ridley, and J. Lowe, 2010: Temporary acceleration of the hydrological cycle in response to a CO2 rampdown. Geophys. Res. Lett., 37, L12705, https://doi.org/10.1029/2010GL043730.
    Xie, S. P., C. Deser, G. A. Vecchi, J. Ma, H. Y. Teng, and A. T. Wittenberg, 2010: Global warming pattern formation: Sea surface temperature and rainfall. J. Climate, 23, 966−986, https://doi.org/10.1175/2009JCLI3329.1.
    Xu, Y. Y., and V. Ramanathan, 2017: Well below 2°C: Mitigation strategies for avoiding dangerous to catastrophic climate changes. Proceedings of the National Academy of Sciences of the United States of America, 114, 10 315−10 323, https://doi.org/10.1073/pnas.1618481114.
    Zelinka, M. D., T. A. Myers, D. T. McCoy, S. Po-Chedley, P. M. Caldwell, P. Ceppi, S. A. Klein, and K. E. Taylor, 2020: Causes of higher climate sensitivity in CMIP6 models. Geophys. Res. Lett., 47, e2019GL085782, https://doi.org/10.1029/2019GL085782.
    Zhang, W. X., T. J. Zhou, L. W. Zou, L. X. Zhang, and X. L. Chen, 2018: Reduced exposure to extreme precipitation from 0.5°C less warming in global land monsoon regions. Nature Communications, 9, 3153, https://doi.org/10.1038/s41467-018-05633-3.
    Zhou, T. J., and F. F. Song, 2014: Representative concentration pathway (RCP) projection of climate change by FGOALS. Flexible Global Ocean-Atmosphere-Land System Model: A Modeling Tool for the Climate Change Research Community, T. J. Zhou et al., Eds., Springer, 267−274.
    Zhou, T. J., F. F. Song, and X. L. Chen, 2013: Historical evolution of global and regional surface air temperature simulated by FGOALS-s2 and FGOALS-g2: How reliable are the model results? Adv. Atmos. Sci., 30, 638−657, https://doi.org/10.1007/s00376-013-2205-1.
    Zhou, T. J., and Coauthors, 2020: Development of climate and earth system models in China: Past achievements and new CMIP6 results. Journal of Meteorological Research, 34, 1−19, https://doi.org/10.1007/s13351-020-9164-0.
  • [1] Haijin DAI, Qiang YAO, 2023: Role of Ocean Dynamics in the Seasonal Hadley Cell: A Response to Idealized Arctic Amplification, ADVANCES IN ATMOSPHERIC SCIENCES, 40, 2211-2223.  doi: 10.1007/s00376-022-2057-7
    [2] Yanchi LIU, Yunqi KONG, Qinghua YANG, Xiaoming HU, 2023: Influence of Surface Types on the Seasonality and Inter-Model Spread of Arctic Amplification in CMIP6, ADVANCES IN ATMOSPHERIC SCIENCES, 40, 2288-2301.  doi: 10.1007/s00376-023-2338-9
    [3] HUANG Wenyu, WANG Bin*, LI Lijuan, DONG Li, LIN Pengfei, YU Yongqiang, ZHOU Tianjun, LIU Li, XU Shiming, XIA Kun, PU Ye, WANG Lu, LIU Mimi, SHEN Si, HU Ning, WANG Yong, SUN Wenqi, and DONG Fang, 2014: Variability of Atlantic Meridional Overturning Circulation in FGOALS-g2, ADVANCES IN ATMOSPHERIC SCIENCES, 31, 95-109.  doi: 10.1007/s00376-013-2155-7
    [4] Yongqi GAO, Helge DRANGE, 2004: The Effect of Diapycnal Mixing on the Ventilation and CFC-11 Uptake in the Southern Ocean, ADVANCES IN ATMOSPHERIC SCIENCES, 21, 755-766.  doi: 10.1007/BF02916372
    [5] Xiao DONG, Jiangbo JIN, Hailong LIU, He ZHANG, Minghua ZHANG, Pengfei LIN, Qingcun ZENG, Guangqing ZHOU, Yongqiang YU, Mirong SONG, Zhaohui LIN, Ruxu LIAN, Xin GAO, Juanxiong HE, Dongling ZHANG, Kangjun CHEN, 2021: CAS-ESM2.0 Model Datasets for the CMIP6 Ocean Model Intercomparison Project Phase 1 (OMIP1), ADVANCES IN ATMOSPHERIC SCIENCES, 38, 307-316.  doi: 10.1007/s00376-020-0150-3
    [6] Fukai LIU, Yiyong LUO, Jian LU, Xiuquan WAN, 2021: The Role of Ocean Dynamics in the Cross-equatorial Energy Transport under a Thermal Forcing in the Southern Ocean, ADVANCES IN ATMOSPHERIC SCIENCES, 38, 1737-1749.  doi: 10.1007/s00376-021-1099-6
    [7] Qizhen SUN, Timo VIHMA, Marius O. JONASSEN, Zhanhai ZHANG, 2020: Impact of Assimilation of Radiosonde and UAV Observations from the Southern Ocean in the Polar WRF Model, ADVANCES IN ATMOSPHERIC SCIENCES, 37, 441-454.  doi: 10.1007/s00376-020-9213-8
    [8] YANG Hui, 2011: The Significant Relationship between the Arctic Oscillation (AO) in December and the January Climate over South China, ADVANCES IN ATMOSPHERIC SCIENCES, 28, 398-407.  doi: 10.1007/s00376-010-0019-y
    [9] Liang ZHAO, Yunwen LIU, Yihui DING, Qingquan LI, Wei DONG, Xinyong SHEN, Wei CHENG, Haoxin YAO, Ziniu XIAO, 2023: The Warm Arctic–Cold Eurasia Pattern and Its Key Region in Winter in CMIP6 Model Simulations, ADVANCES IN ATMOSPHERIC SCIENCES, 40, 2138-2153.  doi: 10.1007/s00376-022-2201-4
    [10] Xiaolei CHEN, Yimin LIU, Guoxiong WU, 2017: Understanding the Surface Temperature Cold Bias in CMIP5 AGCMs over the Tibetan Plateau, ADVANCES IN ATMOSPHERIC SCIENCES, 34, 1447-1460.  doi: 10.1007s00376-017-6326-9
    [11] ZHU Jiawen, ZENG Xiaodong, 2015: Comprehensive Study on the Influence of Evapotranspiration and Albedo on Surface Temperature Related to Changes in the Leaf Area Index, ADVANCES IN ATMOSPHERIC SCIENCES, 32, 935-942.  doi: 10.1007/s00376-014-4045-z
    [12] YANG Yang, REN Rongcai, Ming CAI, RAO Jian, 2015: Attributing Analysis on the Model Bias in Surface Temperature in the Climate System Model FGOALS-s2 through a Process-Based Decomposition Method, ADVANCES IN ATMOSPHERIC SCIENCES, 32, 457-469.  doi: 10.1007/s00376-014-4061-z
    [13] Jiang Hao, Wang Keli, 2001: Analysis of the Surface Temperature on the Tibetan Plateau from Satellite, ADVANCES IN ATMOSPHERIC SCIENCES, 18, 1215-1223.  doi: 10.1007/s00376-001-0035-z
    [14] Jiangbo JIN, Duoying JI, Xiao DONG, Kece FEI, Run GUO, Juanxiong HE, Yi YU, Zhaoyang CHAI, He ZHANG, Dongling ZHANG, Kangjun CHEN, Qingcun ZENG, 2024: CAS-ESM2.0 Dataset for the Carbon Dioxide Removal Model Intercomparison Project (CDRMIP), ADVANCES IN ATMOSPHERIC SCIENCES.  doi: 10.1007/s00376-023-3089-3
    [15] Laura DE LA TORRE, Luis GIMENO, Juan Antonio A\~NEL, Raquel NIETO, 2007: The Role of the Solar Cycle in the Relationship Between the North Atlantic Oscillation and Northern Hemisphere Surface Temperatures, ADVANCES IN ATMOSPHERIC SCIENCES, 24, 191-198.  doi: 10.1007/s00376-007-0191-x
    [16] LIU Peng, QIAN Yongfu, HUANG Anning, 2009: Impacts of Land Surface and Sea Surface Temperatures on the Onset Date of the South China Sea Summer Monsoon, ADVANCES IN ATMOSPHERIC SCIENCES, 26, 493-502.  doi: 10.1007/s00376-009-0493-2
    [17] 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, Dongling ZHANG, 2021: CAS-ESM2.0 Model Datasets for the CMIP6 Flux-Anomaly-Forced Model Intercomparison Project (FAFMIP), ADVANCES IN ATMOSPHERIC SCIENCES, 38, 296-306.  doi: 10.1007/s00376-020-0188-2
    [18] Xiangdong ZHANG, Yunfei FU, Zhe HAN, James E. OVERLAND, Annette RINKE, Han TANG, Timo VIHMA, Muyin WANG, 2022: Extreme Cold Events from East Asia to North America in Winter 2020/21: Comparisons, Causes, and Future Implications, ADVANCES IN ATMOSPHERIC SCIENCES, 39, 553-565.  doi: 10.1007/s00376-021-1229-1
    [19] Xichen LI, Xianyao CHEN, Bingyi WU, Xiao CHENG, Minghu DING, Ruibo LEI, Di QI, Qizhen SUN, Xiaoyu WANG, Wenli ZHONG, Lei ZHENG, Meijiao XIN, Xiaocen SHEN, Chentao SONG, Yurong HOU, 2023: China’s Recent Progresses in Polar Climate Change and Its Interactions with the Global Climate System, ADVANCES IN ATMOSPHERIC SCIENCES, 40, 1401-1428.  doi: 10.1007/s00376-023-2323-3
    [20] Xiaozhen LIN, Chaofan LI, Riyu LU, Adam A. SCAIFE, 2018: Predictable and Unpredictable Components of the Summer East Asia-Pacific Teleconnection Pattern, ADVANCES IN ATMOSPHERIC SCIENCES, 35, 1372-1380.  doi: 10.1007/s00376-018-7305-5

Get Citation+

Export:  

Share Article

Manuscript History

Manuscript received: 02 June 2020
Manuscript revised: 30 July 2020
Manuscript accepted: 20 August 2020
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Surface Temperature Changes Projected by FGOALS Models under Low Warming Scenarios in CMIP5 and CMIP6

    Corresponding author: Shang-Min LONG, smlong@hhu.edu.cn
  • 1. Key Laboratory of Marine Hazards Forecasting, Ministry of Natural Resources, Hohai University, Nanjing 210098, China
  • 2. College of Oceanography, Hohai University, Nanjing 210098, China
  • 3. State key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics and Center for Monsoon System Research, Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100029, China

Abstract: To meet the low warming targets proposed in the 2015 Paris Agreement, substantial reduction in carbon emissions is needed in the future. It is important to know how surface climates respond under low warming targets. The present study investigates the surface temperature changes under the low-forcing scenario of Representative Concentration Pathways (RCP2.6) and its updated version (Shared Socioeconomic Pathways, SSP1-2.6) by the Flexible Global Ocean–Atmosphere–Land System (FGOALS) models participating in phases 5 and 6 of the Coupled Model Intercomparison Project (CMIP5 and CMIP6, respectively). In both scenarios, radiative forcing (RF) first increases to a peak of 3 W m−2 around 2045 and then decreases to 2.6 W m−2 by 2100. Global mean surface air temperature rises in all FGOALS models when RF increases (RF increasing stage) and declines or holds nearly constant when RF decreases (RF decreasing stage). The surface temperature change is distinct in its sign and magnitude between the RF increasing and decreasing stages over the land, Arctic, North Atlantic subpolar region, and Southern Ocean. Besides, the regional surface temperature change pattern displays pronounced model-to-model spread during both the RF increasing and decreasing stages, mainly due to large intermodel differences in climatological surface temperature, ice-albedo feedback, natural variability, and Atlantic Meridional Overturning Circulation change. The pattern of tropical precipitation change is generally anchored by the spatial variations of relative surface temperature change (deviations from the tropical mean value) in the FGOALS models. Moreover, the projected changes in the updated FGOALS models are closer to the multi-model ensemble mean results than their predecessors, suggesting that there are noticeable improvements in the future projections of FGOALS models from CMIP5 to CMIP6.

摘要: 为了实现2015年巴黎协定提出的低增温目标,人类社会需要在未来大量地减少碳排放量,相应的未来情景称之为低增温情景。全球表面温度在低增温背景下如何演变是气候变化中的重要问题。本文研究了在CMIP5中低强迫情景(RCP2.6)和CMIP6中相应升级的低强迫情景(SSP1-2.6)下,4个FGOALS模式预估的未来全球表面温度的变化趋势及相应机制。在RCP2.6和SSP1-2.6两种情景下,外辐射强迫(RF)都是先不断增长,在2045年左右达到3W m-2的峰值,然后在2100年下降到2.6W m-2。在4个FGOALS模式中,全球平均表面气温都是在RF增长的时候(RF上升阶段)不断上升,而在RF下降的时候(RF下降阶段)保持稳定不变或者下降。在RF上升和下降两个不同阶段,表面温度的变化符号和变化强度在陆地、北冰洋、北大西洋副极地地区和南极都有非常显著的差异。此外,表面温度变化的区域分布特征在这两个阶段则都表现出显著的模式间差异性,这主要是因为不同模式在模拟气候态表面温度、冰-反照率反馈机制、自然变率和北大西洋经向翻转环流的变化方面存在非常大的区别。在FGOALS模式中,热带降水的空间分布特征大都受表面温度相对变化(即与热带平均增温强度的偏差)的空间结构所调控。FGOALS模式在低增温情景下得到的预估结果,在CMIP6中相对于在CMIP5中更接近于多模式集合平均的结果。这表明从CMIP5到CMIP6,FGOALS模式对低增温情景下的未来预估模拟有显著提高。

    • The 2015–19 climate report released by the World Meteorological Organization shows that that global mean surface air temperature (GMST) had increased by 1.1°C relative to the preindustrial level (1850–1900) in 2019, leaving a small gap to the 1.5°C warming target proposed in the 2015 Paris Agreement. Recently, the concentrations of greenhouse gases (GHGs), especially CO2, in the atmosphere have risen to record levels, resulting in a global warming trend that will inevitably remain for generations. Global warming has caused significant and severe impacts in terms of ocean temperature and acidity, sea level, ice loss, and extreme events (Church and White, 2011; IPCC, 2013, 2018; Bindoff et al., 2019; Cheng et al., 2019). Therefore, the 2°C and 1.5°C low warming targets are proposed to avert atmospheric GHG concentrations reaching the point of “dangerous anthropogenic interference” with our climate system (Mann, 2009; IPCC, 2018). Indeed, even a 0.5°C decrease from a 2°C to 1.5°C warming level would substantially lower the irreversible damages to climate system and ecosystems (Schaeffer et al., 2012; Schleussner et al., 2017).

      In recent years, concerns about the climate responses under the scenarios anchored by low warming targets (low warming scenarios) have motivated a large body of research and even the release of the special IPCC report on 1.5°C global warming after the 2015 Paris Agreement (IPCC, 2018; Li et al. 2018; Long et al., 2018, 2020; Nangombe et al., 2018; Palter et al., 2018; Qu and Huang, 2018; Zhang et al., 2018; Chen et al., 2019b). Previous studies imply that the spatial distributions and underlying mechanisms of the climate responses may be substantially different between low-forcing and high-forcing scenarios, as the low warming scenarios require much lower or even negative carbon emissions compared to that in the scenarios with moderate or no mitigation efforts (van Vuuren et al., 2011; IPCC, 2013, 2018; Sanderson et al., 2016; Xu and Ramanathan, 2017). This highlights that further studies are needed to deepen our understanding on the dynamics and uncertainty of climate changes under low warming scenarios.

      In phases 5 and 6 of the Coupled Model Intercomparison Project (CMIP5 and CMIP6, respectively), the low-forcing scenario of the Representative Concentration Pathway (RCP2.6) and its updated version (Shared Socioeconomic Pathway, SSP1-2.6) are both categorized as low warming scenarios (Taylor et al., 2012; Eyring et al., 2016). In both RCP2.6 and SSP1-2.6, radiative forcing (RF) is designed to follow nearly the same pathway that first increases to a peak of 3 W m−2 around 2045 and then decreases to 2.6 W m−2 by 2100. Climate responses are substantially different between RF increasing and decreasing stages as the deep ocean warms persistently despite the RF decrease, and thus can significantly shape the trajectories of GMST and surface climate change (Long et al., 2018, 2020). Surface temperature is a key factor in ocean–atmosphere–land interaction, and its future change could greatly affect changes in regional precipitation (Xie et al., 2010; Huang et al., 2013; Long et al., 2016), atmospheric circulation (Ma et al., 2012), and even ocean circulation (Wang et al., 2015; Chen et al., 2019a). However, surface temperatures changes under the low warming scenario in the newly released CMIP6 outputs have not been well studied.

      The present study mainly addresses this issue based on the outputs of the Flexible Global Ocean–Atmosphere–Land System (FGOALS) climate models participating in the CMIPs. The models are developed at the Institute of Atmospheric Physics (IAP)/State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics (LASG), Chinese Academy of Sciences. The FGOALS models in CMIP5 and CMIP6 respectively provide simulations under the RCP2.6 and SSP1-2.6 scenarios (Bao et al., 2013; Lin et al., 2013; He et al., 2019, 2020; Zhou et al., 2020). The state-of-the-art FGOALS models in CMIP6 display substantial improvement in model resolutions, parameterization of physical processes, and tuning method from their predecessor versions in CMIP5, with smaller climate drift and reduced model biases in simulated climatology, seasonal cycles, and climate variability (Guo et al., 2020; Li et al., 2020). It is important to compare the surface temperature changes in the latest two generations of the FGOALS climate models to gain insight into the climate responses and underlying mechanisms under low warming scenarios.

      Therefore, in this study, the surface temperature changes simulated by the FGOALS models under the RCP2.6 and SSP1-2.6 scenarios in CMIP5 (FGOALS-g2 and FGOALS-s2) and CMIP6 (FGOALS-g3 and FGOALS-f3-L) are investigated. We show that there are substantial differences in the patterns of surface temperature change between the RF increasing and decreasing stages over the land, Arctic, North Atlantic (NA) subpolar region, and Southern Ocean. Besides, the pattern of regional surface warming displays striking differences between the two generations of FGOALS models, mainly due to model-to-model differences in the climatological surface temperature field. The surface temperature differences among models can result in large intermodel spread in regional climate responses through different processes like the ice-albedo feedback and Atlantic Meridional Overturning Circulation (AMOC) response. The important implication of this study is that improving the simulation of the surface temperature climatology is helpful in achieving reliable future projections under low warming targets.

      The rest of the paper is organized as follows. Section 2 describes the model outputs and methods. Section 3 presents the global mean responses under the RCP2.6 and SSP1-2.6 scenarios. Sections 4 and 5 investigate the patterns of global and regional surface temperature change, respectively. Section 6 provides a summary, along with some further discussion.

    2.   Model outputs and methods
    • The last two generations of climate system models developed at LASG-IAP are FGOALS2 and FGOALS3, including two parallel subversions in CMIP5 (FGOALS-g2 and FGOALS-s2) and three parallel subversions in CMIP6 (FGOALS-g3, FGOALS-f3-L and FGOALS-f3-H). A detailed description of the model configurations can be found in related literature (Bao et al., 2013; Li et al., 2013, 2020; Zhou and Song, 2014; Guo et al., 2020; Zhou et al., 2020). Only four subversions are analyzed in the present study, as outputs of FGOALS-f3-H are currently not available on the CMIP6 data portal. Each version of the FGOALS models configures a similar coupling framework including oceanic, sea-ice, and land components, with differences mainly in the atmospheric models. Detailed information on these components is provided in Table 1. Note that the FGOALS-g3 (FGOALS-f3-L) model is the updated version of the FGOALS-g2 (FGOALS-s2) model.

      CMIPModelComponent
      OceanSea iceAtmosphereLand
      CMIP5 FGOALS-g2 LICOM2 CICE4 GAMIL2 CLM3
      360 × 196 128 × 60
      30 levels 26 levels
      FGOALS-s2 LICOM2 CSIM5 SAMIL2 CLM3
      360 × 196 128 × 108
      30 levels 26 levels
      CMIP6 FGOALS-g3 LICOM3 CICE4 GAMIL3 CLM4
      360 × 218 188 × 80
      30 levels 26 levels
      FGOALS-f3-L LICOM3 CICE4 FAMIL CLM4
      360 × 218 288 × 180
      30 levels 32 levels

      Table 1.  Model components and corresponding horizontal resolutions of FGOALS models in CMIP5 and CMIP6. The components of all the models are the Finite-volume Atmospheric model (FAMIL), the Spectral Atmospheric Model of IAP/LASG (SAMIL), the LASG/IAP Climate system Ocean Model (LICOM), the Community Land Model (CLM), the Community Sea Ice Model (CSIM), and the Los Alamos sea ice model (CICE). The version number of the component models are labeled after the acronyms.

      The equilibrium climate sensitivity (ECS), defined as the equilibrium temperature under doubled CO2 forcing, is 2.1–4.7 K in CMIP5 models and 1.8–5.6 K in CMIP6 models (Zelinka et al., 2020), which is mainly due to stronger positive cloud feedbacks from decreasing extratropical low-cloud coverage and albedo in models from CMIP6 than CMIP5. In contrast, the ECS is about 3.7 K in FGOALS-g2 and 4.5 K in FGOALS-s2, but decreases to 2.84°C for FGOALS-g3 and 2.98°C for FGOALS-f3 (Zhou et al., 2013, 2020), which might be associated with the differences in model biases and model-simulated internal variability, Arctic climate, and ocean circulation responses.

      The ability of climate models in reproducing the climatology of observations, which is measured by the model bias (i.e., deviation from observation), is an essential metric for evaluating model performance. Figure 1 shows the spatial distribution of annual mean climatology biases in SST and 2 m air temperature (TAS) in FGOALS models for 1979–2005. The ERSST.v5 (Huang et al., 2017) and high-resolution (0.5° × 0.5°) CRU TS4.04 (Harris et al., 2020) data are referenced as observations to calculate the model biases. Generally, the FGOALS models display similar bias patterns in SST and TAS in both CMIPs, with large cold SST bias in the North Pacific, warm SST bias along the eastern coast of subtropical ocean basins and the Southern Ocean, and large cold TAS bias in Eurasia and North America. The maximal and minimal model biases in both SST and TAS also change insignificantly from CMIP5 to CMIP6. However, the global root-mean-square error (RMSE) and global mean value of SST biases are reduced in FGOALS from CMIP5 to CMIP6, with noticeable reduction in the North Pacific, south of Greenland, Southern Hemisphere (SH), eastern subtropical oceans, and Southern Ocean. For TAS bias, there is also a prominent decrease in its magnitude over western Eurasia, North America and Australia. Besides, the updated FGOALS models significantly reduce the TAS biases of their predecessors over regions along the Rocky and Andes Mountains and Himalayas, suggesting an improvement in resolving the effect of topography in CMIP6 models. As a result, the global RMSE and global mean value of TAS biases also decrease, especially for the latter. The evaluation of model biases in FGOALS models suggests that there are significant improvements and changes in the performances of FGOALS-g3 and FGOALS-f3-L from their predecessors.

      Figure 1.  Spatial distribution of annual mean biases in SST (units: °C) and 2 m air temperature (TAS; units: °C) in (a, e) FGOALS-g2, (b, f) FGOALS-g3, (c, g) FGOALS-s2, and (d, h) FGOALS-f3-L, during 1979–2005. The ERSST.v5 and CRU TS4.04 data are referenced as the observations. The maximum value (Max), minimum value (Min), RMSE, and mean value (Mean) of global biases are labeled at the top of each panel.

    • Monthly outputs (1850–2100) of historical simulations and low warming scenarios (RCP2.6 and SSP1-2.6) from the FGOALS models are analyzed. Besides, an additional 15 CMIP6 models (Table 2) are used to compare with their family predecessors in CMIP5. The CMIP5 and CMIP6 multi-model ensemble mean (MME) results are calculated based on the 15 pairs of models without the FGOALS models. The pre-industrial control runs are also used to remove the effect of climate drift in the models. Near-surface air temperature, surface temperature (skin temperature), precipitation, and zonal winds are used in this study. Note that in CMIP outputs, surface skin temperature is equivalent to SST in ice-free ocean. All atmospheric variables are linearly interpolated onto a common grid of 2° latitude × 2° longitude for ease of comparison. Only one member of each model is utilized for the analyses.

      CMIP5CMIP6
      1BCC_CSM1.1(m)BCC-CSM2-MR
      2CanESM2CanESM5
      3CESM1-CAM5CESM2
      4CNRM-CM5CNRM-CM6-1
      5GFDL-CM3GFDL-ESM4
      6GISS-E2-RGISS-E2-1-G
      7HadGEM2-ESHadGEM3-GC31-LL
      8IPSL-CM5A-LRIPSL-CM6A-LR
      9MIROC5MIROC6
      10MIROC-ESMMIROC-ES2L
      11MPI-ESM-LRMPI-ESM1-2-LR
      12MPI-ESM-MRMPI-ESM1-2-HR
      13MRI-CGCM3MRI-ESM2-0
      14NorESM1-MNorESM2-LM
      15NorESM1-MENorESM2-MM

      Table 2.  CMIP models used in the present study.

      As RF first increases to a peak around the year 2045 and then decreases, we separately calculate the linear trends during the RF increasing stage (1850–2050) and RF decreasing stage (2050–2100) to investigate the climate responses during these two distinct periods. The separation point for the trend calculation (2050) is chosen to lead lag the RF inflection point (2045) by only five years because the time scale of the fast response in the ocean mixed layer is 3–5 years (Held et al., 2010). The AMOC index is defined as the maximum value of the meridional stream function at 35°N in the Atlantic.

    3.   Global mean responses
    • Figure 2 displays the temporal evolution of GMST change relative to the preindustrial level (1850–99) through 2100 from the FGOALS models under RCP2.6 and SSP1-2.6. When the RF (gray line) increases, GMST generally rises steadily, despite large multidecadal variations before 1980, in all four models. However, the magnitudes of the increasing trend vary greatly across the models, with FGOALS-g2 (blue line) and FGOALS-g3 (red line) displaying a trend that is about 0.5°C lower than that in their corresponding MME results (black lines) by 2050. Moreover, the GMST changes in FGOALS-g2 and FGOALS-g3 are even close to the minimal value of the range in results from CMIP models after 2050. In contrast, GMST strikingly enhances in FGOALS-s2 (green line), above the MME result of around 0.7°C by 2050, which exceeds the top range of the results from the 15 CMIP5 models (light green shadow). In FGOALS-f3-L, the GMST trajectory is consistent with the CMIP6 MME result, despite some deviations during 1960–2040. It is also worth noting that the FGAOLS models in CMIP6 are closer to the MME results than their predecessors in CMIP5, as shown by both the MME and CMIP models’ ranges.

      Figure 2.  Global mean annual surface air temperature change (units: °C) relative to pre-industrial level (1850–99 mean) in each FGOALS model (colored solid lines) and the CMIP5 and CMIP6 MMEs (black lines), with an 11-year running mean applied. The vertical dashed lines indicate the separation point for trend calculation. Note that the gray lines show the pathway of RF [right-hand axis in (b, d)], and the light green shadows indicate the range of the CMIP models.

      During the RF decreasing stage (2050–2100), as discussed in previous studies (Long et al., 2018, 2020), the ratio between the contributions from fast and slow responses during the RF decreasing stage determines the further trend in GMST. This is because when RF ramps down, the ocean mixed layer would cool fast owing to rapid atmospheric cooling, hence lowering the GMST (fast cooling effect). In contrast, the deep ocean warms slowly but persistently through 2100 and thus reduces the downward heat transfer from the mixed layer, and this would additionally fuel the upper ocean warming and hence increase the GMST (slow warming effect). In both FGOALS-g2 and FGOALS-s2, the fast cooling effect outweighs the slow warming effect, leaving the GMST to slightly and sharply decrease (all significant) during 2050–2100, respectively. In FGOALS-g3 and FGOALS-f3-L, the fast cooling effect and slow warming effect nearly offset each other, leading to an insignificant trend in GMST during the RF decreasing stage, which is also a common feature in most CMIP5 models (Long et al., 2020).

      The increase in GMST meets the 1.5°C warming target in FGOALS-g2 and FGOALS-g3, and is still below the 2°C warming level in FGOALS-f3-L. In FGOALS-s2, the increase in GMST exceeds 2.0°C after 2010 and maximizes at nearly 2.5°C. The large GMST increase in FGOALS-s2 is associated with a lack of the effect from direct aerosol cooling in the atmospheric model (Bao et al., 2013) and large ice-albedo feedback in the sea-ice model (CSIM5), which will be discussed in the next section.

      The abovementioned features in GMST trajectories are more robust over land (red lines in Fig. 3) than ocean (green lines in Fig. 3) in all four models. The land–sea warming contrast (brown lines) also increases before 2050 and then decreases (Figs. 3a and b) or holds nearly constant (Figs. 3c and d) during 2050–2100, consistent with the GMST pathways (black lines). It is worth noting that the range of the CMIP models increases slightly after 2050, suggesting increased intermodel spread and hence model uncertainty during the RF decreasing stage.

      Figure 3.  Global mean annual surface air temperature change (units: °C) over ocean (green lines) and land (red lines) and their difference (land–sea warming contrast, brown lines) relative to the pre-industrial level (1850–99 mean), with an 11-year running mean applied.

      Global mean responses under the low warming scenario suggest that the GMST trajectory may not follow the RF pathway when RF decreases, mainly due to the warming effect from the slow increase in deep ocean temperature (Long et al., 2020). This is also a common feature in FGOALS and other CMIP models (Fig. 2). As regional climate change will largely deviate from the global mean response, we further investigate the spatial pattern of surface temperature change under low warming scenarios in detail.

    4.   Pattern of global surface temperature change
    • Figure 4 shows the linear trend of surface skin temperature (∆TS) under RCP2.6 and SSP1-2.6 during 1850–2050 and 2050–2100. When RF increases, surface warming is generally large over land and polar regions, especially the Arctic. The reduced surface warming or even surface cooling in the NA (i.e., the so-called NA warming hole) and Southern Ocean is robust in all four models, consistent with the MME results (Figs. 5a and b), despite large differences in the warming magnitude and detailed spatial structure, like the location of the NA warming hole. The El Niño-like warming pattern, Indian Ocean Dipole-like warming structure, and reduced subtropical warming in the SH are also prominent in all four models, consistent with the MME results from CMIP5 and CMIP6. These tropical warming patterns also exist in FGOALS-g2 but are not well displayed owing to the small magnitude of warming. Besides, the surface warming is locally enhanced over high mountain regions like the Tibetan Plateau and Rocky and Andes Mountains in all FGOALS models, illustrating the effect of topography in shaping the surface warming structure.

      Figure 4.  Linear trend of annual surface skin temperature (units: °C), which is equivalent to SST over ice-free ocean, during (a–d) 1850–2050 and (e–h) 2050–2100 in the four FGOALS models. The white (black) dots indicate the trend is insignificant (significant) at the 95% level. The global mean value is labeled in the title of each plot.

      Figure 5.  MME linear trends of annual surface temperature (units: °C) during (a, b) 1850–2050 and (c, d) 2050–2100 in (a, c) CMIP5 and (b, d) CMIP6, along with (e–h) their SNRs, defined as the absolute value of MME change divided by the intermodel standard deviation. The magenta and black contours indicate SNRs of 1 and 3, respectively.

      During 2050–2100, corresponding to the significant decreasing trend of GMST (Figs. 2a and b), the surface cooling pattern is prominent in FGOALS-g2 and FGOALS-s2. However, the cooling magnitude differs across regions and is mainly large in the tropics in FGOALS-g2 and in the Northern Hemisphere (NH) mid and high latitudes in FGOALS-s2. Despite the GMST change being negligible during 2050–2100 in these two CMIP6 models, there is still a significant cooling trend in the tropics and warming trend in the NH mid and high latitudes in FGOALS-g3, but overall weak temperature change in FGOALS-f3-L. Besides, there is a broad significant increasing trend in the Southern Ocean in both FGOALS-g3 and FGOALSf-3-L, which is missing in the CMIP5 FGOALS models. This suggest that ∆TS may evolve prominently even under weak GMST change, but with the pattern differing significantly across models. It is worth noting that the ∆TS patterns in the FGOALS models are generally consistent with the MME results during the RF increasing stage, with the global pattern correlation coefficients all exceeding 0.79 (Table 3). However, during the RF decreasing stage, the pattern consistency with the MME results drops dramatically, ranging from −0.11 to 0.5 in the four models, suggesting that there is large model uncertainty in the further changes of surface temperature after 2050. In the CMIP5 and CMIP6 MMEs (Figs. 5c and d), the surface temperature mainly cools over land and warms over the SH oceans, but with the magnitude much reduced compared to that in the FGOALS models during 2050–2100 (Figs. 4e and f). Over the Tibetan Plateau, the surface cooling during the RF decreasing stage is also locally enhanced in most of the FGOALS models (Figs. 4eg), and the CMIPs’ MMEs (Figs. 5c and d), which is similar to the situation during the RF increasing stage and suggests that the surface temperature over that region displays robust responses to RF changes.

      CMIP MME and modelsGlobalTropical oceans
      1850–20502050–21001850–20502050–2100
      CMIP5 MMEFGOALS-g20.86−0.110.73−0.21
      FGOALS-s20.930.380.370.48
      CMIP6 MMEFGOALS-g30.790.500.80−0.04
      FGOALS-f3-L0.93−0.020.59−0.10

      Table 3.  Pattern correlations between FGOALS models and their corresponding MME results. The bold values indicate that the correlations are significant at 95% confidence level.

      To evaluate the role of model uncertainty, which is measured by the intermodel standard deviation (SD), in future projections, we further calculate the rate of model consistency in the sign of MME change and signal-to-noise ratio (SNR) in each grid cell. The former is measured by the rate of models displaying change with the same sign of the MME results and is shown in Figs. 5ad, while the latter is defined as the absolute value of MME change divided by the intermodel SD and is shown in Figs. 5eh. During the RF increasing stage, there is high model consistency in the sign of MME change across the globe, with a model consistency rate below 2/3 only appearing over a very limited area in the NA and Southern Ocean (white dots in Figs. 5a and b). The SNR is generally larger than 3 (black contours in Figs. 5e and f) over most regions during 1850–2050, indicating a robust MME change relative to its intermodel spread, and is only lower than 1 over the NA and Southern Ocean (magenta contours). During 2050–2100, the model consistency largely reduces over most regions, with consistent sign of change among models (black dots) mainly appearing over the Pacific subtropics and land regions with large surface cooling (Figs. 5c and d). Correspondingly, the SNR is much smaller than that during 1850–2050, with a very limited area displaying values larger than 1. The low model consistency and small SNR during the RF decreasing stage suggest that the intermodel spread is much larger than the MME change, which is inevitably associated with the intermodel differences in simulating the natural variability. During 2050–2100 (50 years), as the RF change is weaker and the length of time for trend calculation is much shorter than those during 1850–2050 (200 years), the interference of internal variability in the trend calculation rises consequently. Besides, the warming effect from the deep ocean slow warming also largely offsets the cooling effect from the RF decrease, especially over regions with strong ocean dynamics (Long et al., 2020). All these factors complicate the projections of the surface temperature changes during the RF decreasing stage. As a result, there are large intermodel differences in the ∆TS pattern during 2050–2100, hence lowering the reliability of the MME changes. Indeed, during 1850–2050, the CMIP5 and CMIP6 MME ∆TS patterns are highly similar, suggesting robustness of the projected surface temperature response to the increase in RF. During 2050–2100, despite the GMST change being insignificant in both CMIP5 and CMIP6, the ∆TS pattern diverges substantially over the Arctic, East Asia, North America and NA Ocean.

      Generally, surface temperature responses under low warming scenarios are distinct between the RF increasing and decreasing stages, and vary substantially across models, especially during the RF decreasing stage. Given that the pattern formation mechanisms for ∆TS display large variations in space (Xie et al., 2010), we further investigate the ∆TS pattern in the tropics and three other regions with noticeable local changes (the Arctic, NA subpolar region and Southern Ocean) in detail.

    5.   Patterns of regional surface temperature change
    • Tropical SST plays a key role in regional and global climate, and is an important factor in typhoon/hurricane dynamics, atmospheric convection, and Hadley and Walker cells. The response of tropical SST to global warming is important because it could greatly affect changes in regional precipitation through the so-called “warmer-get-wetter” mechanism (Xie et al., 2010; Huang et al., 2013; Long et al., 2016). Specifically, tropical SST influences atmospheric convection mainly through the relative SST change, i.e., the deviations from tropical mean change, as the tropical convection threshold is mainly determined by the tropical mean SST (Xie et al., 2010; Ma et al., 2012). As a result, precipitation increases over regions with SST warming larger than the tropical mean warming, and vice versa. Note that as the tropical convection threshold also increases as global warming develops (Johnson and Xie, 2010), the relative SST change is defined by subtracting the tropical mean SST change over 30°S–30°N from the local SST change.

      Previous studies suggest that the “warmer-get-wetter” mechanism is prominent in explaining tropical precipitation changes because the dynamic effect from the slowdown of the Walker circulation (Held and Soden, 2006; Vecchi et al., 2006; Vecchi and Soden, 2007; Ma et al., 2012) would largely cancel the thermodynamic effect from the climatological precipitation distribution (Seager et al., 2010; Chadwick et al., 2013). The latter is the so-called “wet-get-wetter” mechanism (Chou and Neelin, 2004; Held and Soden, 2006; Chou et al., 2009). Therefore, we further investigate the patterns of relative SST change in the FGOALS models under RCP2.6 and SSP1-2.6.

      Figure 6 shows the precipitation change (shading) and relative surface temperature change (∆TS*, contours) in the tropical oceans during 1850–2050 and 2050–2100. All four models display an El Niño-like warming structure during 1850–2050, with positive ∆TS* anchoring the increase in precipitation. In contrast, negative ∆TS* prevails over the subtropics and leads to a decrease in precipitation in models except FGOALS-g3. The consistency in the ∆TS* and precipitation change suggests that the “warmer-get-wetter” mechanism works well in FGOALS models.

      Figure 6.  Relative surface temperature change (∆TS*, contours, interval = 0.15°C) and precipitation change (shading) in the tropical oceans during (a–d) 1850–2050 and (e–h) 2050–2100.

      During 2050–2100, the tropical mean SST cooling is noticeable in FGOALS-g2 and FGOALS-g3 (Fig. 4e) but negligible in FGOALS-s2 and FGOALS-f3-L (Fig. 4f). The surface cooling spreads throughout the tropics in FGOALS-g2, especially in the equatorial Pacific, and the SH subtropics in FGOALS-g3, which also leads to a pronounced decrease in precipitation (Figs. 6e and f). However, precipitation still increases over several regions, mainly around the warm pool. This might be associated with the direct effect of the decrease in CO2 during 2050–2100. Previous studies have revealed that, in response to an increase in atmospheric CO2 concentration, global mean precipitation tends to decrease for several decades owing to fast tropospheric adjustment of the atmosphere, which is most prominent over the tropics because of Walker cell adjustment (Mitchell et al., 1987; Allen and Ingram, 2002; Lambert and Webb, 2008; Andrews et al., 2010; Wu et al., 2010; Kamae and Watanabe, 2013). Likewise, a decrease in atmospheric CO2 concentration would drive an increase in global mean precipitation and hence tropical circulation. Therefore, the effects from surface cooling and reduced CO2 jointly shape the pattern of tropical precipitation change during the RF decreasing stage. This is different from the case during the RF increasing stage, as the direct effect of increased CO2 on precipitation change is overwhelmed by the strong surface warming effect. Despite that the SST cools or weakly warms during the RF decreasing stage, the relative SST change still exerts strong control over the spatial structures in precipitation change in all FGOALS models (Figs. 6eh). Precipitation generally increases or decreases slightly over regions with SST warming or SST cooling smaller than the tropical mean, and decreases substantially over regions with large SST cooling.

    • The pronounced Arctic warming (Arctic amplification) and reduced Southern Ocean warming are common regional patterns in all FGOALS models and CMIP MME results during the RF increasing stage (Figs. 4 and 5). The temporal evolutions of the area-weighted mean surface temperature over these regions are presented in Fig. 7.

      Figure 7.  Area-weighted mean annual surface temperature (units: °C) in the (a, b) tropical oceans (20°S–20°N), (c, d) Arctic (60°–90°N), and (e, f) Southern Ocean (65°–50°S), in FGOALS-g2 (blue lines), FGOALS-g3 (red lines), FGOALS-s2 (green lines), and FGOALS-f3-L (magenta lines), with an 11-year running mean applied. The colored dashed lines are the least-squares linear fitting lines for 1850–2050 and 2050–2100.

      Arctic amplification is the most pronounced signal in the long-term climate change during recent decades and under future warming scenarios (Holland and Bitz, 2003; Serreze et al., 2009; Purkey and Johnson, 2010; Serreze and Barry, 2011; IPCC, 2013; Cohen et al., 2014). It involves or interacts with varies processes, including changes in atmospheric and oceanic circulation (Graversen et al., 2008; Chylek et al., 2009; Simmonds and Keay, 2009) and feedbacks associated with temperature in the low latitudes (Pithan and Mauritsen, 2014), cloud (Schweiger et al., 2008), water vapor (Francis and Hunter, 2007), snow, and sea ice (Winton, 2006; Kumar et al., 2010; Screen and Simmonds, 2010). Here, we mainly discuss the role of the ice-albedo feedback, which arises from the fact that ice has a much higher albedo and can reflect more solar radiation than land or water surfaces. An initial surface warming over ice surfaces would lead to ice melting and exposure of the underlying land or water to the atmosphere, and the reduced surface albedo that results would allow a greater absorption of surface heat. This enhances surface warming and leads to more ice melting, forming the positive feedback loop that is key in amplifying the initial Arctic surface temperature anomaly.

      The climatological Arctic surface temperature is about 2°C–3°C colder in FGOALS-g2 and FGOALS-g3 than in FGOALS-s2 and FGOALS-f3-L by 1900 (Figs. 7c and d), suggesting that models with a cold Arctic surface tend to display a large ice thickness that is supposed to result in weak Arctic amplification (Holland and Bitz, 2003). The thick ice would require more heat for melting and hence delay the trigger and development of the ice-albedo feedback. As a result, during 1850–2050, the Arctic warming is much smaller in FGOALS-g2 and FGOALS-g3 than in FGOALS-s2 and FGOALS-f3-L. Indeed, the area-averaged annual mean ice thickness is 3.24 m for FGOALS-g2 and 1.81 m for FGOALS-s2 in the climatology (Song et al., 2014), confirming the relationship between surface temperature and ice thickness.

      It is worth noting that the Arctic warming is the largest in FGOALS-s2, accompanied by the highest Arctic surface temperature in the climatology. As reported in Bao et al. (2013), there is a lack of the direct aerosol effect in FGOALS-s2. As aerosol emissions are largest in the NH, the cooling effect of aerosol would be underestimated in the NH, resulting a warmer Arctic surface and hence thinner ice in FGOALS-s2 than the other models. Therefore, the Arctic amplification is also the most prominent in FGOALS-s2, with a warming close to 7°C relative to the pre-industrial level by 2050. When RF decreases, the Arctic surface temperature displays a sharp decreasing trend of −3°C during 2050–2100, which may be triggered by the shutdown of AMOC in FGOALS-s2, which is discussed in the next section. This suggests that the positive ice-albedo feedback is highly efficient in amplifying the initial temperature anomaly in the sea-ice model of FGOALS-s2. In contrast, the Arctic temperature continues to rise significantly in FGOALS-g3 during the RF decreasing stage, possibly because the thick ice has not melted too much over some regions during the RF increasing stage. In FGOALS-g2 and FGOALS-f3-L, further changes in Arctic surface temperature are relatively small and mainly follow the GMST pathways.

      During the RF increasing stage, there is broad reduced surface warming (i.e., warming trend smaller than the global mean) at 65°S–50°S (Figs. 4ad), and this is accompanied by robust strengthening of the westerlies south of 45°S in all four models (Figs. 8ad), which would drive strong Southern Ocean upwelling and equatorward Ekman transport that jointly suppress the Southern Ocean surface temperature increase (Armour et al., 2016). As a result, the area-weighted-mean Southern Ocean surface temperature increases much slower than most other regions (Figs. 7e and f). During the RF decreasing stage, surface temperature decreases slightly in FGOALS-g2 but continues to increase in the other models, especially in FGOALS-s2 (green line in Fig. 7f). This is distinct from the GMST trajectory during that period. As discussed in previous studies (Held et al., 2010; Long et al., 2014, 2018, 2020), the deep ocean warms persistently despite the decrease in RF, and would additionally increase the surface temperature. In the Southern Ocean high latitudes, as the deep ocean gradually warms (Long et al., 2020), the seawater upwelled to the surface during the RF decreasing stage is also warmer than that during the RF increasing stage. This would weaken the cooling effect from the upwelling and thus fuel the further increase in surface temperature over the Southern Ocean.

      Figure 8.  Linear trend of annual 850 hPa zonal wind (units: m s−1) during (a–d) 1850–2050 and (e–h) 2050–2100 in the four FGOALS models. The white (black) dots indicate the trend is insignificant (significant) at the 95% level.

    • The subpolar NA is a region of complex ocean dynamics. It plays a key role in regional and global climate variability at long time scales, such as Atlantic Multidecadal Oscillation (AMO), as the North Atlantic Deep Water (NADW) forms over there. Under global warming, the NA subpolar region features reduced surface warming relative to the global mean warming level, or even a cooling trend, in the 20th century, i.e., the NA warming hole (Drijfhout et al., 2012). This spatial structure is also evident in all four FGOALS models under low warming scenarios (Figs. 4ad) and is common in CMIP models (IPCC, 2013; Sgubin et al., 2017). However, the NA warming hole appears at different locations in FGOALS models (blue boxes in Figs. 4ad), mainly due to the differences in the NADW formation regions and responses of NA subpolar gyres.

      The NA subpolar region displays the largest SST anomaly in the spatial pattern of the AMO (Li et al., 2020). Mechanisms for the formation of the NA warming hole have been studied in recent decades, including the role of AMOC, subpolar gyre adjustment, and local convection change (Kim and An, 2013; Sgubin et al., 2017; Menary and Wood, 2018; Keil et al., 2020). All these factors are also directly related to the AMO, as the NA subpolar region is also a key region in the AMO. Figure 9 shows the area-weighted mean surface temperature change over the NA warming hole region and corresponding AMO index, defined as the detrended area-weighted mean SST anomaly over (0°–60°N, 80°W–0°), in each model. The warming hole indices display substantial decadal and multidecadal variability and hence a weak warming trend or even cooling trend in the FGOALS models. Indeed, the detrended warming hole indices correlate well with the AMO indices in all four models between 1850 and 2100, with correlation coefficients ranging from 0.53 in FGOALS-f3-L to 0.85 in FGOALS-g2. Therefore, in the FGOALS models, the AMO is important in influencing the surface temperature change over the NA subpolar region or the NA warming hole, which has not been well investigated in previous studies (Drijfhout et al., 2012; Kim and An, 2013; Sgubin et al., 2017; Menary and Wood, 2018) as it is hard to disentangle the causality of AMO and AMOC change.

      Figure 9.  Area-weighted mean annual surface temperature change (unit: °C) in the NA warming hole region (WH index), with an 11-year running mean applied. The regions are marked as blue boxes in Fig. 3, which are (45°–60°N, 65°–30°W) in FGOALS-g2 and FGOALS-g3 and (40°–60°N, 40°–10°W) in FGOALS-s2 and FGOALS-f3-L. The brown lines indicate the AMO index, defined as the detrended area-weighted mean SST over (0°–60°N, 80°W–0°). The change is calculated as the difference between pre-industrial control simulations and simulations from historical and future scenarios. The colored dashed lines are the least-squares linear fitting lines for 1850–2050 and 2050–2100.

      It is worth noting that in FGOALS-g3, the NA warming hole index displays much larger temporal variations than in the other three models. This is because the FGOALS-g3 model greatly overestimates the AMO SST anomalies over the west of Greenland and the Labrador Sea in both the historical and pre-industrial control runs (Li et al., 2020). The large temporal variations in the warming hole index substantially affect the linear trend, as the inflection point (2050) for trend calculation is right in the valley of the AMO cold phase. This leads to a large cooling trend during 1850–2050 and substantial warming trend during 2050–2100 over the NA warming hole region.

      In FGOALS-s2, the NA subpolar surface temperature cools rapidly from 2020 to 2100, mainly due to the rapid weakening of the AMOC (Fig. 10d), which is absent in the other three models. The reduced northward heat transport due to the AMOC slowdown would lead to large cooling over the NH, which may trigger the ice-albedo feedback to cause remarkable Arctic cooling during 2050–2100, consistent with the results in Fig. 7d. Furthermore, the AMOC is nearly shut down in FGOALS-s2, mainly because of the large freshwater input over the NA subpolar regions from the striking Arctic warming-induced ice melting. The AMOC displays large intensity in the climatology (values labeled in the legend for Figs. 10c and d), and hence is more stable under external forcing in the other three models than in FGOALS-s2. Climatologically, the AMOC intensity is directly related to the surface temperature over the NA subpolar region, as cold surface water promotes deep water formation, and vice versa. In addition, from CMIP5 to CMIP6, the FGOALS models largely reduce the warm SST biases south of Greenland (Figs. 1ad) and hence increase the climatological AMOC intensity (Figs. 10c and d). This further highlights that to achieve reliable AMOC projections under future warming scenarios, improving the model ability in reproducing the observed climatological SST over the NA subpolar region is the key. The surface temperature over the NA subpolar regions may also work as a useful “emergency constraint” in correcting AMOC projections in climate models, as the period of direct AMOC observations is rather short.

      Figure 10.  Area-weighted mean annual surface temperature (unit: °C) in the (a, b) NA warming hole region and (c, d) AMOC index (unit: Sv), with an 11-year running mean applied. Note that the curves display the annual mean AMOC index in each year, but not the AMOC change. The colored dashed lines are the least-squares linear fitting lines for 1850–2050 and 2050–2100.

    6.   Summary and discussion
    • The surface temperature changes under low warming scenarios projected by two generations of FGOALS models from CMIP5 and CMIP6 have been investigated in this paper. In these scenarios, RF first increases before 2045 and then decreases, with GMST generally increasing in all FGOALS models during the RF increasing stage and declining or holding nearly constant during the RF decreasing stage. The consistency in GMST change with the MME results is improved in the updated versions of the FGOALS models in CMIP6 from their predecessors, especially from FGOALS-s2 to FGOAMS-f3-L.

      Furthermore, our results show that surface warming patterns are distinct between the RF increasing and decreasing stages, and display large differences across models. The most prominent differences appear over the land, Arctic, NA subpolar region, and Southern Ocean. The pattern consistency between the surface temperature change projected by the FGOALS models and the MME results is high during the RF increasing stage but low during the RF decreasing stage. This is mainly due to the relative short period for trend calculation in the latter, i.e., the natural variability plays an important role in obscuring the externally forced change during that period. Indeed, the CMIP5 and CMIP6 MMEs also display large differences over the NH mid and high latitudes during the RF decreasing stage, suggesting large model uncertainty in projecting the climate responses under weak GMST change. Tropical precipitation change is mainly anchored by the relative sea surface warming pattern during the RF increasing stage, which is a common feature in the FGOALS models. The relative SST change also exerts strong control over the pattern of precipitation change during the RF decreasing stage in all FGOALS models despite them displaying broad SST cooling or weak SST warming in the tropics.

      The important implication of the present study is that by comparing the surface temperature responses in different FGOALS models, the diversity in the projected surface temperature change across models can largely be traced back to the intermodel differences in climatological surface temperature. In the Arctic, models with a warmer surface temperature would result in thinner ice and hence more active ice-albedo feedback and larger Arctic warming during the RF increasing stage than in the other models. The shutdown of AMOC in FGOALS-s2 is suggested to be tightly associated with the striking Arctic amplification-induced freshwater forcing over the NA subpolar region. Besides, the surface temperature over the NA subpolar region determines the magnitude of deep water formation and hence climatological AMOC intensity. Therefore, intermodel differences in NA subpolar climatological SST could lead to distinct AMOC responses and associated climatic effects globally across models in the future. Moreover, in the tropics, a 0.5°C difference in climatological SST could lead to substantially different patterns of precipitation change across the FGOALS models (Fig. 6). Consequently, a more realistic climatological surface temperature field is essential in reducing the intermodel spread of model-projected future changes, especially under low warming scenarios, where the RF is relatively weaker than under medium- and high-emissions scenarios. Further studies are needed to systematically investigate the detailed process involving the surface temperature changes under low warming scenarios in the state-of-the-art climate models from CMIP6.

      Throughout the analyses of projected future surface temperature changes under low warming scenarios, we can see striking differences between both global mean and regional changes from the four FGOALS models during both the RF increasing and decreasing stages. Such differences across models might be associated with complicated process, such as the differences in model biases, cloud simulation, RF data (Nie et al., 2019), AMOC stability, land and sea-ice models, and parameterization schemes. For example, the ECS is substantially reduced in the updated FGOALS models compared to their predecessors, which may relate to the improved simulation in climatological surface temperature, especially in terms of the warm bias over the Arctic and the NA subpolar region and cold bias over land. The model-simulated natural variability is also an important factor influencing the intermodel spread in future projections. As shown in Fig. 9, the magnitude and phase transition of the AMO differ substantially across the FGOALS models, which may involve the intermodel differences in climatological SST over the NA subpolar gyre, storm track, and climate feedbacks, etc. It is worth noting that, despite the locations of the NA warming hole being nearly identical in models from the same group (i.e., FGOALS-g2 and FGOALS-g3, FGOALS-s2 and FGOALs-f3-L), the temporal evolution of the warming hole and AMOC indices are distinct, especially during the RF decreasing stage, and display large multidecadal variability that is dominated by the AMO. As GMST stabilizes during the RF decreasing stage, model-to-model differences in simulating the coherent natural variability mode would lead to substantially different trend changes due to the relatively short period and relatively weak RF change. Therefore, improving the ability of climate models to reproduce the natural variability mode will help to improve the reliability of model-projected future change under low warming scenarios.

      Acknowledgements. We acknowledge the WCRP Working Group on Coupled Modeling, which is responsible for CMIP, and the climate modeling groups for producing and making available their model outputs. All data are available at http://pcmdi9.llnl.gov/. This work was supported by the Natural Science Foundation of China (Grant Nos. 41706026 and 41831175), the National Key Research and Development Program of China (Grant Nos. 2017YFA0604600, 2016YFA0601804 and 2018YFA0605702), the Fundamental Research Funds for the Central Universities (Grant Nos. 2018B04814 and 2018B03114), the open fund of the State Key Laboratory of Satellite Ocean Environment Dynamics, Second Institute of Oceanography (Grant No. QNHX 1808), and the National Basic Research Program of China (Grant No. 2012CB955600).

Reference

Catalog

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return