| 1 |
|
% $Header$ |
| 2 |
|
% $Name$ |
| 3 |
\documentclass[12pt]{article} |
\documentclass[12pt]{article} |
| 4 |
\usepackage{epsfig} |
|
| 5 |
\usepackage{graphics} |
\usepackage[]{graphicx} |
| 6 |
\usepackage{subfigure} |
\usepackage{subfigure} |
| 7 |
|
|
| 8 |
\usepackage[round,comma]{natbib} |
\usepackage[round,comma]{natbib} |
| 37 |
\newlength{\mediumfigwidth}\setlength{\mediumfigwidth}{39pc} |
\newlength{\mediumfigwidth}\setlength{\mediumfigwidth}{39pc} |
| 38 |
%\newlength{\widefigwidth}\setlength{\widefigwidth}{39pc} |
%\newlength{\widefigwidth}\setlength{\widefigwidth}{39pc} |
| 39 |
\newlength{\widefigwidth}\setlength{\widefigwidth}{\textwidth} |
\newlength{\widefigwidth}\setlength{\widefigwidth}{\textwidth} |
| 40 |
\newcommand{\fpath}{.} |
\newcommand{\fpath}{figs} |
| 41 |
|
|
| 42 |
|
% commenting scheme |
| 43 |
|
\newcommand{\ml}[1]{\textsf{\slshape #1}} |
| 44 |
|
|
| 45 |
\title{A Dynamic-Thermodynamic Sea ice Model for Ocean Climate |
\title{A Dynamic-Thermodynamic Sea ice Model for Ocean Climate |
| 46 |
Estimation on an Arakawa C-Grid} |
Estimation on an Arakawa C-Grid} |
| 52 |
\maketitle |
\maketitle |
| 53 |
|
|
| 54 |
\begin{abstract} |
\begin{abstract} |
| 55 |
Some blabla |
|
| 56 |
|
As part of ongoing efforts to obtain a best possible synthesis of most |
| 57 |
|
available, global-scale, ocean and sea ice data, a dynamic and thermodynamic |
| 58 |
|
sea-ice model has been coupled to the Massachusetts Institute of Technology |
| 59 |
|
general circulation model (MITgcm). Ice mechanics follow a viscous plastic |
| 60 |
|
rheology and the ice momentum equations are solved numerically using either |
| 61 |
|
line successive relaxation (LSR) or elastic-viscous-plastic (EVP) dynamic |
| 62 |
|
models. Ice thermodynamics are represented using either a zero-heat-capacity |
| 63 |
|
formulation or a two-layer formulation that conserves enthalpy. The model |
| 64 |
|
includes prognostic variables for snow and for sea-ice salinity. The above |
| 65 |
|
sea ice model components were borrowed from current-generation climate models |
| 66 |
|
but they were reformulated on an Arakawa C-grid in order to match the MITgcm |
| 67 |
|
oceanic grid and they were modified in many ways to permit efficient and |
| 68 |
|
accurate automatic differentiation. This paper describes the MITgcm sea ice |
| 69 |
|
model; it presents example Arctic and Antarctic results from a realistic, |
| 70 |
|
eddy-permitting, global ocean and sea-ice configuration; it compares B-grid |
| 71 |
|
and C-grid dynamic solvers in a regional Arctic configuration; and it presents |
| 72 |
|
example results from coupled ocean and sea-ice adjoint-model integrations. |
| 73 |
|
|
| 74 |
\end{abstract} |
\end{abstract} |
| 75 |
|
|
| 76 |
\section{Introduction} |
\section{Introduction} |
| 77 |
\label{sec:intro} |
\label{sec:intro} |
| 78 |
|
|
|
more blabla |
|
|
|
|
| 79 |
\section{Model} |
\section{Model} |
| 80 |
\label{sec:model} |
\label{sec:model} |
| 81 |
|
|
| 148 |
\frac{\partial{u_{i}}}{\partial{x_{j}}} + |
\frac{\partial{u_{i}}}{\partial{x_{j}}} + |
| 149 |
\frac{\partial{u_{j}}}{\partial{x_{i}}}\right). |
\frac{\partial{u_{j}}}{\partial{x_{i}}}\right). |
| 150 |
\end{equation*} |
\end{equation*} |
| 151 |
The pressure $P$, a measure of ice strength, depends on both thickness |
The maximum ice pressure $P_{\max}$, a measure of ice strength, depends on |
| 152 |
$h$ and compactness (concentration) $c$: \[P = |
both thickness $h$ and compactness (concentration) $c$: |
| 153 |
P^{*}c\,h\,e^{[C^{*}\cdot(1-c)]},\] with the constants $P^{*}$ and |
\begin{equation} |
| 154 |
$C^{*}$. The nonlinear bulk and shear viscosities $\eta$ and $\zeta$ |
P_{\max} = P^{*}c\,h\,e^{[C^{*}\cdot(1-c)]}, |
| 155 |
are functions of ice strain rate invariants and ice strength such that |
\label{eq:icestrength} |
| 156 |
the principal components of the stress lie on an elliptical yield |
\end{equation} |
| 157 |
curve with the ratio of major to minor axis $e$ equal to $2$; they are |
with the constants $P^{*}$ and $C^{*}$. The nonlinear bulk and shear |
| 158 |
given by: |
viscosities $\eta$ and $\zeta$ are functions of ice strain rate |
| 159 |
|
invariants and ice strength such that the principal components of the |
| 160 |
|
stress lie on an elliptical yield curve with the ratio of major to |
| 161 |
|
minor axis $e$ equal to $2$; they are given by: |
| 162 |
\begin{align*} |
\begin{align*} |
| 163 |
\zeta =& \frac{P}{2\Delta} \\ |
\zeta =& \min\left(\frac{P_{\max}}{2\max(\Delta,\Delta_{\min})}, |
| 164 |
\eta =& \frac{P}{2\Delta{e}^2} \\ |
\zeta_{\max}\right) \\ |
| 165 |
|
\eta =& \frac{\zeta}{e^2} \\ |
| 166 |
\intertext{with the abbreviation} |
\intertext{with the abbreviation} |
| 167 |
\Delta = & \left[ |
\Delta = & \left[ |
| 168 |
\left(\dot{\epsilon}_{11}^2+\dot{\epsilon}_{22}^2\right) |
\left(\dot{\epsilon}_{11}^2+\dot{\epsilon}_{22}^2\right) |
| 170 |
2\dot{\epsilon}_{11}\dot{\epsilon}_{22} (1-e^{-2}) |
2\dot{\epsilon}_{11}\dot{\epsilon}_{22} (1-e^{-2}) |
| 171 |
\right]^{-\frac{1}{2}} |
\right]^{-\frac{1}{2}} |
| 172 |
\end{align*} |
\end{align*} |
| 173 |
|
The bulk viscosities are bounded above by imposing both a minimum |
| 174 |
|
$\Delta_{\min}=10^{-11}\text{\,s}^{-1}$ (for numerical reasons) and a |
| 175 |
|
maximum $\zeta_{\max} = P_{\max}/\Delta^*$, where |
| 176 |
|
$\Delta^*=(5\times10^{12}/2\times10^4)\text{\,s}^{-1}$. For stress |
| 177 |
|
tensor compuation the replacement pressure $P = 2\,\Delta\zeta$ |
| 178 |
|
\citep{hibler95} is used so that the stress state always lies on the |
| 179 |
|
elliptic yield curve by definition. |
| 180 |
|
|
| 181 |
|
In the so-called truncated ellipse method the shear viscosity $\eta$ |
| 182 |
|
is capped to suppress any tensile stress \citep{hibler97, geiger98}: |
| 183 |
|
\begin{equation} |
| 184 |
|
\label{eq:etatem} |
| 185 |
|
\eta = \min(\frac{\zeta}{e^2} |
| 186 |
|
\frac{\frac{P}{2}-\zeta(\dot{\epsilon}_{11}+\dot{\epsilon}_{22})} |
| 187 |
|
{\sqrt{(\dot{\epsilon}_{11}+\dot{\epsilon}_{22})^2 |
| 188 |
|
+4\dot{\epsilon}_{12}^2}} |
| 189 |
|
\end{equation} |
| 190 |
|
|
| 191 |
In the current implementation, the VP-model is integrated with the |
In the current implementation, the VP-model is integrated with the |
| 192 |
semi-implicit line successive over relaxation (LSOR)-solver of |
semi-implicit line successive over relaxation (LSOR)-solver of |
| 193 |
\citet{zhang98}, which allows for long time steps that, in our case, |
\citet{zhang98}, which allows for long time steps that, in our case, |
| 327 |
state variables to be advected by ice velocities, namely enthalphy of |
state variables to be advected by ice velocities, namely enthalphy of |
| 328 |
the two ice layers and the thickness of the overlying snow layer. |
the two ice layers and the thickness of the overlying snow layer. |
| 329 |
|
|
|
\section{Funnel Experiments} |
|
|
\label{sec:funnel} |
|
|
|
|
|
\begin{itemize} |
|
|
\item B-grid LSR no-slip |
|
|
\item C-grid LSR no-slip |
|
|
\item C-grid LSR slip |
|
|
\item C-grid EVP no-slip |
|
|
\item C-grid EVP slip |
|
|
\end{itemize} |
|
|
|
|
|
\subsection{B-grid vs.\ C-grid} |
|
|
Comparison between: |
|
|
\begin{itemize} |
|
|
\item B-grid, lsr, no-slip |
|
|
\item C-grid, lsr, no-slip |
|
|
\item C-grid, evp, no-slip |
|
|
\end{itemize} |
|
|
all without ice-ocean stress, because ice-ocean stress does not work |
|
|
for B-grid. |
|
| 330 |
|
|
| 331 |
\subsection{C-grid} |
\subsection{C-grid} |
| 332 |
\begin{itemize} |
\begin{itemize} |
| 373 |
\subsection{Arctic Domain with Open Boundaries} |
\subsection{Arctic Domain with Open Boundaries} |
| 374 |
\label{sec:arctic} |
\label{sec:arctic} |
| 375 |
|
|
| 376 |
The Arctic domain of integration is illustrated in Fig.~\ref{???}. It is |
The Arctic domain of integration is illustrated in Fig.~\ref{fig:arctic1}. It |
| 377 |
carved out from, and obtains open boundary conditions from, the global |
is carved out from, and obtains open boundary conditions from, the |
| 378 |
cubed-sphere configuration of the Estimating the Circulation and Climate of |
global cubed-sphere configuration of the Estimating the Circulation |
| 379 |
the Ocean, Phase II (ECCO2) project \cite{men05a}. The domain size is 420 by |
and Climate of the Ocean, Phase II (ECCO2) project |
| 380 |
384 grid boxes horizontally with mean horizontal grid spacing of 18 km. |
\citet{menemenlis05}. The domain size is 420 by 384 grid boxes |
| 381 |
|
horizontally with mean horizontal grid spacing of 18 km. |
| 382 |
|
|
| 383 |
|
\begin{figure} |
| 384 |
|
%\centerline{{\includegraphics*[width=0.44\linewidth]{\fpath/arctic1.eps}}} |
| 385 |
|
\caption{Bathymetry of Arctic Domain.\label{fig:arctic1}} |
| 386 |
|
\end{figure} |
| 387 |
|
|
| 388 |
There are 50 vertical levels ranging in thickness from 10 m near the surface |
There are 50 vertical levels ranging in thickness from 10 m near the surface |
| 389 |
to approximately 450 m at a maximum model depth of 6150 m. Bathymetry is from |
to approximately 450 m at a maximum model depth of 6150 m. Bathymetry is from |
| 390 |
the National Geophysical Data Center (NGDC) 2-minute gridded global relief |
the National Geophysical Data Center (NGDC) 2-minute gridded global relief |
| 391 |
data (ETOPO2) and the model employs the partial-cell formulation of |
data (ETOPO2) and the model employs the partial-cell formulation of |
| 392 |
\cite{adc97}, which permits accurate representation of the bathymetry. The |
\citet{adcroft97:_shaved_cells}, which permits accurate representation of the bathymetry. The |
| 393 |
model is integrated in a volume-conserving configuration using a finite volume |
model is integrated in a volume-conserving configuration using a finite volume |
| 394 |
discretization with C-grid staggering of the prognostic variables. In the |
discretization with C-grid staggering of the prognostic variables. In the |
| 395 |
ocean, the non-linear equation of state of \cite{jac95}. The ocean model is |
ocean, the non-linear equation of state of \citet{jackett95}. The ocean model is |
| 396 |
coupled to a sea-ice model described hereinabove. |
coupled to a sea-ice model described hereinabove. |
| 397 |
|
|
| 398 |
This particular ECCO2 simulation is initialized from rest using the January |
This particular ECCO2 simulation is initialized from rest using the |
| 399 |
temperature and salinity distribution from the World Ocean Atlas 2001 (WOA01) |
January temperature and salinity distribution from the World Ocean |
| 400 |
[Conkright et al., 2002] and it is integrated for 32 years prior to the |
Atlas 2001 (WOA01) [Conkright et al., 2002] and it is integrated for |
| 401 |
1996-2001 period discussed in the study. Surface boundary conditions are from |
32 years prior to the 1996--2001 period discussed in the study. Surface |
| 402 |
the National Centers for Environmental Prediction and the National Center for |
boundary conditions are from the National Centers for Environmental |
| 403 |
Atmospheric Research (NCEP/NCAR) atmospheric reanalysis [Kistler et al., |
Prediction and the National Center for Atmospheric Research |
| 404 |
2001]. Six-hourly surface winds, temperature, humidity, downward short- and |
(NCEP/NCAR) atmospheric reanalysis [Kistler et al., 2001]. Six-hourly |
| 405 |
long-wave radiations, and precipitation are converted to heat, freshwater, and |
surface winds, temperature, humidity, downward short- and long-wave |
| 406 |
wind stress fluxes using the Large and Pond [1981, 1982] bulk |
radiations, and precipitation are converted to heat, freshwater, and |
| 407 |
formulae. Shortwave radiation decays exponentially as per Paulson and Simpson |
wind stress fluxes using the \citet{large81, large82} bulk formulae. |
| 408 |
[1977]. Additionally the time-mean river run-off from Large and Nurser [2001] |
Shortwave radiation decays exponentially as per Paulson and Simpson |
| 409 |
is applied and there is a relaxation to the monthly-mean climatological sea |
[1977]. Additionally the time-mean river run-off from Large and Nurser |
| 410 |
surface salinity values from WOA01 with a relaxation time scale of 3 |
[2001] is applied and there is a relaxation to the monthly-mean |
| 411 |
months. Vertical mixing follows Large et al. [1994] with background vertical |
climatological sea surface salinity values from WOA01 with a |
| 412 |
diffusivity of 1.5 × 10-5 m2 s-1 and viscosity of 10-3 m2 s-1. A third order, |
relaxation time scale of 3 months. Vertical mixing follows |
| 413 |
direct-space-time advection scheme with flux limiter is employed and there is |
\citet{large94} with background vertical diffusivity of |
| 414 |
no explicit horizontal diffusivity. Horizontal viscosity follows Leith [1996] |
$1.5\times10^{-5}\text{\,m$^{2}$\,s$^{-1}$}$ and viscosity of |
| 415 |
but modified to sense the divergent flow as per Fox-Kemper and Menemenlis [in |
$10^{-3}\text{\,m$^{2}$\,s$^{-1}$}$. A third order, direct-space-time |
| 416 |
press]. Shortwave radiation decays exponentially as per Paulson and Simpson |
advection scheme with flux limiter is employed \citep{hundsdorfer94} |
| 417 |
[1977]. Additionally, the time-mean runoff of Large and Nurser [2001] is |
and there is no explicit horizontal diffusivity. Horizontal viscosity |
| 418 |
applied near the coastline and, where there is open water, there is a |
follows \citet{lei96} but |
| 419 |
relaxation to monthly-mean WOA01 sea surface salinity with a time constant of |
modified to sense the divergent flow as per Fox-Kemper and Menemenlis |
| 420 |
45 days. |
[in press]. Shortwave radiation decays exponentially as per Paulson |
| 421 |
|
and Simpson [1977]. Additionally, the time-mean runoff of Large and |
| 422 |
|
Nurser [2001] is applied near the coastline and, where there is open |
| 423 |
|
water, there is a relaxation to monthly-mean WOA01 sea surface |
| 424 |
|
salinity with a time constant of 45 days. |
| 425 |
|
|
| 426 |
Open water, dry |
Open water, dry |
| 427 |
ice, wet ice, dry snow, and wet snow albedo are, respectively, 0.15, 0.85, |
ice, wet ice, dry snow, and wet snow albedo are, respectively, 0.15, 0.85, |
| 450 |
\item C-grid LSR slip |
\item C-grid LSR slip |
| 451 |
\item C-grid EVP no-slip |
\item C-grid EVP no-slip |
| 452 |
\item C-grid EVP slip |
\item C-grid EVP slip |
| 453 |
|
\item C-grid LSR + TEM (truncated ellipse method, no tensile stress, new flag) |
| 454 |
\item C-grid LSR no-slip + Winton |
\item C-grid LSR no-slip + Winton |
| 455 |
\item speed-performance-accuracy (small) |
\item speed-performance-accuracy (small) |
| 456 |
ice transport through Canadian Archipelago differences |
ice transport through Canadian Archipelago differences |
| 462 |
\begin{itemize} |
\begin{itemize} |
| 463 |
\item advection schemes: along the ice-edge and regions with large |
\item advection schemes: along the ice-edge and regions with large |
| 464 |
gradients |
gradients |
| 465 |
\item C-grid: more transport through narrow straits for no slip |
\item C-grid: less transport through narrow straits for no slip |
| 466 |
conditons, less for free slip |
conditons, more for free slip |
| 467 |
\item VP vs.\ EVP: speed performance, accuracy? |
\item VP vs.\ EVP: speed performance, accuracy? |
| 468 |
\item ocean stress: different water mass properties beneath the ice |
\item ocean stress: different water mass properties beneath the ice |
| 469 |
\end{itemize} |
\end{itemize} |
| 544 |
checkpointing loop. |
checkpointing loop. |
| 545 |
Again, an initial code adjustment is required to support TAFs |
Again, an initial code adjustment is required to support TAFs |
| 546 |
checkpointing capability. |
checkpointing capability. |
| 547 |
The code adjustments are sufficiently simply so as not to cause |
The code adjustments are sufficiently simple so as not to cause |
| 548 |
major limitations to the full nonlinear parent model. |
major limitations to the full nonlinear parent model. |
| 549 |
Once in place, an adjoint model of a new model configuration |
Once in place, an adjoint model of a new model configuration |
| 550 |
may be derived in about 10 minutes. |
may be derived in about 10 minutes. |
| 567 |
We demonstrate the power of the adjoint method |
We demonstrate the power of the adjoint method |
| 568 |
in the context of investigating sea-ice export sensitivities through Fram Strait |
in the context of investigating sea-ice export sensitivities through Fram Strait |
| 569 |
(for details of this study see Heimbach et al., 2007). |
(for details of this study see Heimbach et al., 2007). |
| 570 |
|
%\citep[for details of this study see][]{heimbach07}. %Heimbach et al., 2007). |
| 571 |
The domain chosen is a coarsened version of the Arctic face of the |
The domain chosen is a coarsened version of the Arctic face of the |
| 572 |
high-resolution cubed-sphere configuration of the ECCO2 project |
high-resolution cubed-sphere configuration of the ECCO2 project |
| 573 |
(see Menemenlis et al. 2005). It covers the entire Arctic, |
\citep[see][]{menemenlis05}. It covers the entire Arctic, |
| 574 |
extends into the North Pacific such as to cover the entire |
extends into the North Pacific such as to cover the entire |
| 575 |
ice-covered regions, and comprises parts of the North Atlantic |
ice-covered regions, and comprises parts of the North Atlantic |
| 576 |
down to XXN to enable analysis of remote influences of the |
down to XXN to enable analysis of remote influences of the |
| 581 |
(benchmarks have been performed both on an SGI Altix as well as an |
(benchmarks have been performed both on an SGI Altix as well as an |
| 582 |
IBM SP5 at NASA/ARC). |
IBM SP5 at NASA/ARC). |
| 583 |
|
|
| 584 |
Following a 1-year spinup, the model has been integrated for four years |
Following a 1-year spinup, the model has been integrated for four |
| 585 |
between 1992 and 1995. |
years between 1992 and 1995. It is forced using realistic 6-hourly |
| 586 |
It is forced using realistic 6-hourly NCEP/NCAR atmospheric state variables. |
NCEP/NCAR atmospheric state variables. Over the open ocean these are |
| 587 |
Over the open ocean these are converted into |
converted into air-sea fluxes via the bulk formulae of |
| 588 |
air-sea fluxes via the bulk formulae of Large and Yeager (2004). |
\citet{large04}. Derivation of air-sea fluxes in the presence of |
| 589 |
Derivation of air-sea fluxes in the presence of sea-ice is handled |
sea-ice is handled by the ice model as described in \refsec{model}. |
|
by the ice model as described in Section XXX. |
|
| 590 |
The objective function chosen is sea-ice export through Fram Strait |
The objective function chosen is sea-ice export through Fram Strait |
| 591 |
computed for December 1995 |
computed for December 1995. The adjoint model computes sensitivities |
| 592 |
The adjoint model computes sensitivities to sea-ice export back in time |
to sea-ice export back in time from 1995 to 1992 along this |
| 593 |
from 1995 to 1992 along this trajectory. |
trajectory. In principle all adjoint model variable (i.e., Lagrange |
| 594 |
In principle all adjoint model variable (i.e. Lagrange multipliers) |
multipliers) of the coupled ocean/sea-ice model are available to |
| 595 |
of the coupled ocean/sea-ice model |
analyze the transient sensitivity behaviour of the ocean and sea-ice |
| 596 |
are available to analyze the transient sensitivity behaviour |
state. Over the open ocean, the adjoint of the bulk formula scheme |
| 597 |
of the ocean and sea-ice state. |
computes sensitivities to the time-varying atmospheric state. Over |
| 598 |
Over the open ocean, the adjoint of the bulk formula scheme |
ice-covered parts, the sea-ice adjoint converts surface ocean |
| 599 |
computes sensitivities to the time-varying atmospheric state. |
sensitivities to atmospheric sensitivities. |
| 600 |
Over ice-covered parts, the sea-ice adjoint converts |
|
| 601 |
surface ocean sensitivities to atmospheric sensitivities. |
\reffig{4yradjheff}(a--d) depict sensitivities of sea-ice export |
| 602 |
|
through Fram Strait in December 1995 to changes in sea-ice thickness |
| 603 |
Fig. XXX(a--d) depict sensitivities of sea-ice export through Fram Strait |
12, 24, 36, 48 months back in time. Corresponding sensitivities to |
| 604 |
in December 1995 to changes in sea-ice thickness |
ocean surface temperature are depicted in |
| 605 |
12, 24, 36, 48 months back in time. |
\reffig{4yradjthetalev1}(a--d). The main characteristics is |
| 606 |
Corresponding sensitivities to ocean surface temperature are |
consistency with expected advection of sea-ice over the relevant time |
| 607 |
depicted in Fig. XXX(a--d). |
scales considered. The general positive pattern means that an |
| 608 |
The main characteristics is consistency with expected advection |
increase in sea-ice thickness at location $(x,y)$ and time $t$ will |
| 609 |
of sea-ice over the relevant time scales considered. |
increase sea-ice export through Fram Strait at time $T_e$. Largest |
| 610 |
The general positive pattern means that an increase in |
distances from Fram Strait indicate fastest sea-ice advection over the |
| 611 |
sea-ice thickness at location $(x,y)$ and time $t$ will increase |
time span considered. The ice thickness sensitivities are in close |
| 612 |
sea-ice export through Fram Strait at time $T_e$. |
correspondence to ocean surface sentivitites, but of opposite sign. |
| 613 |
Largest distances from Fram Strait indicate fastest sea-ice advection |
An increase in temperature will incur ice melting, decrease in ice |
| 614 |
over the time span considered. |
thickness, and therefore decrease in sea-ice export at time $T_e$. |
|
The ice thickness sensitivities are in close correspondence to |
|
|
ocean surface sentivitites, but of opposite sign. |
|
|
An increase in temperature will incur ice melting, decrease in ice thickness, |
|
|
and therefore decrease in sea-ice export at time $T_e$. |
|
| 615 |
|
|
| 616 |
The picture is fundamentally different and much more complex |
The picture is fundamentally different and much more complex |
| 617 |
for sensitivities to ocean temperatures away from the surface. |
for sensitivities to ocean temperatures away from the surface. |
| 618 |
Fig. XXX (a--d) depicts ice export sensitivities to |
\reffig{4yradjthetalev10??}(a--d) depicts ice export sensitivities to |
| 619 |
temperatures at roughly 400 m depth. |
temperatures at roughly 400 m depth. |
| 620 |
Primary features are the effect of the heat transport of the North |
Primary features are the effect of the heat transport of the North |
| 621 |
Atlantic current which feeds into the West Spitsbergen current, |
Atlantic current which feeds into the West Spitsbergen current, |
| 624 |
\begin{figure}[t!] |
\begin{figure}[t!] |
| 625 |
\centerline{ |
\centerline{ |
| 626 |
\subfigure[{\footnotesize -12 months}] |
\subfigure[{\footnotesize -12 months}] |
| 627 |
{\includegraphics*[width=0.44\linewidth]{figs/run_4yr_ADJheff_arc_lev1_tim072_cmax2.0E+02.eps}} |
{\includegraphics*[width=0.44\linewidth]{\fpath/run_4yr_ADJheff_arc_lev1_tim072_cmax2.0E+02.eps}} |
| 628 |
%\includegraphics*[width=.3\textwidth]{H_c.bin_res_100_lev1.pdf} |
%\includegraphics*[width=.3\textwidth]{H_c.bin_res_100_lev1.pdf} |
| 629 |
% |
% |
| 630 |
\subfigure[{\footnotesize -24 months}] |
\subfigure[{\footnotesize -24 months}] |
| 631 |
{\includegraphics*[width=0.44\linewidth]{figs/run_4yr_ADJheff_arc_lev1_tim145_cmax2.0E+02.eps}} |
{\includegraphics*[width=0.44\linewidth]{\fpath/run_4yr_ADJheff_arc_lev1_tim145_cmax2.0E+02.eps}} |
| 632 |
} |
} |
| 633 |
|
|
| 634 |
\centerline{ |
\centerline{ |
| 635 |
\subfigure[{\footnotesize |
\subfigure[{\footnotesize |
| 636 |
-36 months}] |
-36 months}] |
| 637 |
{\includegraphics*[width=0.44\linewidth]{figs/run_4yr_ADJheff_arc_lev1_tim218_cmax2.0E+02.eps}} |
{\includegraphics*[width=0.44\linewidth]{\fpath/run_4yr_ADJheff_arc_lev1_tim218_cmax2.0E+02.eps}} |
| 638 |
% |
% |
| 639 |
\subfigure[{\footnotesize |
\subfigure[{\footnotesize |
| 640 |
-48 months}] |
-48 months}] |
| 641 |
{\includegraphics*[width=0.44\linewidth]{figs/run_4yr_ADJheff_arc_lev1_tim292_cmax2.0E+02.eps}} |
{\includegraphics*[width=0.44\linewidth]{\fpath/run_4yr_ADJheff_arc_lev1_tim292_cmax2.0E+02.eps}} |
| 642 |
} |
} |
| 643 |
\caption{Sensitivity of sea-ice export through Fram Strait in December 2005 to |
\caption{Sensitivity of sea-ice export through Fram Strait in December 2005 to |
| 644 |
sea-ice thickness at various prior times. |
sea-ice thickness at various prior times. |
| 649 |
\begin{figure}[t!] |
\begin{figure}[t!] |
| 650 |
\centerline{ |
\centerline{ |
| 651 |
\subfigure[{\footnotesize -12 months}] |
\subfigure[{\footnotesize -12 months}] |
| 652 |
{\includegraphics*[width=0.44\linewidth]{figs/run_4yr_ADJtheta_arc_lev1_tim072_cmax5.0E+01.eps}} |
{\includegraphics*[width=0.44\linewidth]{\fpath/run_4yr_ADJtheta_arc_lev1_tim072_cmax5.0E+01.eps}} |
| 653 |
%\includegraphics*[width=.3\textwidth]{H_c.bin_res_100_lev1.pdf} |
%\includegraphics*[width=.3\textwidth]{H_c.bin_res_100_lev1.pdf} |
| 654 |
% |
% |
| 655 |
\subfigure[{\footnotesize -24 months}] |
\subfigure[{\footnotesize -24 months}] |
| 656 |
{\includegraphics*[width=0.44\linewidth]{figs/run_4yr_ADJtheta_arc_lev1_tim145_cmax5.0E+01.eps}} |
{\includegraphics*[width=0.44\linewidth]{\fpath/run_4yr_ADJtheta_arc_lev1_tim145_cmax5.0E+01.eps}} |
| 657 |
} |
} |
| 658 |
|
|
| 659 |
\centerline{ |
\centerline{ |
| 660 |
\subfigure[{\footnotesize |
\subfigure[{\footnotesize |
| 661 |
-36 months}] |
-36 months}] |
| 662 |
{\includegraphics*[width=0.44\linewidth]{figs/run_4yr_ADJtheta_arc_lev1_tim218_cmax5.0E+01.eps}} |
{\includegraphics*[width=0.44\linewidth]{\fpath/run_4yr_ADJtheta_arc_lev1_tim218_cmax5.0E+01.eps}} |
| 663 |
% |
% |
| 664 |
\subfigure[{\footnotesize |
\subfigure[{\footnotesize |
| 665 |
-48 months}] |
-48 months}] |
| 666 |
{\includegraphics*[width=0.44\linewidth]{figs/run_4yr_ADJtheta_arc_lev1_tim292_cmax5.0E+01.eps}} |
{\includegraphics*[width=0.44\linewidth]{\fpath/run_4yr_ADJtheta_arc_lev1_tim292_cmax5.0E+01.eps}} |
| 667 |
} |
} |
| 668 |
\caption{Same as Fig. XXX but for sea surface temperature |
\caption{Same as \reffig{4yradjheff} but for sea surface temperature |
| 669 |
\label{fig:4yradjthetalev1}} |
\label{fig:4yradjthetalev1}} |
| 670 |
\end{figure} |
\end{figure} |
| 671 |
|
|
| 690 |
|
|
| 691 |
\paragraph{Acknowledgements} |
\paragraph{Acknowledgements} |
| 692 |
We thank Jinlun Zhang for providing the original B-grid code and many |
We thank Jinlun Zhang for providing the original B-grid code and many |
| 693 |
helpful discussions. |
helpful discussions. ML thanks Elizabeth Hunke for multiple explanations. |
| 694 |
|
|
| 695 |
\bibliography{bib/journal_abrvs,bib/seaice,bib/genocean,bib/maths,bib/mitgcmuv,bib/fram} |
\bibliography{bib/journal_abrvs,bib/seaice,bib/genocean,bib/maths,bib/mitgcmuv,bib/fram} |
| 696 |
|
|
| 700 |
%%% mode: latex |
%%% mode: latex |
| 701 |
%%% TeX-master: t |
%%% TeX-master: t |
| 702 |
%%% End: |
%%% End: |
| 703 |
|
|
| 704 |
|
|
| 705 |
|
A Dynamic-Thermodynamic Sea ice Model for Ocean Climate |
| 706 |
|
Estimation on an Arakawa C-Grid |
| 707 |
|
|
| 708 |
|
Introduction |
| 709 |
|
|
| 710 |
|
Ice Model: |
| 711 |
|
Dynamics formulation. |
| 712 |
|
B-C, LSR, EVP, no-slip, slip |
| 713 |
|
parallellization |
| 714 |
|
Thermodynamics formulation. |
| 715 |
|
0-layer Hibler salinity + snow |
| 716 |
|
3-layer Winton |
| 717 |
|
|
| 718 |
|
Idealized tests |
| 719 |
|
Funnel Experiments |
| 720 |
|
Downstream Island tests |
| 721 |
|
B-grid LSR no-slip |
| 722 |
|
C-grid LSR no-slip |
| 723 |
|
C-grid LSR slip |
| 724 |
|
C-grid EVP no-slip |
| 725 |
|
C-grid EVP slip |
| 726 |
|
|
| 727 |
|
Arctic Setup |
| 728 |
|
Configuration |
| 729 |
|
OBCS from cube |
| 730 |
|
forcing |
| 731 |
|
1/2 and full resolution |
| 732 |
|
with a few JFM figs from C-grid LSR no slip |
| 733 |
|
ice transport through Canadian Archipelago |
| 734 |
|
thickness distribution |
| 735 |
|
ice velocity and transport |
| 736 |
|
|
| 737 |
|
Arctic forward sensitivity experiments |
| 738 |
|
B-grid LSR no-slip |
| 739 |
|
C-grid LSR no-slip |
| 740 |
|
C-grid LSR slip |
| 741 |
|
C-grid EVP no-slip |
| 742 |
|
C-grid EVP slip |
| 743 |
|
C-grid LSR no-slip + Winton |
| 744 |
|
speed-performance-accuracy (small) |
| 745 |
|
ice transport through Canadian Archipelago differences |
| 746 |
|
thickness distribution differences |
| 747 |
|
ice velocity and transport differences |
| 748 |
|
|
| 749 |
|
Adjoint sensitivity experiment on 1/2-res setup |
| 750 |
|
Sensitivity of sea ice volume flow through Fram Strait |
| 751 |
|
*** Sensitivity of sea ice volume flow through Canadian Archipelago |
| 752 |
|
|
| 753 |
|
Summary and conluding remarks |
Parent Directory
| 
Revision Log
| 
Revision Graph
| 
Patch