Weak oceanic heat transport as a cause of the instability of glacial climates

The stability of the thermohaline circulation of modern and glacial climates is compared with the help of a two dimensional ocean—atmosphere—sea ice coupled model. It turns out to be more unstable as less freshwater forcing is required to induce a polar halocline catastrophy in glacial climates. The large insulation of the ocean by the extensive sea ice cover changes the temperature boundary condition and the deepwater formation regions moves much further South. The nature of the instability is of oceanic origin, identical to that found in ocean models under mixed boundary conditions. With similar strengths of the oceanic circulation and rates of deep water formation for warm and cold climates, the loss of stability of the cold climate is due to the weak thermal stratification caused by the cooling of surface waters, the deep water temperatures being regulated by the temperature of freezing. Weaker stratification with similar overturning leads to a weakening of the meridional oceanic heat transport which is the major negative feedback stabilizing the oceanic circulation. Within the unstable regime periodic millennial oscillations occur spontaneously. The climate oscillates between a strong convective thermally driven oceanic state and a weak one driven by large salinity gradients. Both states are unstable. The atmosphere of low thermal inertia is carried along by the oceanic overturning while the variation of sea ice is out of phase with the oceanic heat content. During the abrupt warming events that punctuate the course of a millennial oscillation, sea ice variations are shown respectively to damp (amplify) the amplitude of the oceanic (atmospheric) response. This sensitivity of the oceanic circulation to a reduced concentration of greenhouse gases and to freshwater forcing adds support to the hypothesis that the millennial oscillations of the last glacial period, the so called Dansgaard—Oeschger events, may be internal instabilities of the climate system.     

Multiple sea-ice states and abrupt MOC transitions in a general circulation ocean model

Sea ice has been suggested, based on simple models, to play an important role in past glacial–interglacial oscillations via the so-called “sea-ice switch” mechanism. An important requirement for this mechanism is that multiple sea-ice extents exist under the same land ice configuration. This hypothesis of multiple sea-ice extents is tested with a state-of-the-art ocean general circulation model coupled to an atmospheric energy–moisture-balance model. The model includes a dynamic-thermodynamic sea-ice module, has a realistic ocean configuration and bathymetry, and is forced by annual mean forcing. Several runs with two different land ice distributions represent present-day and cold-climate conditions. In each case the ocean model is initiated with both ice-free and fully ice-covered states. We find that the present-day runs converge approximately to the same sea-ice state for the northern hemisphere while for the southern hemisphere a difference in sea-ice extent of about three degrees in latitude between the different runs is observed. The cold climate runs lead to meridional sea-ice extents that are different by up to four degrees in latitude in both hemispheres. While approaching the final states, the model exhibits abrupt transitions from extended sea-ice states and weak meridional overturning circulation, to less extended sea ice and stronger meridional overturning circulation, and vice versa. These transitions are linked to temperature changes in the North Atlantic high-latitude deep water. Such abrupt changes may be associated with Dansgaard–Oeschger events, as proposed by previous studies. Although multiple sea ice states have been observed, the difference between these states is not large enough to provide a strong support for the sea-ice-switch mechanism.     

Simulated decadal oscillations of the Atlantic meridional overturning circulation in a cold climate state

The significance of the Atlantic meridional overturning circulation (MOC) for regional and hemispheric climate change requires a complete understanding using fully coupled climate models. Here we present a persistent, decadal oscillation in a coupled atmosphere–ocean general circulation model. While the present study is limited by the lack of comparisons with paleo-proxy records, the purpose is to reveal a new theoretically interesting solution found in the fully-coupled climate model. The model exhibits two multi-century-long stable states with one dominated by decadal MOC oscillations. The oscillations involve an interaction between anomalous advective transport of salt and surface density in the North Atlantic subpolar gyre. Their time scale is fundamentally determined by the advection. In addition, there is a link between the MOC oscillations and North Atlantic Oscillation (NAO)-like sea level pressure anomalies. The analysis suggests an interaction between the NAO and an anomalous subpolar gyre circulation in which sea ice near and south of the Labrador Sea plays an important role in generating a large local thermal anomaly and a meridional temperature gradient. The latter induces a positive feedback via synoptic eddy activity in the atmosphere. In addition, the oscillation only appears when the Nordic Sea is completely covered by sea ice in winter, and deep convection is active only near the Irminger Sea. Such conditions are provided by a substantially colder North Atlantic climate than today.     

