CO2 climate sensitivity and snow-sea-ice albedo parameterization in an atmospheric GCM coupled to a mixed-layer ocean model

The snow-sea-ice albedo parameterization in an atmospheric general circulation model (GCM), coupled to a simple mixed-layer ocean and run with an annual cycle of solar forcing, is altered from a version of the same model described by Washington and Meehl (1984). The model with the revised formulation is run to equilibrium for 1 × CO₂ and 2 × CO₂ experiments. The 1 ×CO₂ (control) simulation produces a global mean climate about 1° warmer than the original version, and sea-ice extent is reduced. The model with the altered parameterization displays heightened sensitivity in the global means, but the geographical patterns of climate change due to increased carbon dioxide (CO₂) are qualitatively similar. The magnitude of the climate change is affected, not only in areas directly influenced by snow and ice changes but also in other regions of the globe, including the tropics where sea-surface temperature, evaporation, and precipitation over the oceans are greater. With the less-sensitive formulation, the global mean surface air temperature increase is 3.5 °C, and the increase of global mean precipitation is 7.12%. The revised formulation produces a globally averaged surface air temperature increase of 4.04 °C and a precipitation increase of 7.25%, as well as greater warming of the upper tropical troposphere. Sensitivity of surface hydrology is qualitatively similar between the two cases with the larger-magnitude changes in the revised snow and ice-albedo scheme experiment. Variability of surface air temperature in the model is comparable to observations in most areas except at high latitudes during winter. In those regions, temporal variation of the sea-ice margin and fluctuations of snow cover dependent on the snow-ice-albedo formulation contribute to larger-than-observed temperature variability. This study highlights an uncertainty associated with results from current climate GCMs that use highly parameterized snow-sea-ice albedo schemes with simple mixed-layer ocean models.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Future surface mass balance of the Antarctic ice sheet and its influence on sea level change, simulated by a regional atmospheric climate model

A regional atmospheric climate model with multi-layer snow module (RACMO2) is forced at the lateral boundaries by global climate model (GCM) data to assess the future climate and surface mass balance (SMB) of the Antarctic ice sheet (AIS). Two different GCMs (ECHAM5 until 2100 and HadCM3 until 2200) and two different emission scenarios (A1B and E1) are used as forcing to capture a realistic range in future climate states. Simulated ice sheet averaged 2 m air temperature (T_2m) increases (1.8–3.0 K in 2100 and 2.4–5.3 K in 2200), simultaneously and with the same magnitude as GCM simulated T_2m. The SMB and its components increase in magnitude, as they are directly influenced by the temperature increase. Changes in atmospheric circulation around Antarctica play a minor role in future SMB changes. During the next two centuries, the projected increase in liquid water flux from rainfall and snowmelt, together 60–200 Gt year^?1, will mostly refreeze in the snow pack, so runoff remains small (10–40 Gt year^?1). Sublimation increases by 25–50 %, but remains an order of magnitude smaller than snowfall. The increase in snowfall mainly determines future changes in SMB on the AIS: 6–16 % in 2100 and 8–25 % in 2200. Without any ice dynamical response, this would result in an eustatic sea level drop of 20–43 mm in 2100 and 73–163 mm in 2200, compared to the twentieth century. Averaged over the AIS, a strong relation between $\Updelta$ SMB and $\Updelta\hbox{T}_{2{\rm m}}$ of 98 ± 5 Gt w.e. year^?1 K^?1 is found.     

The transport of CO2-induced warming into the ocean: an analysis of simulations by the OSU coupled atmosphere-ocean general circulation model

