{"id":62,"date":"2014-06-06T13:16:21","date_gmt":"2014-06-06T17:16:21","guid":{"rendered":"http:\/\/web.whoi.edu\/famos\/?page_id=62"},"modified":"2014-06-06T13:16:21","modified_gmt":"2014-06-06T17:16:21","slug":"freshwater-heat-and-volume-budgets","status":"publish","type":"page","link":"https:\/\/web.whoi.edu\/famos\/freshwater-heat-and-volume-budgets\/","title":{"rendered":"Freshwater, Heat and Volume Budgets"},"content":{"rendered":"

Arctic Ice and Ocean Heat and Freshwater Fluxes: A FAMOS Model Intercomparison Project<\/strong><\/p>\n

\u00a0Sheldon Bacon and Yevgeny Aksenov
\nNational Oceanography Centre, Southampton, UK; November 2013.<\/p>\n

1.\u00a0 Introduction<\/strong><\/p>\n

It is important to understand how models perform in comparison with each other, and in comparison with the real world.\u00a0 One way to do this is to define a metric that can be objectively quantified in all circumstances (in principle).\u00a0 Over the last few years, we have achieved this in the real world by the combination of a set of Arctic Ocean measurements into a near-complete boundary monitoring array.\u00a0 Together with land, this array completely encircles the Arctic Ocean, such that the principle of mass conservation can be applied, and therefore that real fluxes of heat and freshwater can be calculated, independent of any assumptions about reference values.\u00a0 The procedure is described in Tsubouchi et al. (JGR, 2012).<\/p>\n

Reduced to the very simplest representation, warm waters enter the Arctic;\u00a0 these waters leave the Arctic cooler, and the export includes some sea ice.\u00a0 These waters are also fresher when they leave than when they arrive.\u00a0 The cooling and freshening is caused by surface fluxes.\u00a0 Therefore on some timescale, the net surface fluxes should equal the ice and ocean boundary fluxes.\u00a0 The question of timescales, however, is a difficult one due to the great range of residence times within the Arctic for water masses.\u00a0 In the Fram Strait recirculation, some Atlantic water may only cross the boundary for months, or maybe just weeks, before heading back south.\u00a0 At the opposite extreme, the deep waters of the Canada Basin may not be refreshed for centuries.\u00a0 With model runs of typical duration of order decades, it will be interesting to see how the ocean circulation filters the surface fluxes to produce the boundary fluxes, and to compare the amplitudes and phases of the modelled annual cycles of surface and boundary fluxes.<\/p>\n

We note here three important assumptions that enable the analysis to proceed.\u00a0 First, that an appropriate timescale long enough to allow mass conservation but short enough to be useful is one month; so the outputs will be time series of monthly-mean fluxes (Tsubouchi et al., 2012).\u00a0 Second, that Fury and Hecla Strait is closed;\u00a0 if a model does indeed open this strait, then it should be included;\u00a0 but see Tsubouchi et al. (2012) for the rationale behind its exclusion.\u00a0 Third, while a timescale of one month allows mass conservation, it says nothing about heat or freshwater (or salinity) storage.\u00a0 While there have been publications describing storage changes on decadal timescales, we have as yet no information on the nature of storage on monthly timescales.\u00a0 Therefore we can treat storage as an uncertainty (again, as in Tsubouchi et al., 2012).\u00a0 However, if time and resources permit, we ask that participants in this experiment calculate month-to-month changes in mass, heat and salinity storage.<\/p>\n

Finally, before describing the method, we note that with this approach to model analysis, we very much welcome the inclusion of results both from forced and from coupled models, and also from both global and regional models \u2013\u00a0as long as the latter span the defined region.\u00a0 We also welcome results for multiple runs of the same model, where used (eg) to study the impact of parameter changes.<\/p>\n


\n2.\u00a0 Method<\/strong><\/p>\n

2.1\u00a0 Geography<\/em><\/p>\n

The actual measurement locations are shown in Tsubouchi et al. (2012), which is reproduced here as figure 1.\u00a0 We slightly simplify the definition of the boundary, in order that it conforms to meridians of longitude and parallels of latitude.\u00a0 The Arctic Ocean is enclosed by sections across the four main gateways \u2013 Fram, Davis and Bering Straits, and the Barents Sea Opening (BSO).\u00a0 We define the locations of these gateways as follows.<\/p>\n

Fram Strait:\u00a0 coast-to-coast, Greenland to Svalbard, across 79\u02daN<\/p>\n

Davis Strait:\u00a0 coast-to-coast, Baffin Island to Greenland, across 67\u02daN<\/p>\n

Bering Strait:\u00a0 coast-to-coast, Siberia to Alaska, across 65.75\u02daN<\/p>\n

BSO:\u00a0 coast-to-coast, Norway to Svalbard, up 19\u02daE.<\/p>\n

 <\/p>\n

To calculate the balance for the fresh water and heat in the Arctic, the models are required to produce model masks for the surface flux calculations to enclose the vertical model transects.\u00a0 Note that the end points of the sections should be chosen according the individual model coastlines to cover the straits from land to land. Models which need to use \u201cstaircase\u201d transects should check that their positions do not deviate significantly from the observational transects.\u00a0 If tolerance becomes a problem, we would be happy to discuss this further. Included in the transect model cells should be tracer (ie, temperature and salinity) points. \u00a0For staggered grids, the velocity used in the transport calculations should be interpolated in a consistent manner with the calculations used in the model runs under analysis. \u00a0For regional models with open boundaries in Bering Strait, the model transect in the strait should be away from the model boundary and at least one model row inside the domain. \u00a0We cannot accommodate models which have open boundaries north of Davis Strait.<\/p>\n

Specify the surface flux mask for the area enclosed by the transects (as above). \u00a0For \u201cstaircase\u201d transects, the border points of the mask should follow the \u201cstaircases\u201d.\u00a0 Care should be taken to match the transects (as the boundaries) and border points of the mask.<\/p>\n


\n2.2 \u00a0Surface Fluxes<\/em><\/p>\n

The aim is to produce time series of monthly means of surface heat and freshwater fluxes (using the same time base as for boundary fluxes, \u00a72.3 below).\u00a0 First, specify whether forced or coupled model is used.\u00a0 If forced, then specify the forcing set, any modifications to it, and procedures involved in converting input data to surface fluxes (if applicable).\u00a0 Most models will store surface heat and freshwater fluxes in some form, either as integrals or possibly as components (eg latent, sensible, long-wave, short-wave for heat;\u00a0 or precipitation, evaporation, runoff for freshwater).<\/p>\n

If the model uses a surface salinity restoring (relaxation) term, it should be included in the freshwater balances. \u00a0The currently commonly-used approach for the salinity restoring is to correct surface (top model level\/layer, k=1 below) salinity towards observational (climatological) values by adding a virtual freshwater flux \u201cFWcor\u201d to the surface freshwater forcing of the ocean, ie:<\/p>\n

FWcor=gamma*dz*[Sal(k=1) \u2013 Sal_obs]\/Sal(k=1)<\/p>\n

Here, Sal(k=1) is the simulated top level salinity, Sal_obs is the observed top level salinity, dz is the thickness of the top model box and gamma is the restoring timescale.\u00a0 Following this technique, FWcor is added to the freshwater balance as an additional surface freshwater forcing without changing the volume of the ocean. \u00a0For other restoring techniques, a different approach may need to be discussed.\u00a0 If other freshwater balancing components are applied (this for example concerns regional models which may have \u201cdrying effects\u201d of the ocean due to imbalance in atmospheric fluxes), then they need to be included in the total surface fresh water fluxes.<\/p>\n

2.3\u00a0 Ice and Ocean Fluxes<\/em><\/p>\n

The aim is to produce time series of monthly means of boundary heat and freshwater fluxes (using the same time base as for surface fluxes, \u00a72.2 above).\u00a0 While the main outputs must be vertical-area integrals around the whole circuit, integrals for each of the four straits should also be stored, as diagnostics.\u00a0 We next describe the calculations required for heat flux (first) and freshwater (second).\u00a0 Terms are standard:\u00a0 velocity (v, m\/s), density (rho, kg\/m3), specific heat capacity (cp, J\/kg\/\u02daC), latent heat capacity (L, J\/kg), temperature (temp, \u02daC), potential temperature (theta, \u02daC), salinity (S, no units);\u00a0 note that vertical area integrals are dxdz, where z is depth (m) and x the along-section coordinate (m), from the sea bed up to the surface, with due allowance made for sea ice, as described in context below.<\/p>\n

Note that for linear free surface models, sea surface heights are not included in the volume transport calculations, but only in the freshwater and heat transport calculations. For non-linear free surface models, the sea surface heights are included in the volume transport calculations as well.<\/p>\n

2.3.1\u00a0 Heat fluxes<\/p>\n

(i)\u00a0 For the liquid (seawater, sensible heat) component, calculate the area integral of:
\ncp[seawater] x rho x v x theta
\n(ii)\u00a0 For the solid (sea ice, latent heat) component, calculate the area integral of:
\nv[sea-ice] x rho[sea-ice] x L<\/p>\n

 <\/p>\n

(iii)\u00a0 For the solid (sea ice, sensible heat) component, this calculation will depend on sea ice model implementation, ie, whether the sea ice is (a) embedded, or (b) \u201clevitating\u201d.\u00a0 If embedded, then include the ice temperature in the liquid integral, so add the area integral of:<\/p>\n

cp[sea-ice] x rho[sea-ice] v[sea-ice] x temp[sea-ice]<\/p>\n

If the sea ice is levitating, we need the difference between the ocean temperature where the ice \u201cshould be\u201d and modelled ice temperature, so, for the cross-sectional area of the sea ice, add the area integral of:<\/p>\n

cp[sea-ice] x rho[sea-ice] x v[sea-ice] x temp[sea-ice] minus cp[seawater] x rho x v x theta<\/p>\n

(iv)\u00a0 Add all three components (one liquid plus two solid) together.\u00a0 Note the possibility of variable cp both for seawater and sea ice.<\/p>\n

2.3.2\u00a0 Freshwater fluxes<\/p>\n

(i)\u00a0 For each month, calculate total vertical boundary area mean salinity Sbar, including both ice and seawater, and also including snow.<\/p>\n

(ii)\u00a0 The subsequent calculation again differs depending on whether sea ice is embedded or levitating.\u00a0 The embedded case is straightforward.\u00a0 For the liquid component, compute:<\/p>\n

the area integral of (v x S) plus area integral of (v[sea-ice] x S[sea-ice]), all divided by Sbar<\/p>\n

For the levitating case, for the solid component, for the cross-sectional area of sea ice, calculate the area integral of (v[sea-ice] x S[sea-ice] minus [v x S]).\u00a0 Add the result to the area integral of (v x S) for the liquid component.<\/p>\n

2.4 \u00a0Time Period<\/em><\/p>\n

We hope participants will analyse the whole of their model runs, assuming these to be generally multi-decadal.\u00a0 If this is not possible, we wish to focus first on the modern period (after 1950), and second on 2000-2010, as this period is best for observations.\u00a0 We are happy to discuss this.<\/p>\n

2.5\u00a0 Storage<\/em><\/p>\n

Seawater mass storage, heat storage and salinity storage should be calculated as the net changes of ocean mass, temperature and salinity for both liquid and solid components (therefore including sea ice), from month to month.\u00a0 It may be that we need to correspond about these issues.<\/p>\n

3. \u00a0Outcomes<\/strong><\/p>\n

When ready, please send your time series products to Bacon & Aksenov.\u00a0 We aim for collation of results in the hope that initial results, at least, might be ready for presentation to the Royal Society Sea Ice meeting being organised by Daniel Feltham and to be held in London in September 2014.\u00a0 We expect a subsequent multi-author publication to appear in Proc. Roy. Soc. A.<\/p>\n

If all goes well, we would like to discuss next steps at the FAMOS 2014 meeting.<\/p>\n

 <\/p>\n

Figure 1:\u00a0 Bathymetric configuration in Davis, Fram, and Bering Straits, and the Barents Sea Opening (BSO), showing CTD stations (red cross), model grid points (green cross), mooring locations (blue diamond), and station numbers (including model grid points). Bathymetric contour intervals (CI) are shown for each strait; the CI for the Arctic panel is 1000 m.\u00a0 From Tsubouchi et al. (2012).<\/p>\n

Reference<\/strong><\/p>\n

Tsubouchi, T., S. Bacon, A. C. Naveira Garabato, Y. Aksenov, S. W. Laxon, E. Fahrbach, A. Beszczynska-M\u00f6ller, E. Hansen, C. M. Lee and R. B. Ingvaldsen, 2012: The Arctic Ocean in summer:\u00a0 a quasi-synoptic inverse estimate of boundary fluxes and water mass transformation.\u00a0 J. Geophys. Res.<\/em>, 117<\/strong>, C01024, <\/em>doi:10.1029\/2011JC007174.<\/p>\n","protected":false},"excerpt":{"rendered":"

Arctic Ice and Ocean Heat and Freshwater Fluxes: A FAMOS Model Intercomparison Project \u00a0Sheldon Bacon and Yevgeny Aksenov National Oceanography Centre, Southampton, UK; November 2013. 1.\u00a0 Introduction It is important to understand how models perform in comparison with each other, and in comparison with the real world.\u00a0 One way to do this is to define.<\/p>\n","protected":false},"author":42,"featured_media":0,"parent":0,"menu_order":0,"comment_status":"closed","ping_status":"closed","template":"","meta":[],"_links":{"self":[{"href":"https:\/\/web.whoi.edu\/famos\/wp-json\/wp\/v2\/pages\/62"}],"collection":[{"href":"https:\/\/web.whoi.edu\/famos\/wp-json\/wp\/v2\/pages"}],"about":[{"href":"https:\/\/web.whoi.edu\/famos\/wp-json\/wp\/v2\/types\/page"}],"author":[{"embeddable":true,"href":"https:\/\/web.whoi.edu\/famos\/wp-json\/wp\/v2\/users\/42"}],"replies":[{"embeddable":true,"href":"https:\/\/web.whoi.edu\/famos\/wp-json\/wp\/v2\/comments?post=62"}],"version-history":[{"count":1,"href":"https:\/\/web.whoi.edu\/famos\/wp-json\/wp\/v2\/pages\/62\/revisions"}],"predecessor-version":[{"id":63,"href":"https:\/\/web.whoi.edu\/famos\/wp-json\/wp\/v2\/pages\/62\/revisions\/63"}],"wp:attachment":[{"href":"https:\/\/web.whoi.edu\/famos\/wp-json\/wp\/v2\/media?parent=62"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}