A global hybrid coupled model based on atmosphere-SST feedbacks

A global hybrid coupled model is developed, with the aim of studying the effects of ocean-atmosphere feedbacks on the stability of the Atlantic meridional overturning circulation. The model includes a global ocean general circulation model and a statistical atmosphere model. The statistical atmosphere model is based on linear regressions of data from a fully coupled climate model on sea surface temperature both locally and hemispherically averaged, being the footprint of Atlantic meridional overturning variability. It provides dynamic boundary conditions to the ocean model for heat, freshwater and wind-stress. A basic but consistent representation of ocean-atmosphere feedbacks is captured in the hybrid coupled model and it is more than 10 times faster than the fully coupled climate model. The hybrid coupled model reaches a steady state with a climate close to the one of the fully coupled climate model, and the two models also have a similar response (collapse) of the Atlantic meridional overturning circulation to a freshwater hosing applied in the northern North Atlantic.     

Simulation of the Atlantic meridional overturning circulation in an atmosphere-ocean global coupled model. Part II: weakening in a climate change experiment: a feedback mechanism

Most state-of-the art global coupled models simulate a weakening of the Atlantic meridional overturning circulation (MOC) in climate change scenarios but the mechanisms leading to this weakening are still being debated. The third version of the CNRM (Centre National de Recherches Météorologiques) global atmosphere-ocean-sea ice coupled model (CNRM-CM3) was used to conduct climate change experiments for the Intergovernmental Panel on Climate Change Fourth Assessment Report (IPCC AR4). The analysis of the A1B scenario experiment shows that global warming leads to a slowdown of North Atlantic deep ocean convection and thermohaline circulation south of Iceland. This slowdown is triggered by a freshening of the Arctic Ocean and an increase in freshwater outflow through Fram Strait. Sea ice melting in the Barents Sea induces a local amplification of the surface warming, which enhances the cyclonic atmospheric circulation around Spitzberg. This anti-clockwise circulation forces an increase in Fram Strait outflow and a simultaneous increase in ocean transport of warm waters toward the Barents Sea, favouring further sea ice melting and surface warming in the Barents Sea. Additionally, the retreat of sea ice allows more deep water formation north of Iceland and the thermohaline circulation strengthens there. The transport of warm and saline waters toward the Barents Sea is further enhanced, which constitutes a second positive feedback.     

The ocean’s role in setting the mean position of the Inter-Tropical Convergence Zone

Through study of observations and coupled climate simulations, it is argued that the mean position of the Inter-Tropical Convergence Zone (ITCZ) north of the equator is a consequence of a northwards heat transport across the equator by ocean circulation. Observations suggest that the hemispheric net radiative forcing of climate at the top of the atmosphere is almost perfectly symmetric about the equator, and so the total (atmosphere plus ocean) heat transport across the equator is small (order 0.2 PW northwards). Due to the Atlantic ocean’s meridional overturning circulation, however, the ocean carries significantly more heat northwards across the equator (order 0.4 PW) than does the coupled system. There are two primary consequences. First, atmospheric heat transport is southwards across the equator to compensate (0.2 PW southwards), resulting in the ITCZ being displaced north of the equator. Second, the atmosphere, and indeed the ocean, is slightly warmer (by perhaps 2 °C) in the northern hemisphere than in the southern hemisphere. This leads to the northern hemisphere emitting slightly more outgoing longwave radiation than the southern hemisphere by virtue of its relative warmth, supporting the small northward heat transport by the coupled system across the equator. To conclude, the coupled nature of the problem is illustrated through study of atmosphere–ocean–ice simulations in the idealized setting of an aquaplanet, resolving the key processes at work.     

A look at the ocean in the EC-Earth climate model