The OSU global coupled atmosphere-ocean general circulation model has been used to investigate a 2xCO₂-induced climate change. A previous analysis of the simulated 2xCO₂–1xCO₂ temperature differences showed that the CO₂-induced warming penetrated into the ocean and thereby caused a delay in the equilibration of the climate system with an estimatede-folding time of 50–75 years. The objective of the present study is to determine by what pathways and through which physical processes the simulated ocean general circulation produces the penetration of the CO₂-induced warming into the ocean.A global-mean oceanic heat budget analysis shows that the ocean gains heat at a rate of 3 W/m² due to the CO₂ doubling, and that this heat penetrates downward into the ocean predominantly through the reduction in the convective overturning. A zonal-mean oceanic heat budget analysis shows that the surface warming increases from the tropics toward the midlatitudes of both hemispheres and gradually penetrated into the deeper ocean, with a greater penetration in the subtropics and midlatitudes than in the equatorial region. The zonal-mean heat budget analysis also shows that the CO₂-induced warming of the ocean occurs predominantly through the down-ward transport of heat, with the meridional heat flux being only of secondary importance. In the tropics the penetration of the CO₂-induced heating is minimized by the upwelling of cold water. In the subtropics the heating is transported down-ward more readily by the downwelling existing there. In the high latitudes the suppressed convection plays the dominant role in the downward penetration of the CO₂-induced heating. The latter result should be considered as tentative, however, as the ocean component of the coupled model employed a prescribed surface salinity field and did not include the mechanism of brine rejection when sea water freezes into sea ice.     

Verification of model simulated mass balance, flow fields and tabular calving events of the Antarctic ice sheet against remotely sensed observations

The Antarctic ice sheet (AIS) has the greatest potential for global sea level rise. This study simulates AIS ice creeping, sliding, tabular calving, and estimates the total mass balances, using a recently developed, advanced ice dynamics model, known as SEGMENT-Ice. SEGMENT-Ice is written in a spherical Earth coordinate system. Because the AIS contains the South Pole, a projection transfer is performed to displace the pole outside of the simulation domain. The AIS also has complex ice-water-granular material-bedrock configurations, requiring sophisticated lateral and basal boundary conditions. Because of the prevalence of ice shelves, a ‘girder yield’ type calving scheme is activated. The simulations of present surface ice flow velocities compare favorably with InSAR measurements, for various ice-water-bedrock configurations. The estimated ice mass loss rate during 2003–2009 agrees with GRACE measurements and provides more spatial details not represented by the latter. The model estimated calving frequencies of the peripheral ice shelves from 1996 (roughly when the 5-km digital elevation and thickness data for the shelves were collected) to 2009 compare well with archived scatterometer images. SEGMENT-Ice’s unique, non-local systematic calving scheme is found to be relevant for tabular calving. However, the exact timing of calving and of iceberg sizes cannot be simulated accurately at present. A projection of the future mass change of the AIS is made, with SEGMENT-Ice forced by atmospheric conditions from three different coupled general circulation models. The entire AIS is estimated to be losing mass steadily at a rate of ~120 km³/a at present and this rate possibly may double by year 2100.     

Air-sea Interaction of Typhoon Sinlaku （2002） Simulated by the Canadian MC2 Model

Three experiments for the simulation of typhoon Sinlaku (2002) over the western North Pacific are performed in this study by using the Canadian Mesoscale Compressible Community (MC2) atmospheric model. The objective of these simulations is to investigate the air-sea interaction during extreme weather conditions, and to determine the sensitivity of the typhoon evolution to the sea surface temperature (SST)cooling induced by the typhoon. It is shown from the three experiments that the surface heat fluxes have a substantial influence on the slow-moving cyclone over its lifetime. When the SST in the East China coastal ocean becomes 1℃ cooler in the simulation, less latent heat and sensible heat fluxes from the underlying ocean to the cyclone tend to reduce the typhoon intensity. The cyclone is weakened by 7 hPa at the time of its peak intensity. The SST cooling also has impacts on the vertical structure of the typhoon by weakening the warm core and drying the eye wall. With a finer horizontal resolution of (1/6)°×(1/6)°, the model produces higher surface wind, and therefore more surface heat fluxes are emitted from the ocean surface to the cyclone, in the finer-resolution MC2 grid compared with the relatively lower resolution of 0.25°×0.25°MC2 grid.     

Sensitivity of a Greenland ice sheet model to ice flow and ablation parameters: consequences for the evolution through the last climatic cycle

