温盐环流与全球增暖的数值模拟 (一)纬向平均温盐环流的模拟

本文是用简单海一气耦合模型模拟温盐环流在全球增暖事件中作用的研究工作的第一部分。为了建立一个简单海一气耦合模型，我们首先根据Ｗｒｉｇｈｔ和Ｓｔｏｋｅｒ等人的设计复制出一个包括大西洋、太平洋和南大洋在内的二维温盐环流模式，从等温、等盐和无运动的初始状态出发，在给定的年平均海表强迫下将模式积分了４０００年，模拟出了和原作相似的温盐环流。对模拟结果的分析表明，相对于北太平洋而言，北大西洋北部的高盐、低温特点（后者是由两大洋在地理上的差别决定的）是形成当代温盐环流的主要原因；从与温盐环流相联系的海表热通量来看，北大西洋北部是向大气提供热量的主要源地；模式温盐环流对于海表盐度通量的敏感性试验的结果表明，对于纬圈平均的二维模式而言，要想模拟出合理的温盐环流就必须人为地提高北大西洋北部的海表盐度，文章分析了这种作法的物理根据；模式中的对流过程对于温盐环流的维持是至关重要的，对比有无季节循环的试验结果可以看出，虽然温度场的明显的季节变化只出现在模式的最上面两层，但由于引进季节循环后冬季高纬海洋的对流活动加强，后者直接影响到温盐环流，使更多的深海热量上传并向大气释放。这是使海洋温跃层得以保持合理．厚度的一个重要原因。     

Simulation of the atmospheric circulation in the Weddell Sea region using the limited-area model REMO

Summary ?To analyse the applicability of a limited-area atmosphere model to the Southern Ocean, a one-year simulation for 1985 is performed using the REgional MOdel REMO at 55-km horizontal grid-spacing implemented for the Antarctic regions of the Weddell, Bellingshausen and Amundsen Seas. To evaluate the performance of REMO, a comparison of model results to observations and to reanalysis/analysis data sets is carried out. REMO is initialized and driven at the lateral and lower boundaries by data of the European Centre for Medium-Range Weather Forecasts (ECMWF) reanalysis (ERA15). Overall, REMO is an appropriate tool for further climate studies in Antarctic regions. It reproduces reasonably well basic spatial patterns and the seasonal cycle of the atmospheric circulation. However, the simulated mean sea level pressure (MSLP) is predominantly lower than the MSLP provided by observations and by ERA. Considerable temperature differences in the lower troposphere over sea ice in winter cause discrepancies between the REMO and ERA pressure fields in the mid-troposphere too. The precipitation rate P of the REMO simulation agrees qualitatively well with main features of the observed climatological spatial distribution described in literature. The seasonal cycle of P in the inner Weddell Sea reflects the Antarctic semi-annual oscillation. Concerning the forcing fields, the ERA sea ice surface temperatures in winter are generally higher than satellite derived surface temperatures. Although the differences are 10 to 15 K in the southern Weddell Sea, this deficiency of the ERA data hardly influences the mean large-scale circulation. Received October 10, 2001; revised April 22, 2002; accepted May 12, 2002     

ARE THERE INTERANNUAL-TO-DECADAL SCALE OSCILLATIONS ASSOCIATED WITH SEA ICE-THERMOHALINE CIRCULATION INTERACTIONS IN A SIMPLE COUPLED ATMO-SPHERE-OCEAN-SEA ICE MODEL?