EC-Earth is a newly developed global climate system model. Its core components are the Integrated Forecast System (IFS) of the European Centre for Medium Range Weather Forecasts (ECMWF) as the atmosphere component and the Nucleus for European Modelling of the Ocean (NEMO) developed by Institute Pierre Simon Laplace (IPSL) as the ocean component. Both components are used with a horizontal resolution of roughly one degree. In this paper we describe the performance of NEMO in the coupled system by comparing model output with ocean observations. We concentrate on the surface ocean and mass transports. It appears that in general the model has a cold and fresh bias, but a much too warm Southern Ocean. While sea ice concentration and extent have realistic values, the ice tends to be too thick along the Siberian coast. Transports through important straits have realistic values, but generally are at the lower end of the range of observational estimates. Exceptions are very narrow straits (Gibraltar, Bering) which are too wide due to the limited resolution. Consequently the modelled transports through them are too high. The strength of the Atlantic meridional overturning circulation is also at the lower end of observational estimates. The interannual variability of key variables and correlations between them are realistic in size and pattern. This is especially true for the variability of surface temperature in the tropical Pacific (El Ni?o). Overall the ocean component of EC-Earth performs well and helps making EC-Earth a reliable climate model.     

海冰模式CICE4.0与LASG/IAP气候系统模式的耦合试验

利用美国Los Alamos国家实验室发展的最新海冰模式(CICE4.0)替代了LASG/IAP气候系统模式(FGOALS_g1.1)中的海冰模式(CSIM4), 形成新的耦合模式。在此基础上, 利用新的耦合模式对20世纪中后期的全球气候进行了模拟, 来检验CICE4.0对耦合模式中海冰和海洋模拟结果的改进。结果表明CICE4.0对于FGOALS_g1.1的极地气候模拟有一定改进作用, 主要表现在:(1) 南北极海冰边缘碎冰区显著减少; (2) 南大洋海表温度和海冰的模拟明显改善, 分布特征与观测非常吻合。但是新耦合模式也存在如下不足: (1) 北大西洋海冰相对偏多, 北大西洋经圈翻转环流大大减弱, 这主要是由于北大西洋海表面温度的冷误差造成的; (2) 南北极大气环流场的模拟无明显改善。此外, 本文还比较了采用不同短波辐射方案对于耦合模拟结果的影响, 结果表明, 相对于CCSM3短波辐射方案, Delta-Eddington方案模拟的海表面温度偏冷, 海冰厚度偏厚, 北大西洋经圈翻转环流略有偏弱。     

大气环流模式（SAMIL）海气耦合前后性能的比较

基于耦合器框架，中国科学院大气物理研究所大气科学和地球流体力学数值模拟国家重点实验室大气环流谱模式 （SAMIL）最近成功地实现了与海洋、海冰等气候分量模式的耦合，形成了“非通量调整”的海－陆－气－冰直接耦合的气候模式系统（FGOALS-s）。在耦合系统中，由于海温、海冰等的分布由预报模式驱动，大气与海洋、海冰之间引入了相互作用过程，这样大气环流的模拟特征与耦合前会有不同。为分析耦合系统的性能，作者对耦合前后的模拟结果进行了分析比较，重点是大气模拟特征的差异。结果表明，耦合前、后大气环流的基本特征相似，都能成功地模拟出主要的环流系统分布及季节变化，但是由于海温和海冰的模拟存在系统性的偏差，使得耦合后的大气环流受到明显影响。例如耦合后热带海温偏冷，南大洋、北太平洋和北大西洋等中纬度地区的海温偏高，导致海温等值线向高纬海域的伸展较弱，海温经向梯度减小。耦合后海冰在北极区域范围偏大，在南极周边地区则偏小。海温、海冰分布模拟的偏差影响到中、高纬低层大气的温度。热带海温偏低，使得赤道地区降水偏弱，凝结潜热减少，热带对流层中高层温度比耦合前要低，大气温度的经向梯度减小。经向温度梯度的改变，直接影响到对平均经圈环流及西风急流强度的模拟。尽管耦合系统中海温、海冰的模拟存在偏差，但在亚洲季风区，耦合后季风环流及降水等的分布都比耦合前单独大气模式的结果合理，表明通过海[CD*2]气相互作用可减少耦合前季风区的模拟误差，改善季风模拟效果。比较发现，海温、海冰模拟的偏差，除与海洋模式中经向热输送偏弱、海冰模式中海冰处理等有关外，也与大气模式中总云量模拟偏低有关。大气模式本身的误差，特别是云、辐射过程带来的误差，对耦合结果具有极为重要的影响。完全耦合后，这些误差通过与海洋、海冰的反馈作用而放大。因此，对于FGOALS-s而言，要提高耦合系统的整体性能，除改进各气候分量模式的模拟性能外，需要重点改进大气模式中的云、辐射过程。     