Sensitivity experiments are conducted to test the influence of poorly known model parameters on the simulation of the Greenland ice sheet by means of a three dimensional numerical model including the mechanical and thermal processes within the ice. Two types of experiments are performed: steady-state climatic conditions and simulations over the last climatic cycle with a climatic forcing derived from the GRIP record. The experiments show that the maximum altitude of the ice sheet depends on the ice flow parameters (deformation and sliding law coefficients, geothermal flux) and that it is low when the ice flow is fast. On the other hand, the maximum altitude is not sensitive to the ablation strength and consequently during the climatic cycle it is driven by changes in accumulation rate. The ice sheet extension shows the opposite sensitivity: it is barely affected by ice flow velocity and the ice covered area is smaller for large ablation coefficients. For colder climates, when there is no ablation, the ice sheet extension depends on the sea level. An interesting result is that the variations with time of the altitude at the ice divide (Summit) do not depend on the parameters we tested. The present modelled ice sheets resulting from the climatic cycle experiments are compared with the present measured ice sheet in order to find the set of parameters that gives the best fit between modelled and measured geometry. It seems that, compared to the parameter set most commonly used, higher ablation rate coefficents must be used. Received: 19 September 1995 / Accepted: 30 May 1996     

The sea ice mass budget of the Arctic and its future change as simulated by coupled climate models

Arctic sea ice mass budgets for the twentieth century and projected changes through the twenty-first century are assessed from 14 coupled global climate models. Large inter-model scatter in contemporary mass budgets is strongly related to variations in absorbed solar radiation, due in large part to differences in the surface albedo simulation. Over the twenty-first century, all models simulate a decrease in ice volume resulting from increased annual net melt (melt minus growth), partially compensated by reduced transport to lower latitudes. Despite this general agreement, the models vary considerably regarding the magnitude of ice volume loss and the relative roles of changing melt and growth in driving it. Projected changes in sea ice mass budgets depend in part on the initial (mid twentieth century) ice conditions; models with thicker initial ice generally exhibit larger volume losses. Pointing to the importance of evolving surface albedo and cloud properties, inter-model scatter in changing net ice melt is significantly related to changes in downwelling longwave and absorbed shortwave radiation. These factors, along with the simulated mean and spatial distribution of ice thickness, contribute to a large inter-model scatter in the projected onset of seasonally ice-free conditions.     

Evaluation of the North Atlantic Oscillation as simulated by a coupled climate model

 The realism of the Hadley Centre’s coupled climate model (HadCM2) is evaluated in terms of its simulation of the winter North Atlantic Oscillation (NAO), a major natural mode of the Northern Hemisphere atmosphere that is currently the subject of considerable scientific interest. During 1400 y of a control integration with present-day radiative forcing levels, HadCM2 exhibits a realistic NAO associated with spatial patterns of sea level pressure, synoptic activity, temperature and precipitation anomalies that are very similar to those observed. Spatially, the main model deficiency is that the simulated NAO has a teleconnection with the North Pacific that is stronger than observed. In a temporal sense the simulation is compatible with the observations if the recent observed trend (from low values in the 1960s to high values in the early 1990s) in the winter NAO index (the pressure difference between Gibraltar and Iceland) is ignored. This recent trend is, however, outside the range of variability simulated by the control integration of HadCM2, implying that either the model is deficient or that external forcing is responsible for the variation. It is shown, by analysing two ensembles, each of four HadCM2 integrations that were forced with historic and possible future changes in greenhouse gas and sulphate aerosol concentrations, that a small part of the recent observed variation may be a result of anthropogenic forcing. If so, then the HadCM2 experiments indicate that the anthropogenic effect should reverse early next century, weakening the winter pressure gradient between Gibraltar and Iceland. Even combining this anthropogenic forcing and internal variability cannot explain all of the recent observed variations, indicating either some model deficiency or that some other external forcing is partly responsible. Received: 20 August 1998 / Accepted: 12 May 1999     

Influence of the North American Dipole on ENSO onset as simulated by a coupled ocean-Atmosphere model