1　INTRODUCTIONIn order to gain further insight into the nature of decadal- scale climate variability at highlatitudes( e.g.,Mysak et al.,1 990 ;Deser and black- mon,1 993) ,there have been a number ofrecent model studies of sea ice- thermohaline circulation interactions which exhibitoscillationson this timescale( Yang and Neelin,1 993;Zhang et al.,1 995 ;Yang and Huang,1 996 ) .Acommon feature of these studies is that the ocean models are integrated using mixedboundary conditions( MBC…     

Simulation of the asian monsoon by IAP AGCM coupled with an advanced land surface model (IAP94)

ＳｉｍｕｌａｔｉｏｎｏｆｔｈｅＡｓｉａｎＭｏｎｓｏｏｎｂｙＩＡＰＡＧＣＭＣｏｕｐｌｅｄｗｉｔｈａｎＡｄｖａｎｃｅｄＬａｎｄＳｕｒｆａｃｅＭｏｄｅｌ（ＩＡＰ９４）ＺｅｎｇＱｉｎｇｃｕｎ（曾庆存）,ＤａｉＹｏｎｇｊｉｕ（戴永久）ａｎｄＸｕｅＦｅｎｇ（薛峰...     

A statistical�Cdynamical scheme for reconstructing ocean forcing in the Atlantic. Part II: methodology, validation and application to high-resolution ocean models

A novel statistical?Cdynamical scheme has been developed to reconstruct the sea surface atmospheric variables necessary to force an ocean model. Multiple linear regressions are first built over a so-called learning period and over the entire Atlantic basin from the observed relationship between the surface wind conditions, or predictands, and the anomalous large scale atmospheric circulations, or predictors. The latter are estimated in the extratropics by 500?hPa geopotential height weather regimes and in the tropics by low-level wind classes. The transfer function further combined to an analog step is then used to reconstruct all the surface variables fields over 1958?C2002. We show that the proposed hybrid scheme is very skillful in reproducing the mean state, the seasonal cycle and the temporal evolution of all the surface ocean variables at interannual timescale. Deficiencies are found in the level of variance especially in the tropics. It is underestimated for 2-m temperature and humidity as well as for surface radiative fluxes in the interannual frequency band while it is slightly overestimated at higher frequency. Decomposition in empirical orthogonal function (EOF) shows that the spatial and temporal coherence of the forcing fields is however very well captured by the reconstruction method. For dynamical downscaling purposes, reconstructed fields are then interpolated and used to carry out a high-resolution oceanic simulation using the NATL4 (1/4°) model integrated over 1979?C2001. This simulation is compared to a reference experiment where the original observed forcing fields are prescribed instead. Mean states between the two experiments are virtually undistinguishable both in terms of surface fluxes and ocean dynamics estimated by the barotropic and the meridional overturning streamfunctions. The 3-dimensional variance of the simulated ocean is well preserved at interannual timescale both for temperature and salinity except in the tropics where it is underestimated. The main modes of interannual variability assessed through EOF are correctly reproduced for sea surface temperature, barotropic streamfunction and mixed layer depth both in terms of spatial structure and temporal evolution. Collectively, our results provide evidence that the statistical?Cdynamical scheme presented in this two-part study is an efficient and promising tool to infer oceanic changes (in particular those related to the wind-driven circulation) due to modifications in the large-scale atmospheric circulation. As a prerequisite, we have here validated the method for present-day climate; we encourage its use for climate change studies with some adaptations though.     

How Well Do CMIP6 and CMIP5 Models Simulate the Climatological Seasonal Variations in Ocean Salinity?

This paper includes a comprehensive assessment of 40 models from the Coupled Model Intercomparison Project phase 5 (CMIP5) and 33 models from the CMIP phase 6 (CMIP6) to determine the climatological and seasonal variation of ocean salinity from the surface to 2000 m. The general pattern of the ocean salinity climatology can be simulated by both the CMIP5 and CMIP6 models from the surface to 2000-m depth. However, this study shows an increased fresh bias in the surface and subsurface salinity in the CMIP6 multimodel mean, with a global average of ?0.44 g kg^?1 for the sea surface salinity (SSS) and ?0.26 g kg^?1 for the 0–1000-m averaged salinity (S1000) compared with the CMIP5 multimodel mean (?0.25 g kg^?1 for the SSS and ?0.07 g kg^?1 for the S1000). In terms of the seasonal variation, both CMIP6 and CMIP5 models show positive (negative) anomalies in the first (second) half of the year in the global average SSS and S1000. The model-simulated variation in SSS is consistent with the observations, but not for S1000, suggesting a substantial uncertainty in simulating and understanding the seasonal variation in subsurface salinity. The CMIP5 and CMIP6 models overestimate the magnitude of the seasonal variation of the SSS in the tropics in the region 20°S–20°N but underestimate the magnitude of the seasonal change in S1000 in the Atlantic and Indian oceans. These assessments show new features of the model errors in simulating ocean salinity and support further studies of the global hydrological cycle.     

大气环流的年代际变化 II.GCM数值模拟研究

类似大气环流模式比较计划（AMIP）的数值模拟,将实际观测的海表水温（SST）资料引入模式进行40多年的数值积分,得到长时间的大气环流模拟结果。分析数值模拟结果发现,无论是大气中的主要涛动和遥相关型,还是重要大气环流系统都极为清楚地存在着年代际变化特征,包括10～20年准周期振荡和可能的30年以上的准周期振荡;而且上述主要环流系统的形势及其年代际变化大都与实际观测资料所给出的结果相一致。顺便分析中国东部气候的模拟结果,年代际变化特征（包括60年代的气候突变）也很清楚,并同大气环流变化配合十分合理。结果也表明,同研究季节和年际变化一样,大气环流模式（AGCM）数值模拟也是研究大气环流和气候年代际变化的有效方法。     

CWRF对中国近海热带气旋活动季节和年际变化及环境场影响的模拟

利用欧洲中期天气预报中心的ERA-Interim再分析资料驱动CWRF模式对1982-2016年中国近海的热带气旋活动进行了模拟,分析了CWRF对热带气旋频数季节、年际变化和路径的模拟能力,并探讨了环境场模拟对热带气旋模拟的影响.结果 表明:CWRF能够合理模拟热带气旋频数的季节和年际变化,但模拟的频数较观测总体偏低,...     

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.     

A mechanism for Atlantic multidecadal variability in the Kiel Climate Model

Atlantic Multidecadal Variability (AMV) is investigated in a millennial control simulation with the Kiel Climate Model (KCM), a coupled atmosphere–ocean–sea ice model. An oscillatory mode with approximately 60 years period and characteristics similar to observations is identified with the aid of three-dimensional temperature and salinity joint empirical orthogonal function analysis. The mode explains 30 % of variability on centennial and shorter timescales in the upper 2,000 m of the North Atlantic. It is associated with changes in the Atlantic Meridional Overturning Circulation (AMOC) of ±1–2 Sv and Atlantic Sea Surface Temperature (SST) of ±0.2 °C. AMV in KCM results from an out-of-phase interaction between horizontal and vertical ocean circulation, coupled through Irminger Sea convection. Wintertime convection in this region is mainly controlled by salinity anomalies transported by the Subpolar Gyre (SPG). Increased (decreased) dense water formation in this region leads to a stronger (weaker) AMOC after 15 years, and this in turn leads to a weaker (stronger) SPG after another 15 years. The key role of salinity variations in the subpolar North Atlantic for AMV is confirmed in a 1,000 year long simulation with salinity restored to model climatology: No low frequency variations in convection are simulated, and the 60 year mode of variability is absent.     

大气季节内振荡的耦合模式数值模拟

分析GOALS/LASG海气耦合模式10年积分200hPa纬向风场的逐日输出结果,引用1980～1989年期间逐日的NCEP/NCAR再分析资料作为实测对照,结果显示该耦合模式抓住了热带大气低频振荡(IO)的基本时空分布特征,模拟IO的强度较多数大气模式强而接近真实,但空间一致性仍不清晰,典型周期不够显著。NCEP资料与耦合模式都反映模拟IO的季节变化与其年际变化有关,模拟较强IO的年份表现IO的季节变化特征也较真实。模拟IO的年际变化与热带东太平洋的SST呈明显的负相关变化。SST暖异常的年份,IO活动较弱。IO变化滞后于SST异常60天左右的相关性最显著。对比单独积分GOALS/LASG的大气模式的结果,发现二者的主要差别在于耦合模式再现IO的季节性特征更真实,反映了海气耦合对IO变化的调制作用。利用海气耦合模式,理解IO对流活动与上层海洋的相互作用过程,是真实描述IO必要的手段。     

Climatology and interannual variability simulated by the ARPEGE-OPA coupled model

A 10-year simulation with a coupled ocean-atmosphere general circulation model (CGCM) is presented. The model consists of the climate version of the Météo-France global forecasting model, ARPEGE, coupled to the LODYC oceanic model, OPA, by the CERFACS coupling package OASIS. The oceanic component is dynamically active over the tropical Pacific, while climatological time-dependent sea surface temperatures (SSTs) are prescribed outside of the Pacific domain. The coupled model shows little drift and exhibits a very regular seasonal cycle. The climatological mean state and seasonal cycle are well simulated by the coupled model. In particular, the oceanic surface current pattern is accurately depicted and the location and intensity of the Equatorial Undercurrent (EUC) are in good agreement with available data. The seasonal cycle of equatorial SSTs captures quite realistically the annual harmonic. Some deficiencies remain including a weak zonal equatorial SST gradient, underestimated wind stress over the Pacific equatorial band and an additional inter-tropical convergence zone (ITCZ) south of the equator in northern winter and spring. Weak interannual variability is present in the equatorial SST signal with a maximum amplitude of 0.5°C.     

东亚夏季风的自然变率——NCARCam3模拟结果分析

利用季节循环的全球观测海表温度及海冰驱动NCARCam3全球大气环流模式的100a模拟结果,通过定义东亚夏季风指数,分析了模拟的大气内部变化中东亚夏季风的变化特征。结果表明：模拟的东亚夏季风自然变率主要表现为3—7a较显著的年际周期,并具有较明显的年代际变化特征。在弱夏季风年代,亚洲大陆海平面气压增强,日本附近及东亚沿海地区海平面气压降低;500hPa位势高度上,欧洲地区为负高度距平,里海附近地区为正高度距平,日本及其以东太平洋为负高度距平,易形成类似欧亚（EU）型的遥相关波列。在强夏季风年代,其环流异常分布基本与弱夏季风年代相反。模拟的东亚夏季风变化与夏季大气内部500hPa高度场上EU型遥相关波列的关系密切。     

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.     

A world ocean model for greenhouse sensitivity studies: resolution intercomparison and the role of diagnostic forcing

We have developed an improved version of a world ocean model with the intention of coupling to an atmospheric model. This article documents the simulation capability of this 1° global ocean model, shows improvements over our earlier 5° version, and compares it to features simulated with a 0.5° model. These experiments use a model spin-up methodology whereby the ocean model can subsequently be coupled to an atmospheric model and used for order 100-year coupled model integrations. With present-day computers, 1° is a reasonable compromise in resolution that allows for century-long coupled experiments. The 1° ocean model is derived from a 0.5°-resolution model developed by A. Semtner (Naval Postgraduate School) and R. Chervin (National Center for Atmospheric Research) for studies of the global eddy-resolving world ocean circulation. The 0.5° bottom topography and continental outlines have been altered to be compatible with the 1° resolution, and the Arctic Ocean has been added. We describe the ocean simulation characteristics of the 1° version and compare the result of weakly constraining (three-year time scale) the three-dimensional temperature and salinity fields to the observations below the thermocline (710 m) with the model forced only at the top of the ocean by observed annual mean wind stress, temperature, and salinity. The 1° simulations indicate that major ocean circulation patterns are greatly improved compared to the 5° version and are qualitatively reproduced in comparison to the 0.5° version. Using the annual mean top forcing alone in a 100-year simulation with the 1° version preserves the general features of the major observed temperature and salinity structure with most climate drift occurring mainly beneath the thermocline in the first 50–75 years. Because the thermohaline circulation in the 1° version is relatively weak with annual mean forcing, we demonstrate the importance of the seasonal cycle by performing two sensitivity experiments. Results show a dramatic intensification of the meridional overturning circulation (order of magnitude) with perpetual winter surface temperature forcing in the North Atlantic and strong intensification (factor of three) with perpetual early winter temperatures in that region. These effects are felt throughout the Atlantic (particularly an intensified and northward-shifted Gulf Stream outflow). In the Pacific, the temperature gradient strengthens in the thermocline, thus helping counter the systematic error of a thermocline that is too diffuse.Partial support is provided by the Office of Health and Environmental Research of the US Department of Energy The National Center for Atmospheric Research is sponsored by the National Science Foundation     

A zonally averaged,three-basin ocean circulation model for climate studies

A two-dimensional, three-basin ocean model suitable for long-term climate studies is developed. The model is based on the zonally averaged form of the primitive equations written in spherical coordinates. The east-west density difference which arises upon averaging the momentum equations is taken to be proportional to the meridional density gradient. Lateral exchanges of heat and salt between the basins are explicitly resolved. Moreover, the model includes bottom topography and has representations of the Arctic Ocean and of the Weddell and Ross seas. Under realistic restoring boundary conditions, the model reproduces the global conveyor belt: deep water is formed in the Atlantic between 60 and 70°N at a rate of about 17 Sv (1 Sv=10⁶ m3 s^–1) and in the vicinity of the Antarctic continent, while the Indian and Pacific basins show broad upwelling. Superimposed on this thermohaline circulation are vigorous wind-driven cells in the upper thermocline. The simulated temperature and salinity fields and the computed meridional heat transport compare reasonably well with the observational estimates. When mixed boundary conditions (i.e., a restoring condition on sea-surface temperature and flux condition on sea-surface salinity) are applied, the model exhibits an irregular behavior before reaching a steady state characterized by self-sustained oscillations of 8.5-y period. The conveyor-belt circulation always results at this stage. A series of perturbation experiments illustrates the ability of the model to reproduce different steady-state circulations under mixed boundary conditions. Finally, the model sensitivity to various factors is examined. This sensitivity study reveals that the bottom topography and the presence of a submarine meridional ridge in the zone of the Drake Passage play a crucial role in determining the properties of the model bottom-water masses. The importance of the seasonality of the surface forcing is also stressed.     

Circulation patterns, daily precipitation in Portugal and implications for climate change simulated by the second Hadley Centre GCM

Based on principal component analysis (PCA) and a k-means clustering algorithm, daily mean sea level pressure (MSLP) fields over the northeastern Atlantic and Western Europe, simulated by the Hadley Centre's second generation coupled ocean-atmosphere GCM (HADCM2) control run (HADCM2CON), are validated by comparison with the observed daily MSLP fields. It is clear that HADCM2 reproduces daily MSLP fields and its seasonal variability over the region very well, despite suffering from some deficiencies, such as the systematic displacement of the atmospheric centres of action. Four daily circulation patterns, previously identified from the observed daily MSLP fields over the area and well related to daily precipitation in Portugal, were also well classified from the daily MSLP fields simulated by HADCM2. The model can also simulate rather successfully the relationships between the four daily circulation patterns and daily precipitation in southern Portugal. However, compared with observations, daily precipitation intensities simulated by the model are too weak in southern Portugal. Nevertheless, HADCM2 represents a considerable improvement relative to the UKTR experiment. The results described here imply that it is doubtful whether regional precipitation scenarios provided by HADCM2 can be directly applied in impact studies and that a downscaling technique, based on daily circulation patterns, might be successful in reproducing local and regional precipitation characteristics. Moreover, the four circulation patterns can also be clearly identified in the two perturbed experiments, one under greenhouse gases forcing only (HADCM2GHG) and the other under additional forcing of sulphate aerosol (HADCM2SUL), although changes in the frequencies of occurrence of certain circulation patterns are found. Nevertheless, the observed links between regional precipitation in southern Portugal and large-scale atmospheric circulation seem likely to hold in the model's perturbed climate. It is therefore credible to use those links to downscale large-scale atmospheric circulation from GCM simulations to obtain future precipitation scenarios in southern Portugal. Received: 21 August 1998 / Accepted: 28 May 1999     

Dynamical contribution to sea surface salinity variations in the eastern Gulf of Guinea based on numerical modelling

In this study, we analyse the seasonal variability of the sea surface salinity (SSS) for two coastal regions of the Gulf of Guinea from 1995 to 2006 using a high resolution model (1/12°) embedded in a Tropical Atlantic (1/4°) model. Compared with observations and climatologies, our model demonstrates a good capability to reproduce the seasonal and spatial variations of the SSS and mixed layer depth. Sensitivity experiments are carried out to assess the respective impacts of precipitations and river discharge on the spatial structure and seasonal variations of the SSS in the eastern part of the Gulf of Guinea. In the Bight of Biafra, both precipitations and river runoffs are necessary to observe permanent low SSS values but the river discharge has the strongest impact on the seasonal variations of the SSS. South of the equator, the Congo river discharge alone is sufficient to explain most of the SSS structure and its seasonal variability. However, mixed layer budgets for salinity reveal the necessity to take into account the horizontal and vertical dynamics to explain the seasonal evolution of the salinity in the mixed layer. Indeed evaporation, precipitations and runoffs represent a relatively small contribution to the budgets locally at intraseasonal to seasonal time scales. Horizontal advection always contribute to spread the low salinity coastal waters offshore and thus decrease the salinity in the eastern Gulf of Guinea. For the Bight of Biafra and the Congo plume region, the strong seasonal increase of the SSS observed from May/June to August/September, when the trade winds intensify, results from a decreasing offshore spread of freshwater associated with an intensification of the salt input from the subsurface. In the Congo plume region, the subsurface salt comes mainly from advection due to a strong upwelling but for the Bight of Biafra, entrainment and vertical mixing also play a role. The seasonal evolution of horizontal advection in the Bight of Biafra is mainly driven by eddy correlations between salinity and velocities, but it is not the case in the Congo plume.     

Modelling the sea ice-ocean seasonal cycle in Hudson Bay, Foxe Basin and Hudson Strait, Canada

The seasonal cycle of water masses and sea ice in the Hudson Bay marine system is examined using a three-dimensional coastal ice-ocean model, with 10 km horizontal resolution and realistic tidal, atmospheric, hydrologic and oceanic forcing. The model includes a level 2.5 turbulent kinetic energy equation, multi-category elastic-viscous-plastic sea-ice rheology, and two layer sea ice with a single snow layer. Results from a two-year long model simulation between August 1996 and July 1998 are analyzed and compared with various observations. The results demonstrate a consistent seasonal cycle in atmosphere-ocean exchanges and the formation and circulation of water masses and sea ice. The model reproduces the summer and winter surface mixed layers, the general cyclonic circulation including the strong coastal current in eastern Hudson Bay, and the inflow of oceanic waters into Hudson Bay. The maximum sea-ice growth rates are found in western Foxe Basin, and in a relatively large and persistent polynya in northwestern Hudson Bay. Sea-ice advection and ridging are more important than local thermodynamic growth in the regions of maximum sea-ice cover concentration and thickness that are found in eastern Foxe Basin and southern Hudson Bay. The estimate of freshwater transport to the Labrador Sea confirms a broad maximum during wintertime that is associated with the previous summers freshwater moving through Hudson Strait from southern Hudson Bay. Tidally driven mixing is shown to have a strong effect on the modeled ice-ocean circulation.