Key features of the IPSL ocean atmosphere model and its sensitivity to atmospheric resolution

This paper presents the major characteristics of the Institut Pierre Simon Laplace (IPSL) coupled ocean–atmosphere general circulation model. The model components and the coupling methodology are described, as well as the main characteristics of the climatology and interannual variability. The model results of the standard version used for IPCC climate projections, and for intercomparison projects like the Paleoclimate Modeling Intercomparison Project (PMIP 2) are compared to those with a higher resolution in the atmosphere. A focus on the North Atlantic and on the tropics is used to address the impact of the atmosphere resolution on processes and feedbacks. In the North Atlantic, the resolution change leads to an improved representation of the storm-tracks and the North Atlantic oscillation. The better representation of the wind structure increases the northward salt transports, the deep-water formation and the Atlantic meridional overturning circulation. In the tropics, the ocean–atmosphere dynamical coupling, or Bjerknes feedback, improves with the resolution. The amplitude of ENSO (El Niño-Southern oscillation) consequently increases, as the damping processes are left unchanged.     

Impact of freshwater discharge from the Greenland ice sheet on North Atlantic climate variability

Using a coupled ocean–atmosphere general circulation model, we investigated the impact of Greenland ice sheet melting on North Atlantic climate variability. The positive-degree day (PDD) method was incorporated into the model to control continental ice melting (PDD run). Models with and without the PDD method produce a realistic pattern of North Atlantic sea surface temperature (SST) variability that fluctuates from decadal to multidecadal periods. However, the interdecadal variability in PDD run is significantly dominated in the longer time scale compared to that in the run without PDD method. The main oscillatory feature in these experiments likely resembles the density-driven oscillatory mode. A reduction in the ocean density over the subpolar Atlantic results in suppression of the Atlantic Meridional Overturning Circulation (AMOC), leading to a cold SST due to a weakening of northward heat transport. The decreased surface evaporation associated with the cold SST further reduces the ocean density and thus, simultaneously acts as a positive feedback mechanism. The southward meridional current associated with the suppressed AMOC causes a positive tendency in the ocean density through density advection, which accounts for the phase transition of this oscillatory mode. The Greenland ice melting process reduces the mean meridional current and meridional density gradient because of additional fresh water flux, which suppress the delayed negative feedback due to meridional density advection. As a result, the oscillation period becomes longer and the transition is more delayed.     

Oceanic Origin of A Recent La Nina-Like Trend in the Tropical Pacific

Global ocean temperature has been rising since the late 1970s at a speed unprecedented during the past century of recordkeeping.This accelerated warming has profound impacts not only on the marine ecosystem and oceanic carbon uptake but also on the global water cycle and climate.During this rapid warming period,the tropical Pacific displays a pronounced La Nin a-like trend,characterized by an intensification of west-east SST gradient and of atmospheric zonal overturning circulation,namely the Walker circulation.This La Nin a-like trend differs from the El Nin o-like trend in warm climate projected by most climate models,and cannot be explained by responses of the global water cycle to warm climate.The results of this study indicate that the intensification of the zonal SST gradient and the Walker circulation are associated with recent strengthening of the upper-ocean meridional overturning circulation.     

Climate sensitivity to changes in solar insolation in a simple coupled climate model

A simple coupled ocean, atmosphere and sea-ice model is presented. The idealised model consists of a zonally averaged land and ocean strip of constant angular width extending from pole to pole. The meridional energy transport in the ocean is modelled by contributions from the large scale thermohaline overturning cells and from horizontal diffusive fluxes. The atmospheric meridional energy transports are parametrised as diffusive fluxes in addition to advective transports in the Hadley domain. This parametrisation resolves the equatorward moisture transport as well as the poleward transport of potential energy in the upper branch of the Hadley circulation. The model reproduces the annual averaged meridional energy transports in the climate system with a small number of free model parameters. The basic feedbacks between the three climatic components are studied by investigating the model's sensitivity towards reductions in the solar insolation. It is found that the meridional energy transport in the ocean does not amplify the ice albedo feedback. This has important implications for modelling the climate sensitivity in atmosphere-only models, as these would exaggerate the sensitivity to changes in the solar insolation if their parametrisations of the meridional energy transport are constrained by surface temperatures. The role of the dependence of the atmospheric transports on the meridional temperature gradient is shown to have a significant influence on the sensitivity on the coupled model, and the inclusion of seasonal cycles greatly increase the models sensitivity. The Hadley circulation does significantly alter the strength of the ice-albedo feedback in the coupled model. The idealised configuration of the model makes it a useful tool for studying the feedbacks in the ocean-atmosphere-sea ice system in the context of the "Snowball Earth" hypothesis.     