北美偶极子(NAD)是热带北大西洋西部和北美东北部的南北向海平面气压异常偶极型模态.以往的观测研究表明,NAD可以有效地影响ENSO事件的爆发.本文利用全球耦合模式FGOALS-g2,评估了NAD与ENSO的关系.结果表明,该模式能较好地重现NAD模态.进一步的分析验证了冬季NAD可以通过强迫冬末春初副热带东北太平洋上空的反气旋和暖海温的出现,在随后的冬季触发El Ni?o事件.此外,在同化NAD实验中,发生El Ni?o事件的概率增加了将近一倍.相比之下,NAO未能在副热带东北太平洋上空引起表面风和海温的异常,因而不能有效地激发次年冬季ENSO事件.     

Near future (2016-40) summer precipitation changes over China as projected by a regional climate model (RCM) under the RCP8.5 emissions scenario: Comparison between RCM downscaling and the driving GCM

Multi-decadal high resolution simulations over the CORDEX East Asia domain were performed with the regional climate model RegCM3 nested within the Flexible Global Ocean-Atmosphere-Land System model, Grid-point Version 2 (FGOALS-g2). Two sets of simulations were conducted at the resolution of 50 km, one for present day (1980–2005) and another for near-future climate (2015–40) under the Representative Concentration Pathways 8.5 (RCP8.5) scenario. Results show that RegCM3 adds value with respect to FGOALS-g2 in simulating the spatial patterns of summer total and extreme precipitation over China for present day climate. The major deficiency is that RegCM3 underestimates both total and extreme precipitation over the Yangtze River valley. The potential changes in total and extreme precipitation over China in summer under the RCP8.5 scenario were analyzed. Both RegCM3 and FGOALS-g2 results show that total and extreme precipitation tend to increase over northeastern China and the Tibetan Plateau, but tend to decrease over southeastern China. In both RegCM3 and FGOALS-g2, the change in extreme precipitation is weaker than that for total precipitation. RegCM3 projects much stronger amplitude of total and extreme precipitation changes and provides more regional-scale features than FGOALS-g2. A large uncertainty is found over the Yangtze River valley, where RegCM3 and FGOALS-g2 project opposite signs in terms of precipitation changes. The projected change of vertically integrated water vapor flux convergence generally follows the changes in total and extreme precipitation in both RegCM3 and FGOALS-g2, while the amplitude of change is stronger in RegCM3. Results suggest that the spatial pattern of projected precipitation changes may be more affected by the changes in water vapor flux convergence, rather than moisture content itself.     

Comparison of the last interglacial climate simulated by a coupled global model of intermediate complexity and an AOGCM

The climate at the Last Interglacial Maximum (125 000 years before present) is investigated with the atmosphere-ocean general circulation model ECHAM-1/LSG and with the climate system model of intermediate complexity CLIMBER-2. Comparison of the results of the two models reveals broad agreement in most large-scale features, but also some discrepancies. The fast turnaround time of CLIMBER-2 permits one to perform a number of sensitivity experiments to (1) investigate the possible reasons for these differences, in particular the impact of different freshwater fluxes to the ocean, (2) analyze the sensitivity of the results to changes in the definition of the modern reference run concerning CO₂ levels (preindustrial versus “present”), and (3) estimate the role of vegetation in the changed climate. Interactive vegetation turns out to be capable of modifying the initial climate signals significantly, leading especially to warmer winters in large parts of the Northern Hemisphere, as indicated by various paleodata. Differences due to changes in the atmospheric CO₂ content and due to interactive vegetation are shown to be at least of the same order of magnitude as differences between the two completely different models, demonstrating the importance of careful experimental design.     

The interactive climate and vegetation along the Pole-Equator belts simulated by a global coupled model

The interaction between climate and vegetation along four Pole-Equator-Pole （PEP） belts were explored using a global two-way coupled model, AVIM-GOALS, which links the ecophysiological processes at the land surface with the general circulation model （GCM）. The PEP belts are important in linking the climate change with the variation of sea and land, including terrestrial ecosystems. Previous PEP belts studies have mainly focused on the paleoclimate variation and its reconstruction. This study analyzes and discusses the interaction between modern climate and vegetation represented by leaf area index （LAI） and net primary production （NPP）. The results show that the simulated LAI variation, corresponding to the observed LAI variation, agrees with the peak-valley variation of precipitation in these belts. The annual mean NPP simulated by the coupled model is also consistent with PIK NPP data in its overall variation trend along the four belts, which is a good example to promote global ecological studies by coupling the climate and vegetation models. A large discrepancy between the simulated and estimated LAI emerges to the south of 15°N along PEP 3 and to the south of 18°S in PEP 1S, and the discrepancy for the simulated NPP and PIK data in the two regions is relatively smaller in contrast to the LAI difference. Precipitation is a key factor affecting vegetation variation, and the overall trend of LAI and NPP corresponds more obviously to precipitation variation than temperature change along most parts of these PEP belts.     

