共查询到20条相似文献，搜索用时 15 毫秒
1.
《Deep Sea Research Part II: Topical Studies in Oceanography》1999,46(6-7):1385-1425
In high-latitude oceans with seasonal ice cover, the ice and the low-salinity mixed layer form an interacting barrier for the heat flux from the ocean to the atmosphere. The presence of a less dense surface layer allows ice to form, and the ice cover reduces the heat loss to the atmosphere. The ice formation weakens the stability at the base of the mixed layer, leading to stronger entrainment and larger heat flux from below. This heat transport retards, and perhaps stops, the growth of the ice cover. As much heat is then entrained from below as is lost to the atmosphere. This heat loss further reduces the stability, and unless a net ice melt occurs, the mixed layer convects. Two possibilities exist: (1) A net ice melt, sufficient to retain the stability, will always occur and convection will not take place until all ice is removed. The deep convection will then be thermal, deepening the mixed layer. (2) The ice remains until the stability at the base of the mixed layer disappears. The mixed layer then convects, through haline convection, into the deep ocean. Warm water rises towards the surface and the ice starts to melt, and a new mixed layer is reformed. The present work discusses the interactions between ice cover and entrainment during winter, when heat loss to the atmosphere is present. One crucial hypothesis is introduced: “When ice is present and the ocean loses sensible heat to the atmosphere and to ice melt, the buoyancy input at the sea surface due to ice melt is at a minimum”. Using a one-dimensional energy-balance model, applied to the artificial situation, where ice melts directly on warmer water, it is found that this corresponds to a constant fraction of the heat loss going to ice melt. It is postulated that this partitioning holds for the ice cover and the mixed layer in the high-latitude ocean. When a constant fraction of heat goes to ice melt, at least one deep convection event occurs, before the ice cover can be removed by heat entrained from below. After one or several convection events the ice normally disappears and a deep-reaching thermal convection is established. Conditions appropriate for the Weddell Sea and the Greenland Sea are examined and compared with field observations. With realistic initial conditions no convection occurs in the warm regime of the Weddell Sea. A balance between entrained heat and atmospheric heat loss is established and the ice cover remains throughout the winter. At Maud Rise convection may occur, but late in winter and normally no polynya can form before the summer ice melt. In the central Greenland Sea the mixed layer generally convects early in winter and the ice is removed by melting from below as early as February or March. This is in agreement with existing observations. 相似文献
2.
《Ocean Modelling》2009,28(3-4):114-129
A newly developed global Finite Element Sea Ice–Ocean Model (FESOM) is presented. The ocean component is based on the Finite Element model of the North Atlantic (FENA) but has been substantially updated and extended. In addition to a faster realization of the numerical code, state-of-the-art parameterizations of subgrid-scale processes have been implemented. A Redi/GM scheme is employed to parameterize the effects of mesoscale eddies on lateral tracer distribution. Vertical mixing and convection are parameterized as a function of the Richardson number and the Monin–Obukhov length. A finite element dynamic-thermodynamic sea ice–model has been developed and coupled to the ocean component. Sea ice thermodynamics have been derived from the standard AWI sea ice model featuring a prognostic snow layer but neglecting internal heat storage. The dynamic part offers the viscous-plastic and elastic-viscous-plastic rheologies. All model components are discretized on a triangular/tetrahedral grid with a continuous, conforming representation of model variables. The coupled model is run in a global configuration and forced with NCEP daily atmospheric reanalysis data for 1948–2007. Results are analysed with a slight focus on the Southern Hemisphere. Many aspects of sea ice distribution and hydrography are found to be in good agreement with observations. As in most coarse-scale models, Gulf Stream transport is underestimated, but transports of the Kuroshio and the Antarctic Circumpolar Current appear realistic. The seasonal cycles of Arctic and Antarctic sea ice extents and Antarctic sea ice thickness are well captured; long- and short-term variability of ice coverage is found to be reproduced realistically in both hemispheres. The coupled model is now ready to be used in a wide range of applications. 相似文献
3.
利用Modular Ocean Model version4(MOM4)海洋环流模式设计了2个全球海洋-海冰耦合的数值实验,以分析比较Pacanowski and Philander(PP)和K-Profile Parameterization(KPP)两种不同垂直混合方案在全球海洋上层模拟中的表现。实验结果表明,PP和KPP方案在中纬和高纬海域模拟海温差别较大,后者模拟结果好于前者;在低纬海域差别较小,但赤道断面PP方案模拟结果较好;KPP方案能很好地模拟太平洋赤道潜流,而PP方案模拟的赤道潜流位置偏浅。 相似文献
4.
《Ocean Modelling》2011,36(4):304-313
We implemented an explicit forcing of the complete lunisolar tides into an ocean model which is part of a coupled atmosphere–hydrology–ocean–sea ice model. An ensemble of experiments with this climate model shows that the model is significantly affected by the induced tidal mixing and nonlinear interactions of tides with low frequency motion. The largest changes occur in the North Atlantic where the ocean current system gets changed on large scales. In particular, the pathway of the North Atlantic Current is modified resulting in improved sea surface temperature fields compared to the non-tidal run. These modifications are accompanied by a more realistic simulation of the convection in the Labrador Sea. The modification of sea surface temperature in the North Atlantic region leads to heat flux changes of up to 50 W/m2. The climate simulations indicate that an improvement of the North Atlantic Current has implications for the simulation of the Western European Climate, with amplified temperature trends between 1950 and 2000, which are closer to the observed trends. 相似文献
5.
The sensitivity of the North Atlantic gyre circulation to high latitude buoyancy forcing is explored in a global, non-eddy resolving ocean general circulation model. Increased buoyancy forcing strengthens the deep western boundary current, the northern recirculation gyre, and the North Atlantic Current, which leads to a more realistic Gulf Stream path. High latitude density fluxes and surface water mass transformation are strongly dependent on the choice of sea ice and salinity restoring boundary conditions. Coupling the ocean model to a prognostic sea ice model results in much greater buoyancy loss in the Labrador Sea compared to simulations in which the ocean is forced by prescribed sea ice boundary conditions. A comparison of bulk flux forced hindcast simulations which differ only in their sea ice and salinity restoring forcings reveals the effects of a mixed thermohaline boundary condition transport feedback whereby small, positive temperature and salinity anomalies in subpolar regions are amplified when the gyre spins up as a result of increased buoyancy loss and convection. The primary buoyancy flux effects of the sea ice which cause the simulations to diverge are ice melt, which is less physical in the diagnostic sea ice model, and insulation of the ocean, which is less physical with the prognostic sea ice model. Increased salinity restoring ensures a more realistic net winter buoyancy loss in the Labrador Sea, but it is found that improvements in the Gulf Stream simulation can only be achieved with the excessive buoyancy loss associated with weak salinity restoring. 相似文献
6.
The resolution of the sea-ice component of a coarse-resolution global ocean general circulation model (GCM) has been enhanced to about 22 km in the Southern Ocean. The ocean GCM is designed for long-term integrations suitable for investigations of the deep-ocean equilibrium response to changes in southern hemisphere high-latitude processes. The space and time scales of the high-resolution sea-ice component are commensurate with those of the resolution of satellite passive-microwave sea-ice data. This provides the opportunity for a rigorous evaluation of simulated sea-ice characteristics. It is found that the satellite-derived continuous high ice concentration of the interior winter ice pack can only be captured when vertical oceanic mixing is modified in a way that less local, intermittent convection occurs. Furthermore, the width and the variability of the coastal polynyas around the Antarctic continent and its ice shelves are best captured when some form of ice-shelf melting is accounted for. The width of the wintertime ice edge is reasonably reproduced, while its variability remains underestimated, closely following the coarse-grid pattern of the ocean model due to its high dependence on ocean temperature. Additional variability besides daily winds, e.g. in form of idealized tidal currents, improves the temporal and spatial ice-edge variability, while leads in the interior ice pack become more abundant, more in line with the fine-scale satellite-derived texture. The coast- or ice-shelf line is described on the fine grid based on satellite passive-microwave data. This method requires parts of a coarse coastal ocean grid cell to be covered by an inert layer of “fast ice” or “ice shelf”. Reasonable long-term global deep-ocean properties can only be achieved when these areas are not inert, i.e. are exposed to heat flux and ice growth, or when the vertical mixing parameterization allows for excessive open-ocean convection. The model area exposed to cold high-latitude atmospheric conditions thus being most decisive for a realistic representation of the long-term deep-ocean properties, suggests that high-latitude coastlines are definitely in need of being represented at high resolution, including ice sheets and their effects on the heat and freshwater flux for the ocean. 相似文献
7.
A coding error in the s-Coordinate Primitive Equation Model (SPEM) has led to misleading statements about the behaviour of the Mellor–Yamada level 2 parameterization of vertical mixing. It has been claimed that the scheme removes static instability only very slowly and preserves statically unstable stratifications for an unrealistic long time. This note corrects this statement by demonstrating that the Mellor–Yamada mixing scheme, if implemented correctly, tends to overestimate rather than underestimate vertical mixing in seasonally ice-covered seas. Similar to other mixing schemes with the same behaviour, this leads to spurious open ocean deep convection, an unrealistic homogenization of the water column, and a significant reduction of sea ice volume. 相似文献
8.
《Deep Sea Research Part II: Topical Studies in Oceanography》1999,46(6-7):1335-1355
A three-dimensional hydrodynamic ocean model coupled to a thermohydrodynamic model for young sea ice is applied to study shallow haline convection in the central Greenland Sea, with an emphasis on sub-mesoscale ice–ocean interactions. Two types of young sea ice are distinguished; i.e., frazil and pancake ice, both acting different on surface heat, salt, and momentum fluxes. Two scenarios are considered: (a) continued frazil-ice production during steady winds, and (b) the same scenario but with the intermittent formation of pancake ice during a short intervening period of low winds. Brine release due to new-ice production creates shallow convection in both cases. Under conditions of continued frazil-ice production, ice streaks develop at the sea surface, finally becoming oriented roughly parallel to the wind. These streaks are the result of convective plumes that induce organized patterns of convergent and divergent surface currents. Frazil-ice is herded into convergence zones where it becomes as thick as 6 m within 24 h. The studies suggest a strong relationship between streak spacing and the penetration depth of convection, given by an aspect ratio in the range of 2–3. After pancake ice has been formed, however, the organized ice streaks vanish, developing into complex patterns of pancake ice. This finding is in agreement with recent field observations in the Greenland Sea Is-Odden ice tongue. With the existence of pancake ice, moreover, the surface-averaged buoyancy flux decreases and is determined from the integral of local sub-mesoscale ice–ocean interactions. 相似文献
9.
Mikitoshi Hirabara Hiroyuki Tsujino Hideyuki Nakano Goro Yamanaka 《Journal of Oceanography》2012,68(5):771-796
An experiment using a global ocean–ice model with an interannual forcing data set was conducted to understand the variability in the Southern Ocean. A winter-persisting polynya in the Weddell Sea (the Weddell Polynya, WP) was simulated. The process of WP breaking out after no-WP years was explored using the successive WPs found in the late 1950s. The results suggested that the anomalously warm deep water, saline surface layer, and a cyclonic wind stress over the Maud polynya region in early winter are essential for the surface layer to be dense enough to trigger deep convections which maintain a winter-persisting polynya; also, the reanalyzed surface air temperature (SAT) over the observed polynya region is too high for an ocean–ice model’s bulk formula to yield sufficient upward heat fluxes to induce WP formation. Therefore the Weddell Polynya, a series of WPs observed from satellite in the mid-1970s, is reproduced by replacing the SAT with a climatological one. Subsequent to the successive WP events, density anomalies excited in the Weddell Sea propagate northward in the Atlantic deep basins. The Antarctic Circumpolar Current (ACC) is enhanced through the increased meridional density gradient. The enhanced ACC and its meandering over the abyssal ridges excite buoyancy anomalies near the bottom at the southwestern end of the South Pacific basin. The buoyancy signals propagate northward and eventually arrive in the northern North Pacific. 相似文献
10.
《Ocean Modelling》2002,4(3-4):291-311
Coupled general circulation models (GCMs) have had weak El Niño/Southern Oscillation variability that has been attributed to a diffuse thermocline in the modeled equatorial Pacific Ocean. Consequently, there have been many attempts to improve the thermocline by developing new or improved ocean vertical mixing schemes. This paper investigates the influence of gradient Richardson Number-based vertical mixing scheme profiles in a tropical Pacific Ocean GCM. It has been common for vertical mixing schemes to be assessed in tropical Pacific Ocean models that have a limited latitudinal domain bounded by zonal walls with sponge layers. However, recent work has shown that warm surface water can accumulate in these models and stop them from achieving the observed sharp equatorial thermocline. The present model employs a parameterized wall heat transport scheme that prevents warm surface water from accumulating. Thus we are able assess the influence of vertical mixing profiles in an ocean model that does not allow warm surface water to accumulate and influence the thermocline.In this paper we evaluate the equatorial performance of three different Richardson number (Ri)-based vertical mixing profiles: an integer power (IP) profile based on the observations of Peters, Gregg and Toole; a form of the Pacanowski and Philander profile modified to have low background mixing; and the Max Planck Institute profile. With the accumulation of warm surface water prevented, each of these profiles is able to achieve a sharp thermocline. When compared with observations, the IP profile achieves a better upwelling velocity distribution. We also examine the influence on equatorial performance of very high mixing coefficients at low Richardson number, and of low background mixing coefficients. 相似文献
11.
The near-inertial waves (NIWs) are important for energy cascade in the ocean. They are usually significantly reinforced by strong winds, such as typhoon. Due to relatively coarse resolutions in contemporary climate models, NIWs and associated ocean mixing need to be parameterized. In this study, a parameterization for NIWs proposed by Jochum in 2013 (J13 scheme), which has been widely used, is compared with the observations in the South China Sea, and the observations are treated as model outputs. Under normal conditions, the J13 scheme performs well. However, there are noticeable discrepancies between the J13 scheme and observations during typhoon. During Typhoon Kalmaegi in 2014, the inferred value of the boundary layer is deeper in the J13 scheme due to the weak near-inertial velocity shear in the vertical. After typhoon, the spreading of NIWs beneath the upper boundary layer is much faster than the theoretical prediction of inertial gravity waves, and this fast process is not rendered well by the J13 scheme. In addition, below the boundary layer, NIWs and associated diapycnal mixing last longer than the direct impacts of typhoon on the sea surface. Since the energy dissipation and diapycnal mixing below the boundary layer are bounded to the surface winds in the J13 scheme, the prolonged influences of typhoon via NIWs in the ocean interior are missing in this scheme. Based on current examination, modifications to the J13 scheme are proposed, and the modified version can reduce the discrepancies in the temporal and vertical structures of diapycnal mixing. 相似文献
12.
In the Weddell Sea during the winters of 1974–1976 a significant opening in the sea-ice cover occurred in the vicinity of a large bathymetric feature — the Maud Rise seamount. The event is commonly referred to as the Weddell Polynya. Aside from such a large-scale, relatively persistent polynya in the Weddell Sea, transient, small-scale polynya can also appear in the sea-ice cover at various times throughout the winter and at various locations with respect to the Maud Rise. The underlying causes for the occurrence of such transient polynya have not been unambiguously identified. We hypothesize that variations in the mean ocean currents are one major contributor to such variability in the sea-ice cover. Analysis of the sea-ice equations with certain idealized patterns of ocean currents serving as forcing is shown to lead to Ekman transports of sea ice favorable to the initiation of transient polynya. Aside from the actual spatial pattern of the idealized ocean currents, many other factors need also be taken into account when looking at such transient polynya. Two other such factors discussed are variations in the sea-ice thickness field and the treatment of the sea-ice rheology. Simulations of a sea-ice model coupled to a dynamical ocean model show that the interaction of (dynamical) oceanic currents with large-scale topographic features, such as the Maud Rise, does lead to the formation of transient polynya, again through Ekman transport effects. This occurs because the seamount has a dynamic impact on the three-dimensional oceanic flow field all the way up through the water column, and hence on the near surface ocean currents that are in physical contact with the sea ice. Further simulations of a sea-ice model coupled to a dynamic ocean model and forced with atmospheric buoyancy fluxes show that transient polynya can be enhanced when atmospheric cooling provides a positive feedback mechanism allowing preferential open-ocean convection to occur. The convection, which takes hold at sites where transient polynya have been initiated by sea-ice–ocean stress interaction, has an enhancing effect arising from the convective access to warmer, deeper waters. To investigate all of these effects in a hierarchical manner we use a primitive equation coupled sea-ice–ocean numerical model configured in a periodic channel domain with specified atmospheric conditions. We show that oceanic flow variability can account for temporal variability in small-scale, transient polynya and thus point to a plausible mechanism for the initiation of large-scale, sustained polynya such as the Weddell Polynya event of the mid 1970s. 相似文献
13.
《Ocean Modelling》2002,4(2):137-172
A new sea ice model, GELATO, was developed at Centre National de Recherches Météorologiques (CNRM) and coupled with OPA global ocean model. The sea ice model includes elastic–viscous–plastic rheology, redistribution of ice floes of different thicknesses, and it also takes into account leads, snow cover and snow ice formation. Climatologies of atmospheric surface parameters are used to perform a 20-year global ocean–sea ice simulation, in order to compute surface heat fluxes from diagnosed sea ice or ocean surface temperature. A surface salinity restoring term is applied only to ocean grid cells with no sea ice to avoid significant surface salinity drifts, but no correction of sea surface temperature is introduced. In the Arctic the use of an ocean model substantially improves the representation of sea ice, and particularly of the ice edge in all seasons, as advection of heat and salt can be more accurately accounted for than in the case of, for example, a sea ice–ocean mixed layer model. In contrast, in the Antarctic, a region where ocean convective processes bear a much stronger influence in shaping sea ice characteristics, a better representation of convection and probably of sea ice (for example, of frazil sea ice, brine rejection) would be needed to improve the simulation of the annual cycle of the sea ice cover. The effect of the inclusion of several ice categories in the sea ice model is assessed by running a sensitivity experiment in which only one category of sea ice is considered, along with leads. In the Arctic, such an experiment clearly shows that a multicategory sea ice model better captures the position of the sea ice edge and yields much more realistic sea ice concentrations in most of the region, which is in agreement with results from Bitz et al. [J. Geophys. Res. 106 (C2) (2001) 2441–2463]. 相似文献
14.
Numerical simulations of the Mediterranean sea outflow: impact of the entrainment parameterization in an isopycnic coordinate ocean model 总被引：1，自引：0，他引：1
Gravity current entrainment is essential in determining the properties of the interior ocean water masses that result from marginal sea overflows. Although the individual entraining billows will be unresolvable in large-scale ocean models for the foreseeable future, some large-scale simulations are now being carried out that do resolve the intermediate scale environment which may control the rate of entrainment. Hallberg [Mon. Wea. Rev. 128 (2000) 1402] has recently developed an implicit diapycnal mixing scheme for isopycnic coordinate ocean models that includes the Richardson number dependent entrainment parameterization of Turner [J. Fluid Mech. 173 (1986) 431], and which may be capable of representing the gravity current evolution in large-scale ocean models. The present work uses realistic regional simulations with the Miami Isopycnic Coordinate Ocean Model (MICOM) to evaluate ability of this scheme to simulate the entrainment that is observed to occur in the bottom boundary currents downstream of the Mediterranean outflow. These simulations are strikingly similar to the observations, indicating that this scheme does produce realistic mixing between the Mediterranean outflow and the North Atlantic Central Water. Sensitivity studies identify the critical Richardson number below which vigorous entrainment occurs as a particularly important parameter. Some of these experiments also show meddies detaching from the Mediterranean undercurrent at locations that appear to be highly influenced by topographic features. 相似文献
15.
《Ocean Modelling》2010,35(3-4):137-149
Passive microwave satellite observations of ice extent and concentration form the foundation of sea ice model evaluations, due to their wide spatial coverage and decades-long availability. Observations related to other model quantities are somewhat more limited but increasing as interest in high-latitude processes intensifies. Sea ice thickness, long judged a critical quantity in the physical system, is now being scrutinized more closely in sea ice model simulations as more expansive measurements become available. While albedo is often the first parameter chosen by modelers to adjust simulated ice thickness, this paper explores a set of less prominent parameters to which thickness is also quite sensitive. These include parameters associated with sea ice conductivity, mechanical redistribution, oceanic heat flux, and ice–ocean dynamic stress, in addition to shortwave radiation. Multiple combinations of parameter values can produce the same mean ice thickness using the Los Alamos Sea Ice Model, CICE. One of these “tuned” simulations is compared with a variety of observational data sets in both hemispheres. While deformed ice area compares well with the limited observations available for ridged ice, thickness measurements differ such that the model cannot agree with all of them simultaneously. Albedo and ice–ocean dynamic parameters that affect the turning of the ice relative to the ocean currents have the largest effect on ice thickness, of the parameters tested here. That is, sea ice thickness is highly sensitive to changes in external forcing by the atmosphere or ocean, and therefore serves as a sensitive diagnostic for high-latitude change. 相似文献
16.
《Deep Sea Research Part I: Oceanographic Research Papers》2001,48(5):1255-1277
Isopycnic and cartesian model simulations for the overflow and spreading of dense water are compared with each other and with independent transient-tracer observations. This case study is performed for Adriatic dense water overflowing into the deep eastern Mediterranean with chlorofluoromethane (CFC-12) observations used to test the model simulations. The realism of both types of model simulation depends on the representation of diapycnal mixing. In the cartesian model, convective adjustment and mixing dilute the overflow of Adriatic dense water and lead to unrealistic vertical homogenization. Incorporating a modified convection scheme emphasizing the sinking of dense fluid, rather than its mixing, leads to a more realistic penetration of the dense overflow. In the isopycnic model, there is an improved simulation of the overflow, which leads to the density contrast of the deep Mediterranean waters being maintained. However, there is too low a CFC-12 concentration at mid-depths unless explicit diapycnal mixing is incorporated. In each model, the different spreading of dense water is associated with a different bottom pressure torque and depth-integrated transport, and hence with contrasting tracer distributions throughout the water column. 相似文献
17.
一种海冰热力过程参数化方案 总被引：5，自引：1，他引：4
由于海冰热力过程太复杂,难于精确计算冰面和水面热量收支方程中的每一项,并且净热量收支比每项小得多,因此难以确定海冰热力一动力模式中的净热通量。本文根据渤海的水文气象观测,详细分析了太阳短波辐射、长波辐射、云量、感热和潜热等对海冰热力增长函数的贡献,给出了一种海冰热力过程的参数化方案。并选取２个典型的个例进行了对比研究。模拟结果表明,该参数化方案能较好地模拟渤海海冰的热力过程。 相似文献
18.
夏季黄海表面冷水对大气边界层及海雾的影响 总被引：3，自引：0，他引：3
海表面温度(SST)是海气界面上的1个物理量,受到海洋潮汐、海底地形等因素影响,并对海洋大气边界层有着重要的影响.夏季的黄海,由于黄海冷水团的存在和陆架锋的影响,或是潮汐混合的作用导致海水的垂直混合,使海表面温度的分布产生复杂的结构.通过对卫星观测的海表面温度数据分析,发现在夏季黄海有几个SST冷中心的存在:辽东半岛以及山东半岛的顶端、朝鲜半岛的西侧、山东半岛南侧、江苏外海和黄海南部等.本文利用一系列船舶观测资料、卫星遥感数据、再分析数据分析等,并运用数值模拟研究黄海的冷中心对其上大气的影响.在冷区之上,大气稳定度增加,抑制了近海面大气的垂直混合,使海表面风速减弱.通过对船测数据的分析,在冷区位置有海雾多发区的存在,黄海南部冷区上的海雾发生频率达到15％以上.Weather Research and Forecasting(WRF)模式的数值模拟表明,冷中心降低上空的温度,使海表面风速减弱,形成厚度达500m的逆温层,为海雾的形成创造了有利的条件.与船测数据结果所不同的是黄海南部冷中心之上的海雾发生频率可以达到30％,去掉冷区影响的试验表明冷区较冷的海表面温度最多可以使海雾的发生频率增加15％以上. 相似文献
19.
潮致混合对海洋环流的调整起着重要作用。陆架环流的数值模拟中如果不考虑潮汐作用，往往不能得到与观测相符的垂向温盐结构。本文基于调和分析方法，建立了一套潮致混合参数化方案。该方案通过对垂向混合系数进行调和分析，从而得到随时间变化的潮致混合系数。将该方案用于黄海冷水团数值模拟的结果显示，其能够得到与在数值模式开边界直接加入潮汐强迫相当的冷水团温盐结构。和直接引入潮汐强迫相比，这一潮致混合参数化方案的优势在于，它能够大大节省数值模拟计算机时，因此有望显著提高大规模高分辨率的海洋环流及气候模式的模拟能力和计算效率。 相似文献
20.
《Ocean Modelling》2003,5(2):157-170
Model results from a regional model (BRIOS) of the Southern Ocean that includes ice shelf cavities and the interaction between ocean and ice shelves are used to derive a simple parameterization for ice shelf melting and the corresponding fresh water flux in large-scale ocean climate models. The parameterization assumes that the heat loss and fresh water gain due to the ice shelves are proportional to the difference in freezing temperature at the ice shelf edge base and the oceanic temperature on the shelf/slope area of the adjacent ocean as well as an effective area of interaction. This area is proportional to the along-shelf width of ice shelf and an effective cross-shelf distance, which turns out to be rather uniform (5–15 km) for a variety of different ice shelves. The proposed parameterization is easy to implement and valid for a wide range of circumstances. An application of the proposed scheme in a global ice ocean model (CLIO) supports our hypothesis that it can be used successfully and improves both the ocean and sea ice component of the model. This parameterization should also be used in models of the climate system that include a coupling between an ice sheet and an oceanic component. 相似文献