NAO–ocean circulation interactions in a coupled general circulation model

The interplay between the North Atlantic Oscillation (NAO) and the large scale ocean circulation is inspected in a twentieth century simulation conducted with a state-of-the-art coupled general circulation model. Significant lead–lag covariance between oceanic and tropospheric variables suggests that the system supports a damped oscillatory mode involving an active ocean–atmosphere coupling, with a typical NAO-like space structure and a 5 years timescale, qualitatively consistent with a mid-latitude delayed oscillator paradigm. The two essential processes governing the oscillation are (1) a negative feedback between ocean gyre circulation and the high latitude SST meridional gradient and (2) a positive feedback between SST and the NAO. The atmospheric NAO pattern appears to have a weaker projection on the ocean meridional overturning, compared to the gyre circulation, which leads to a secondary role for the thermohaline circulation in driving the meridional heat transport, and thus the oscillatory mode.     

Simulation and Improvements of Oceanic Circulation and Sea Ice by the Coupled Climate System Model FGOALS-f3-L

This study documents simulated oceanic circulations and sea ice by the coupled climate system model FGOALS-f3-L developed at the State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics, Institute of Atmospheric Physics, Chinese Academy of Sciences, under historical forcing from phase 6 of the Coupled Model Intercomparison Project (CMIP6). FGOALS-f3-L reproduces the fundamental features of global oceanic circulations, such as sea surface temperature (SST), sea surface salinity (SSS), mixed layer depth (MLD), vertical temperature and salinity, and meridional overturning circulations. There are notable improvements compared with the previous version, FGOALS-s2, such as a reduction in warm SST biases near the western and eastern boundaries of oceans and salty SSS biases in the tropical western Atlantic and eastern boundaries, and a mitigation of deep MLD biases at high latitudes. However, several obvious biases remain. The most significant biases include cold SST biases in the northwestern Pacific (over 4°C), freshwater SSS biases and deep MLD biases in the subtropics, and temperature and salinity biases in deep ocean at high latitudes. The simulated sea ice shows a reasonable distribution but stronger seasonal cycle than observed. The spatial patterns of sea ice are more realistic in FGOALS-f3-L than its previous version because the latitude–longitude grid is replaced with a tripolar grid in the ocean and sea ice model. The most significant biases are the overestimated sea ice and underestimated SSS in the Labrador Sea and Barents Sea, which are related to the shallower MLD and weaker vertical mixing.     

Impact of a realistic river routing in coupled ocean–atmosphere simulations of the Last Glacial Maximum climate

The presence of large ice sheets over North America and North Europe at the Last Glacial Maximum (LGM) strongly impacted Northern hemisphere river pathways. Despite the fact that such changes may significantly alter the freshwater input to the ocean, modified surface hydrology has never been accounted for in coupled ocean–atmosphere general circulation model simulations of the LGM climate. To reconstruct the LGM river routing, we use the ICE-5G LGM topography. Because of the uncertainties in the extent of the Fennoscandian ice sheet in the Eastern part of the Kara Sea, we consider two more realistic river routing scenarios. The first scenario is characterised by the presence of an ice dammed lake south of the Fennoscandian ice sheet, and corresponds to the ICE-5G topography. This lake is fed by the Ob and Yenisei rivers. In the second scenario, both these rivers flow directly into the Arctic Ocean, which is more consistent with the latest QUEEN ice sheet margin reconstructions. We study the impact of these changes on the LGM climate as simulated by the IPSL_CM4 model and focus on the overturning thermohaline circulation. A comparison with a classical LGM simulation performed using the same model and modern river basins as designed in the PMIP2 exercise leads to the following conclusions: (1) The discharge into the North Atlantic Ocean is increased by 2,000 m³/s between 38° and 54°N in both simulations that contain LGM river routing, compared to the classical LGM experiment. (2) The ice dammed lake is shown to have a weak impact, relative to the classical simulation, both in terms of climate and ocean circulation. (3) In contrast, the North Atlantic deep convection and meridional overturning are weaker than during the classical LGM run if the Ob and Yenisei rivers flow directly into the Arctic Ocean. The total discharge into the Arctic Ocean is increased by 31,000 m³/s, relative to the classical LGM simulation. Consequentially, northward ocean heat transport is weaker, and sea ice more extensive, in better agreement with existing proxy data.     