Comparison of Surface Wind Stress Anomalies over the Tropical Pacific Simulated by an AGCM and by a Simple Atmospheric Model

Ｃｏｍｐａｒｉｓｏｎ　ｏｆ　Ｓｕｒｆａｃｅ　Ｗｉｎｄ　Ｓｔｒｅｓｓ　Ａｎｏｍａｌｉｅｓ　ｏｖｅｒ　ｔｈｅ　Ｔｒｏｐｉｃａｌ　Ｐａｃｉｆｉｃ　Ｓｉｍｕｌａｔｅｄ　ｂｙ　ａｎ　ＡＧＣＭ　Ａｎｄ　ｂｙ　ａ　Ｓｉｍｐｌｅ　Ａｔｍｏｓｐｈｅｒｉｃ　ＭｏｄｅｌＮ...     

夏季亚洲-太平洋涛动的耦合模式模拟

亚洲-太平洋涛动是夏季欧亚大陆东部(15°—50°N,60°—120°E)与北太平洋上空(15°—50°N,180°—120°W)温度场反相变化的现象。亚洲-太平洋涛动指数由对流层上层(500—200 hPa)温度定义,反映了亚洲-太平洋纬向热力差异。基于一个全球海-气耦合模式FGOALS_gl的20世纪气候模拟试验结果,讨论了其对20世纪亚洲-太平洋涛动指数变化的模拟能力。结果表明,较之ERA-40再分析资料(1960—1999年),模式很好地刻画出上层温度场的平均态和主导模态的空间型。从趋势上看,模式对北太平洋上空温度的年代际变化和趋势模拟较好,但未能模拟出亚洲东部陆地上空的降温趋势。从频谱分析结果看,模拟的亚洲-太平洋涛动指数2—3,a的年际变率与再分析资料相当,5-7 a周期的变率较弱。模式能够较好地模拟出与亚洲-太平洋涛动指数相关的亚洲季风区气候异常。在20世纪模拟中,外强迫因子会改变耦合系统的年际变率,在自然因子强迫下亚洲-太平洋涛动指数的功率谱向低频方向增强,人为强迫因子的作用则相反。自然强迫因子和人为强迫因子在不同时期对亚洲-太平洋涛动年际和年代际变率的作用不同。在年际变率中人为强迫因子能够控制亚洲-太平洋涛动的变率使其不致过大;在年代际变率中人为强迫因子会增强自然强迫下亚洲-太平洋涛动的变率。模式上层温度的主导模态受ENSO调制,可能影响亚洲-太平洋涛动的年际变率。因此,模式对ENSO模拟能力的缺陷是制约模式对流层上层温度及亚洲-太平洋涛动指数变率的重要因素。     

The influence of volcanic, solar and CO2 forcing on the temperatures in the Dalton Minimum (1790–1830): a model study

The Dalton Minimum (1790–1830) was a period with reduced solar irradiance and strong volcanic eruptions. Additionally, the atmospheric CO₂ concentrations started to rise from the background level of previous centuries. In this period most empirical climate reconstructions indicate a minimum in global or hemispheric temperatures. Here, we analyse several simulations starting in 1755 with the coupled atmosphere-ocean model ECHO-G driven by different forcing combinations to investigate which external forcing could have contributed most strongly to the reduced temperatures during the Dalton Minimum. Results indicate that on global and hemispheric scales, the volcanic forcing is largely responsible for the temperature drop in this period, especially during its second half, whereas changes in solar forcing and the increasing atmospheric CO₂ concentrations were of minor importance. At regional scales, especially the extratropical, the impact of volcanic forcing is much less discernible due to the large regional variability, a finding that agrees with empirical temperature reconstructions.     