Arctic climate change in 21st century CMIP5 simulations with EC-Earth

The Arctic climate change is analyzed in an ensemble of future projection simulations performed with the global coupled climate model EC-Earth2.3. EC-Earth simulates the twentieth century Arctic climate relatively well but the Arctic is about 2 K too cold and the sea ice thickness and extent are overestimated. In the twenty-first century, the results show a continuation and strengthening of the Arctic trends observed over the recent decades, which leads to a dramatically changed Arctic climate, especially in the high emission scenario RCP8.5. The annually averaged Arctic mean near-surface temperature increases by 12 K in RCP8.5, with largest warming in the Barents Sea region. The warming is most pronounced in winter and autumn and in the lower atmosphere. The Arctic winter temperature inversion is reduced in all scenarios and disappears in RCP8.5. The Arctic becomes ice free in September in all RCP8.5 simulations after a rapid reduction event without recovery around year 2060. Taking into account the overestimation of ice in the twentieth century, our model results indicate a likely ice-free Arctic in September around 2040. Sea ice reductions are most pronounced in the Barents Sea in all RCPs, which lead to the most dramatic changes in this region. Here, surface heat fluxes are strongly enhanced and the cloudiness is substantially decreased. The meridional heat flux into the Arctic is reduced in the atmosphere but increases in the ocean. This oceanic increase is dominated by an enhanced heat flux into the Barents Sea, which strongly contributes to the large sea ice reduction and surface-air warming in this region. Increased precipitation and river runoff lead to more freshwater input into the Arctic Ocean. However, most of the additional freshwater is stored in the Arctic Ocean while the total Arctic freshwater export only slightly increases.     

加利福尼亚附近的海温强迫及其与北太平洋年代际振荡的可能联系

本文利用1958～2018年期间海表面温度（SST）异常和湍流热通量异常变化的关系,探讨了其与北太平洋年代际振荡（PDO）相关的年际和年代际时间尺度上在不同海域的海气相互作用特征。结果表明:在年际尺度上,黑潮—亲潮延伸区（KOE）表现为显著大气强迫海洋,赤道中东太平洋表现为显著海洋强迫大气;在年代际尺度上,PDO北中心表现为大气强迫海洋,加利福尼亚附近则表现为显著海洋强迫大气。进一步分析表明:加利福尼亚附近区域是北太平洋准12年振荡的关键区域之一,与PDO准十年的周期类似,加利福尼亚附近的冷（暖）海温对应其上有反气旋（气旋）型环流,赤道中太平洋海水上翻和北太平洋东部副热带区域经向风应力的变化是北太平洋准12年振荡的另外两个重要环节。     

Coupled atmosphere–ocean data assimilation experiments with a low-order climate model

A simple idealized atmosphere–ocean climate model and an ensemble Kalman filter are used to explore different coupled ensemble data assimilation strategies. The model is a low-dimensional analogue of the North Atlantic climate system, involving interactions between large-scale atmospheric circulation and ocean states driven by the variability of the Atlantic meridional overturning circulation (MOC). Initialization of the MOC is assessed in a range of experiments, from the simplest configuration consisting of forcing the ocean with a known atmosphere to performing fully coupled ensemble data assimilation. “Daily” assimilation (that is, at the temporal frequency of the atmospheric observations) is contrasted with less frequent assimilation of time-averaged observations. Performance is also evaluated under scenarios in which ocean observations are limited to the upper ocean or are non-existent. Results show that forcing the idealized ocean model with atmospheric analyses is inefficient at recovering the slowly evolving MOC. On the other hand, daily assimilation rapidly leads to accurate MOC analyses, provided a comprehensive set of oceanic observations is available for assimilation. In the absence of sufficient observations in the ocean, the assimilation of time-averaged atmospheric observations proves to be more effective for MOC initialization, including the case where only atmospheric observations are available.