The role of the individual air-sea flux components in CO2-induced changes of the ocean's circulation and climate

 In this study we investigate the role of heat, freshwater and momentum fluxes in changing the oceanic climate and thermohaline circulation as a consequence of increasing atmospheric CO₂ concentration. Two baseline integrations with a fully coupled ocean atmosphere general circulation model with either fixed or increasing atmospheric CO₂ concentrations have been performed. In a set of sensitivity experiments either freshwater (precipitation, evaporation and runoff from the continents) and/or momentum fluxes were no longer simulated, but prescribed according to one of the fully coupled baseline experiments. This approach gives a direct estimate of the contribution from the individual flux components. The direct effect of surface warming and the associated feedbacks in ocean circulation are the dominant processes in weakening the Atlantic thermohaline circulation in our model. The relative contribution of momentum and freshwater fluxes to the total response turned out to be less than 25%, each. Changes in atmospheric water vapour transport lead to enhanced freshwater input into middle and high latitudes, which weakens the overturning. A stronger export of freshwater from the Atlantic drainage basin to the Indian and Pacific ocean, on the other hand, intensifies the Atlantic overturning circulation. In total the modified freshwater fluxes slightly weaken the Atlantic thermohaline circulation. The contribution of the modified momentum fluxes has a similar magnitude, but enhances the formation of North Atlantic deep water. Salinity anomalies in the Atlantic as a consequence of greenhouse warming stem in almost equal parts from changes in net freshwater fluxes and from changes in ocean circulation caused by the surface warming due to atmospheric heat fluxes. Important effects of the momentum fluxes are a poleward shift of the front between Northern Hemisphere subtropical and subpolar gyres and a southward shift in the position of the Antarctic circumpolar current, with a clear signal in sea level. Received: 3 May 1999 / Accepted: 11 December 1999     

ENSO dynamics and seasonal cycle in the tropical Pacific as simulated by the ECHAM4/OPYC3 coupled general circulation model

 The new version of the atmospheric general circulation model (AGCM), ECHAM4, at the Max Planck Institute for Meteorology, Hamburg, has been coupled to the OPYC3 isopycnic global ocean general circulation and sea ice model in a multi-century present-day climate simulation. Non-seasonal constant flux adjustment for heat and freshwater was employed to ensure a long-term annual mean state close to present-day climatology. This study examines the simulated upper ocean seasonal cycle and interannual variability in the tropical Pacific for the first 100 years. The coupled model’s seasonal cycle of tropical Pacific SSTs is satisfactory with respect to both the warm pool variation and the Central and Eastern Pacific, with significant errors only in the cold tongue around April. The cold phase cold tongue extent and strength is as observed, and for this the heat flux adjustment does not play a decisive role. A well-established South Pacific convergence zone is characteristic for the new AGCM version. Apart from extending the southeast trades seasonal maximum to midbasin, wind stress pattern and strength are captured. Overall the subsurface structure is consistent with the observed, with a pronounced thermocline at about 150 m depth in the west and rising to the surface from 160 °W to 100 °W. The current system is better resolved than in some previous global models and, on the whole, has the expected shape. The equatorial undercurrent is correctly positioned but the core is only half as strong as observed. The north equatorial current and counter-current also have reduced maximum speeds but the April minimum is captured. As with the companion publication from Roeckner et al. this study finds pronounced tropical Eastern and Central Pacific interannual variability. Simulated and observed NINO3 sea surface temperature (SST) variability is represented by a single, rather broadband, maximum of power spectral density, centered on about 28 months for the simulation and four years for the observations. For simulation and observations, SST, windstress, and upper ocean heat content each exhibit a single dominant large-scale amplitude and phase pattern, suggesting that the model captures the essential dynamics. The amplitude of the essentially standing oscillation in SST in the NINO3 region attains the observed strength, but is weaker at the eastern boundary. Anomalies of upper ocean heat content show off-equatorial westward and equatorial eastward propagation, the latter’s arrival in the east of the basin coinciding with the SST anomalies. Equatorial wind stress anomalies near the date line provide the appropriate forcing and clearly form a response to the anomalous SST. Received: 14 June 1996 / Accepted: 11 November 1997