# Estimation of Open Boundary Conditions Based on an Isopycnic-Coordinate Internal Tidal Model with Adjoint Assimilation Method

Authors: Yanqiu Gao; Guangzhen Jin; Haibo Chen; Xianqing Lv
Journal: Mathematical Problems in Engineering
 (2013)

---

## Abstract

The isopycnic-coordinate internal tidal model with adjoint assimilation method is developed into a three-layer version. Two groups of ideal experiments are carried out in order to investigate the estimation of spatially varying open boundary conditions (OBCs). In group 1, different independent point schemes (IPSs) are used to invert 6 kinds of prescribed distributions of OBCs. Results show that, after assimilation, the cost functions and their gradient norms are reduced by about 2 orders of magnitude and by about 1 order of magnitude, respectively; the mean absolute errors (MAEs) in OBCs and the vector differences of horizontal current are reduced by 1 order of magnitude and by more than 23.53% compared with the values before assimilation, respectively. The results demonstrate that the three-layered model has a good ability in estimating the spatially varying OBCs; the use of IPSs can effectively improve the estimation precision; fewer independent points are feasible when the distribution is simpler, and appropriately more independent points are required when the distribution is more complex. In group 2, by using the optimal IPSs in group 1, the model is also able to successfully invert the OBCs on a real topography in the Luzon Strait area. The results are important to the study of the internal tide in the South China Sea.


---

## Body

## 1. Introduction

Internal tides play an important role in the dissipation of barotropic tides and in mixing the deep ocean [1, 2] and consequently contribute to the large-scale ocean circulation. Numerical simulation is an effective mean to study internal tides [3, 4]. The open boundary conditions (OBCs) are very important in regional ocean modeling and have a critical impact on the simulation results of the internal tides. The treatment of OBCs has been investigated by many researchers [5–10].The adjoint assimilation method has been widely used for parameter estimation. Its basic idea is as follows. First, a forward model is defined by a series of control equations and several initial and boundary conditions. Second, an adjoint model is derived based on the forward model. The cost function keeps decreasing while the parameters are optimized by the adjoint model, achieving an organic combination of numerical model and observations. By assimilating the data in the interior region, the adjoint assimilation method can optimize the parameters automatically and reach a global optimization. In a 3D shallow water flow system (TRIWAQ), Heemink et al. [11] used the adjoint approach to estimate the harmonic constants at the open boundaries, the friction parameter, the viscosity parameter, and the water depth by assimilating the tide gauge data as well as the altimeter data. Taillandier et al. [12] developed a four-dimensional variational data-assimilation approach to estimate the initial and open boundary conditions that force a coastal model according to interior observations. Ayoub [13] carried out experiments to estimate the OBCs in a North Atlantic circulation model using the adjoint method. Gejadze and Copeland [14] studied the open boundary control problem for free-surface barotropic Navier-Stokes equations with adjoint data assimilation technology. Zhang and Lu [15] optimized the OBCs and simulated the M2 tide and tidal current in the Bohai and Yellow Seas by assimilating satellite altimetry. Kazantsev [16] used a variational data assimilation technique to control boundary conditions on rigid boundaries for a shallow-water model. Chen et al. [17] constructed an adjoint assimilation model for simulation of internal tides and simply tested with a series of ideal experiments in which several prescribed spatial distributions of OBCs were successfully inverted by assimilating the model-generated pseudoobservations. Cao et al. [18] and Guo et al. [19] studied the inversion of OBCs by assimilating the T/P altimeter data with adjoint assimilation method in a 2D tidal model. Cao et al. [20] optimized OBCs in a 3D internal tidal model with the adjoint method around Hawaii.In this paper, the isopycnic-coordinate internal tidal model constructed by Chen et al. [17] is developed into a three-layer version in order to investigate the estimation of spatially varying OBCs. Two groups of ideal experiments are carried out. One group is used to verify the efficiency of the model on an ideal topography and to test the ability of the model in estimation of the spatially varying OBCs. In the other group, we estimate the OBCs on a real topography in the Luzon Strait area using the optimal independent point schemes (IPSs) obtained in the former group. This paper is organized as follows. The isopycnic-coordinate internal tidal model with adjoint assimilation method, the optimization method, and the experimental design are briefly introduced in Section 2. The experimental results are discussed in detail in Section 3. Finally, we make a summary and draw some conclusions in Section 4.


## 2. Model and Experimental Design

The isopycnic-coordinate numerical model for simulation of internal tides with adjoint assimilation method is composed of a forward model and its adjoint model. The forward model is used to simulate the internal tides and the adjoint model is used to optimize the modal control variables. Based on this assimilation model, two groups of ideal experiments are carried out to investigate the estimation of spatially varying OBCs.


### 2.1. Isopycnic-Coordinate Internal Tidal Model with Adjoint Assimilation Method

The isopycnic-coordinate internal tidal model employs spherical coordinates in horizontal direction and isopycnic coordinates in the vertical one. This model is based on the nonlinear, time-dependent, and free-surface hydrostatic equations, and is subject to the hydrostatic approximations. The internal mode equations are as follows:upper layer:(1)∂q1∂t+1Rcosφ∂(q1u1)∂λ+1Rcosφ∂(q1v1cosφ)∂φ=0,(2)∂u1∂t+u1Rcosφ∂u1∂λ+v1R∂u1∂φ-u1v1tanφR-fv1+gRcosφ∂∂λ∑m=13(qmρm-hm)-Ah1Δu1+Av1ρ-1q1h-1(u1-u2)=0,(3)∂v1∂t+u1Rcosφ∂v1∂λ+v1R∂v1∂φ+u12tanφR+fu1+gR∂∂φ∑m=13(qmρm-hm)-Ah1Δv1+Av1ρ-1q1h-1(v1-v2)=0,the second layer:(4)∂q2∂t+1Rcosφ∂(q2u2)∂λ+1Rcosφ∂(q2v2cosφ)∂φ=0,(5)∂u2∂t+u2Rcosφ∂u2∂λ+v2R∂u2∂φ-u2v2tanφR-fv2+gRcosφ∂∂λ∑m=13(qmρ2-hm)-Ah2Δu2-Av1ρ-1q2h-1(u1-u2)+Av2ρ-2q2h-2(u2-u3)+gRcosφ∑m=12(1ρ2-1ρm)∂qm∂φ=0,(6)∂v2∂t+u2Rcosφ∂v2∂λ+v2R∂v2∂φ+u22tanφR+fu2+gR∂∂φ∑k=13(qkρ2-hk)-Ah2Δv2-Av1ρ-1q2h-1(v1-v2)+Av2ρ-2q2h-2(v2-v3)+gR∑m=12(1ρ2-1ρm)∂qm∂φ=0,lower layer:(7)∂q3∂t+1Rcosφ∂(q3u3)∂λ+1Rcosφ∂(q3v3cosφ)∂φ=0,(8)∂u3∂t+u3Rcosφ∂u3∂λ+v3R∂u3∂φ-u3v3tanφR-fv3+gRcosφ∂∂λ∑m=13(qmρ3-hm)-Ah3Δu3-Av2ρ-2q3h-2(u2-u3)+κρ3u3q3(u3)2+(v3)2=0,(9)∂v3∂t+u3Rcosφ∂v3∂λ+v3R∂v3∂φ+u32tanφR+fu3+gR∂∂φ∑m=13(qmρ3-hm)-Ah3Δv3-Av2ρ-2q3h-2(v2-v3)+κρ3v3q3(u3)2+(v3)2=0.
Here, t is the time, λ and φ are the east longitude and north latitude, respectively,  um(λ,φ,t) and vm(λ,φ,t) are horizontal velocities in λ and φ in the mth (m=1,2,3, surface for m=1) layer, respectively, and  qm(λ,φ,t) is the time-varying layer mass and
(10)q1=ρ1(h1+η1-η2),q2=ρ2(h2+η2-η3),q3=ρ3(h3+η3),
where hm and ηm are the undisturbed layer thickness and interface elevation above the undisturbed level, respectively. R is the radius of the earth, g is the gravitational acceleration, f is the Coriolis parameter, Ahm is the horizontal eddy viscosity coefficient in the mth layer, respectively, Δ is the Laplace operator, and
(11)Δ(u,v)=1R2cos2φ∂2(u,v)∂λ2+1R2cosφ∂∂φ[cosφ∂(u,v)∂φ].Avm(m=1,2) and κ are the interface and bottom friction coefficients, respectively.  ρ-m=(ρm+ρm+1)/2 and h-m=(hm+hm+1)/2(m=1,2), where ρm is the potential density which is constant in each layer. In the forward model, q, u, and v are the main output and called the state variables in this paper.The barotropic currents are defined by(12)U=(q1u1+q2u2+q3u3)Q,V=(q1v1+q2v2+q3v3)Q,whereQ=q1+q2+q3.Integrating the continuity equations and the momentum equations in the vertical, we get the external mode as follows:(13)∂Q∂t+1Rcosφ∂(QU)∂λ+1Rcosφ∂(QVcosφ)∂φ=0,(14)∂U∂t-fV+Mλ+RU+κρ3u3Q(u3)2+(v3)2=0,(15)∂V∂t+fU+Mφ+RV+κρ3v3Q(u3)2+(v3)2=0,
where
(16)Mλ=gRcosφ∂∂λ(Qρ-h)+gQRcosφ∑k=23[qk∑m=12(1ρk-1ρm)∂qm∂λ],Mφ=gR∂∂φ(Qρ-h)+gQR∑k=23[qk∑m=12(1ρk-1ρm)∂qm∂φ],RU=1Q∑m=13qm(umRcosφ∂um∂λ+vmR∂um∂φ-umvmtanφR)-1Q∑m=13qmAhmΔum-∑m=13um∂∂t(qmQ),RV=1Q∑m=13qm(umRcosφ∂vm∂λ+vmR∂vm∂φ+um2tanφR)-1Q∑m=13qmAhmΔvm-∑m=13vm∂∂t(qmQ).Hereh(λ,φ)  (=h1+h2+h3) is the undisturbed water depth; ρ(λ,φ,t) is the vertically averaged density and ρ=Q/(h+η1). Note that the second terms of Mλ and Mφ represent the resultant effect of the horizontal density gradients in all layers on the sea surface and are the main causes of the surface manifestation of internal tides.Based on the governing equations (1)–(9) of the forward model, its adjoint model can be constructed. Similar to the forward model, the adjoint model also consists of the internal and external modes. The details of the adjoint model in a generalized form can be found in the work of Chen et al. [17].The cost function measuring the data misfit between the model output and observations is minimized through optimizing the control parameters. The cost function is defined as(17)J(q,u,v;p)=12∫Σ[Ku∑m=13(um-u^m)2+Kv∑m=13(vm-v^m)2]dσ,
where m(m=1,2,3) is layer index, ∑ denotes the computing area of both time and space, p represents the generalized control parameters, um and vm are the model simulated, and u^m and v^m are the observations. Here Ku and Kv are the weight matrices and theoretically should be the inverse of the observation error covariance matrix. In this paper, since only ideal experiments with model-generated pseudoobservations are performed, by assuming that the errors in the data are uncorrelated and equally weighted, and Ku and Kv are reduced to unit matrices [21].The Lagrangian function is defined by(18)L(q,u,v;qa,ua,va;p)=J(q,u,v;p)+∑m=13∫Σ[qam·Eq.(3m-2)+uamqm·Eq.(3m-1)+vamqm·Eq.(3m)]dσ,
where qam, uam, and vam are called adjoint variables of state variables qm, um, and vm, respectively.According to the typical theory of Lagrangian multiplier method, we have the following first-order derivates of Lagrangian function with respect to all the variables and parameters:(19)∂L∂qam=0,∂L∂uam=0,∂L∂vam=0;(20)∂L∂qm=0,∂L∂um=0,∂L∂vm=0;(21)∂L∂p=0.
Equation (19) returns the governing equations (1)–(9). The adjoint equations can be derived from (20). From (21) we can obtain the gradients of the cost function with respect to control parameters.In the model, OBC at thenth time step along the open boundary at a certain open boundary grid point (I,J) is described as
(22)ηI,Jn=aI,Jcos(ωnΔt)+bI,Jsin(ωnΔt),
where η is water level, ω denotes the frequency of M2 constituent, Δt is the time step length, and aI,J and bI,J are the Fourier coefficients as well as the tunable parameters of the model. The gradients of cost function with respect to aI,J and bI,J can be deduced from (21) which yields
(23)∂J∂aI,J+∑nTI,Jncos(ωnΔt)=0,∂J∂bI,J+∑nTI,Jnsin(ωnΔt)=0,
where TI,Jn refers to Chen et al. [17].So the inverse problem about the OBCs is transformed into the one about Fourier coefficientsa and b.


### 2.2. Experimental Design

Based on the three-layered assimilation model, two groups of ideal experiments (i.e., groups 1 and 2) are carried out to investigate the influence of IPS on the OBCs estimation. In order to test the ability in estimating the spatially varying OBCs in a nutshell, onlya (one of Fourier coefficients) along east boundary is estimated in this paper. Meanwhile, other open boundary conditions are set to be known.In group 1, on an ideal topography, six kinds of spatial distributions of OBCs are inverted with this three-layered model. Different IPSs are employed in accordance with the spatial distributions of OBCs. The optimal schemes are selected after comparing the results of group 1 and are used to estimate the OBCs in group 2 in which a real topography in the Luzon Strait is employed.The calculation regions (from116°E to 124°E and from 18°N to 23°N) of groups 1 and 2 are shown in Figures 1 and 2, respectively. In both figures, the white stars denote the 15 observation locations and the white dots denote the 27 open boundary points where the barotropic tidal force is installed. The horizontal resolution in this model is 10′×10′, and there are 49×31 grids totally in this area. The two undisturbed interfaces of this three-layered model are placed at the depth of 200 m and 400 m, respectively, and the potential densities in three layers are 1023.42 kg/m3, 1025.02 kg/m3, and 1026.77 kg/m3, respectively, based on the temperature-salinity data of the world oceans data sets (WOA05). The Coriolis parameter is set to be local value, the angular frequency of M2 tide is 1.4050789025×10-4s-1, and the whole-time step is 496.863 s (1/90 of the period of M2 tide) for both external and internal modes. The time of simulation is equal to 20 periods of M2 tide. The horizontal eddy viscosity coefficient is chosen to be Ah=1000. The coefficients of bottom and interface friction are taken as κ=0.003 and Av=0.005, respectively. Every observation location has “observations” of tidal current, if possible, in the upper two layers. The Flather conditions are installed along four boundaries and, without loss of generality, the unknown open boundary is supposed to be east boundary.Figure 1
Plan view of calculation region of group 1.Figure 2
Plan view of calculation region of group 2.The 6 kinds of prescribed distributions of OBCs in group 1 are shown in Figure3. First of all, direct inversion without IPS for each kind of distribution is carried out as a comparison experiment. Next, different number of independent points (Table 1) is employed according to the complexity of OBC distributions. The OBCs at each point along the east boundary can be derived from linear interpolation (the Cressman interpolation is employed in this paper) of values at independent points.Table 1
The distributions of OBCs and the IPSs.
Distributions
The number of independent points
1
2
—
—
2
2
4
—
3
3
5
—
4
5
9
—
5
4
7
13
6
6
9
17Figure 3
Six kinds of prescribed OBCs (the black solid lines) and their inverted distributions.For each distribution there is an optimal IPS in group 1. The optimal IPSs are selected to invert the 6 kinds of prescribed distributions of OBCs in group 2.


## 2.1. Isopycnic-Coordinate Internal Tidal Model with Adjoint Assimilation Method

The isopycnic-coordinate internal tidal model employs spherical coordinates in horizontal direction and isopycnic coordinates in the vertical one. This model is based on the nonlinear, time-dependent, and free-surface hydrostatic equations, and is subject to the hydrostatic approximations. The internal mode equations are as follows:upper layer:(1)∂q1∂t+1Rcosφ∂(q1u1)∂λ+1Rcosφ∂(q1v1cosφ)∂φ=0,(2)∂u1∂t+u1Rcosφ∂u1∂λ+v1R∂u1∂φ-u1v1tanφR-fv1+gRcosφ∂∂λ∑m=13(qmρm-hm)-Ah1Δu1+Av1ρ-1q1h-1(u1-u2)=0,(3)∂v1∂t+u1Rcosφ∂v1∂λ+v1R∂v1∂φ+u12tanφR+fu1+gR∂∂φ∑m=13(qmρm-hm)-Ah1Δv1+Av1ρ-1q1h-1(v1-v2)=0,the second layer:(4)∂q2∂t+1Rcosφ∂(q2u2)∂λ+1Rcosφ∂(q2v2cosφ)∂φ=0,(5)∂u2∂t+u2Rcosφ∂u2∂λ+v2R∂u2∂φ-u2v2tanφR-fv2+gRcosφ∂∂λ∑m=13(qmρ2-hm)-Ah2Δu2-Av1ρ-1q2h-1(u1-u2)+Av2ρ-2q2h-2(u2-u3)+gRcosφ∑m=12(1ρ2-1ρm)∂qm∂φ=0,(6)∂v2∂t+u2Rcosφ∂v2∂λ+v2R∂v2∂φ+u22tanφR+fu2+gR∂∂φ∑k=13(qkρ2-hk)-Ah2Δv2-Av1ρ-1q2h-1(v1-v2)+Av2ρ-2q2h-2(v2-v3)+gR∑m=12(1ρ2-1ρm)∂qm∂φ=0,lower layer:(7)∂q3∂t+1Rcosφ∂(q3u3)∂λ+1Rcosφ∂(q3v3cosφ)∂φ=0,(8)∂u3∂t+u3Rcosφ∂u3∂λ+v3R∂u3∂φ-u3v3tanφR-fv3+gRcosφ∂∂λ∑m=13(qmρ3-hm)-Ah3Δu3-Av2ρ-2q3h-2(u2-u3)+κρ3u3q3(u3)2+(v3)2=0,(9)∂v3∂t+u3Rcosφ∂v3∂λ+v3R∂v3∂φ+u32tanφR+fu3+gR∂∂φ∑m=13(qmρ3-hm)-Ah3Δv3-Av2ρ-2q3h-2(v2-v3)+κρ3v3q3(u3)2+(v3)2=0.
Here, t is the time, λ and φ are the east longitude and north latitude, respectively,  um(λ,φ,t) and vm(λ,φ,t) are horizontal velocities in λ and φ in the mth (m=1,2,3, surface for m=1) layer, respectively, and  qm(λ,φ,t) is the time-varying layer mass and
(10)q1=ρ1(h1+η1-η2),q2=ρ2(h2+η2-η3),q3=ρ3(h3+η3),
where hm and ηm are the undisturbed layer thickness and interface elevation above the undisturbed level, respectively. R is the radius of the earth, g is the gravitational acceleration, f is the Coriolis parameter, Ahm is the horizontal eddy viscosity coefficient in the mth layer, respectively, Δ is the Laplace operator, and
(11)Δ(u,v)=1R2cos2φ∂2(u,v)∂λ2+1R2cosφ∂∂φ[cosφ∂(u,v)∂φ].Avm(m=1,2) and κ are the interface and bottom friction coefficients, respectively.  ρ-m=(ρm+ρm+1)/2 and h-m=(hm+hm+1)/2(m=1,2), where ρm is the potential density which is constant in each layer. In the forward model, q, u, and v are the main output and called the state variables in this paper.The barotropic currents are defined by(12)U=(q1u1+q2u2+q3u3)Q,V=(q1v1+q2v2+q3v3)Q,whereQ=q1+q2+q3.Integrating the continuity equations and the momentum equations in the vertical, we get the external mode as follows:(13)∂Q∂t+1Rcosφ∂(QU)∂λ+1Rcosφ∂(QVcosφ)∂φ=0,(14)∂U∂t-fV+Mλ+RU+κρ3u3Q(u3)2+(v3)2=0,(15)∂V∂t+fU+Mφ+RV+κρ3v3Q(u3)2+(v3)2=0,
where
(16)Mλ=gRcosφ∂∂λ(Qρ-h)+gQRcosφ∑k=23[qk∑m=12(1ρk-1ρm)∂qm∂λ],Mφ=gR∂∂φ(Qρ-h)+gQR∑k=23[qk∑m=12(1ρk-1ρm)∂qm∂φ],RU=1Q∑m=13qm(umRcosφ∂um∂λ+vmR∂um∂φ-umvmtanφR)-1Q∑m=13qmAhmΔum-∑m=13um∂∂t(qmQ),RV=1Q∑m=13qm(umRcosφ∂vm∂λ+vmR∂vm∂φ+um2tanφR)-1Q∑m=13qmAhmΔvm-∑m=13vm∂∂t(qmQ).Hereh(λ,φ)  (=h1+h2+h3) is the undisturbed water depth; ρ(λ,φ,t) is the vertically averaged density and ρ=Q/(h+η1). Note that the second terms of Mλ and Mφ represent the resultant effect of the horizontal density gradients in all layers on the sea surface and are the main causes of the surface manifestation of internal tides.Based on the governing equations (1)–(9) of the forward model, its adjoint model can be constructed. Similar to the forward model, the adjoint model also consists of the internal and external modes. The details of the adjoint model in a generalized form can be found in the work of Chen et al. [17].The cost function measuring the data misfit between the model output and observations is minimized through optimizing the control parameters. The cost function is defined as(17)J(q,u,v;p)=12∫Σ[Ku∑m=13(um-u^m)2+Kv∑m=13(vm-v^m)2]dσ,
where m(m=1,2,3) is layer index, ∑ denotes the computing area of both time and space, p represents the generalized control parameters, um and vm are the model simulated, and u^m and v^m are the observations. Here Ku and Kv are the weight matrices and theoretically should be the inverse of the observation error covariance matrix. In this paper, since only ideal experiments with model-generated pseudoobservations are performed, by assuming that the errors in the data are uncorrelated and equally weighted, and Ku and Kv are reduced to unit matrices [21].The Lagrangian function is defined by(18)L(q,u,v;qa,ua,va;p)=J(q,u,v;p)+∑m=13∫Σ[qam·Eq.(3m-2)+uamqm·Eq.(3m-1)+vamqm·Eq.(3m)]dσ,
where qam, uam, and vam are called adjoint variables of state variables qm, um, and vm, respectively.According to the typical theory of Lagrangian multiplier method, we have the following first-order derivates of Lagrangian function with respect to all the variables and parameters:(19)∂L∂qam=0,∂L∂uam=0,∂L∂vam=0;(20)∂L∂qm=0,∂L∂um=0,∂L∂vm=0;(21)∂L∂p=0.
Equation (19) returns the governing equations (1)–(9). The adjoint equations can be derived from (20). From (21) we can obtain the gradients of the cost function with respect to control parameters.In the model, OBC at thenth time step along the open boundary at a certain open boundary grid point (I,J) is described as
(22)ηI,Jn=aI,Jcos(ωnΔt)+bI,Jsin(ωnΔt),
where η is water level, ω denotes the frequency of M2 constituent, Δt is the time step length, and aI,J and bI,J are the Fourier coefficients as well as the tunable parameters of the model. The gradients of cost function with respect to aI,J and bI,J can be deduced from (21) which yields
(23)∂J∂aI,J+∑nTI,Jncos(ωnΔt)=0,∂J∂bI,J+∑nTI,Jnsin(ωnΔt)=0,
where TI,Jn refers to Chen et al. [17].So the inverse problem about the OBCs is transformed into the one about Fourier coefficientsa and b.


## 2.2. Experimental Design

Based on the three-layered assimilation model, two groups of ideal experiments (i.e., groups 1 and 2) are carried out to investigate the influence of IPS on the OBCs estimation. In order to test the ability in estimating the spatially varying OBCs in a nutshell, onlya (one of Fourier coefficients) along east boundary is estimated in this paper. Meanwhile, other open boundary conditions are set to be known.In group 1, on an ideal topography, six kinds of spatial distributions of OBCs are inverted with this three-layered model. Different IPSs are employed in accordance with the spatial distributions of OBCs. The optimal schemes are selected after comparing the results of group 1 and are used to estimate the OBCs in group 2 in which a real topography in the Luzon Strait is employed.The calculation regions (from116°E to 124°E and from 18°N to 23°N) of groups 1 and 2 are shown in Figures 1 and 2, respectively. In both figures, the white stars denote the 15 observation locations and the white dots denote the 27 open boundary points where the barotropic tidal force is installed. The horizontal resolution in this model is 10′×10′, and there are 49×31 grids totally in this area. The two undisturbed interfaces of this three-layered model are placed at the depth of 200 m and 400 m, respectively, and the potential densities in three layers are 1023.42 kg/m3, 1025.02 kg/m3, and 1026.77 kg/m3, respectively, based on the temperature-salinity data of the world oceans data sets (WOA05). The Coriolis parameter is set to be local value, the angular frequency of M2 tide is 1.4050789025×10-4s-1, and the whole-time step is 496.863 s (1/90 of the period of M2 tide) for both external and internal modes. The time of simulation is equal to 20 periods of M2 tide. The horizontal eddy viscosity coefficient is chosen to be Ah=1000. The coefficients of bottom and interface friction are taken as κ=0.003 and Av=0.005, respectively. Every observation location has “observations” of tidal current, if possible, in the upper two layers. The Flather conditions are installed along four boundaries and, without loss of generality, the unknown open boundary is supposed to be east boundary.Figure 1
Plan view of calculation region of group 1.Figure 2
Plan view of calculation region of group 2.The 6 kinds of prescribed distributions of OBCs in group 1 are shown in Figure3. First of all, direct inversion without IPS for each kind of distribution is carried out as a comparison experiment. Next, different number of independent points (Table 1) is employed according to the complexity of OBC distributions. The OBCs at each point along the east boundary can be derived from linear interpolation (the Cressman interpolation is employed in this paper) of values at independent points.Table 1
The distributions of OBCs and the IPSs.
Distributions
The number of independent points
1
2
—
—
2
2
4
—
3
3
5
—
4
5
9
—
5
4
7
13
6
6
9
17Figure 3
Six kinds of prescribed OBCs (the black solid lines) and their inverted distributions.For each distribution there is an optimal IPS in group 1. The optimal IPSs are selected to invert the 6 kinds of prescribed distributions of OBCs in group 2.


## 3. Experimental Results and Discussions

The forward model is run with certain prescribed OBCs, and the model-generated results of the currents (i.e., currents in the upper two layers) at these observation locations are taken as the pseudo-observations. Then, an initial value (taken as zero here) ofa is assigned to run the forward model. The difference between the simulated values and pseudo-observations plays the role of the external force of the adjoint model. The optimized OBCs can be obtained through the backward integration of the adjoint equations. The inverse integral time of the adjoint equations is equal to a period of M2 tide. With the procedure above repeated, the OBCs will be optimized continually and the difference between simulated values and pseudo-observations will be diminished. Meanwhile, the difference between the prescribed and the estimated OBCs is also decreased. Different iteration numbers are tried in group 1, and 100 iterations are good choice based on the following considerations. (1) The misfit between “observations” and the simulated results is very small and approximately constant when it declines to a certain value after 100 iterations; (2) the cost functions are almost no longer falling after 100 iterations; (3) the calculation amount is acceptable for the 100 iterations. In group 2, the iteration number is also taken as 100.


### 3.1. Results and Discussions of Group 1

In this section, six kinds of spatial distributions of OBCs are inverted on an ideal topography. The results are shown in Figure3. The prescribed distributions are all inverted successfully. The inverted distributions obtained without IPS have some obvious fluctuations, while those obtained with IPSs show good consistency with the prescribed distributions. The results demonstrate that the model has a good ability in estimating the spatially varying OBCs. Moreover, the IPSs can further improve the estimation precision.The mean absolute errors (MAEs) between the prescribed OBCs and the values show that utilization of fewer independent points is feasible when the distribution of OBCs is simpler, and appropriately more independent points can lead to higher precision when the distribution of OBCs is more complex. The result that the use of IPSs can indeed effectively improve the precision of parameter estimation is consistent with the results in [22–24]. The MAEs between the prescribed OBCs and the values obtained with the optimal IPSs corresponding to the 6 kinds of spatial distributions are shown in Table 2. The MAEs in OBCs are reduced by 1 order of magnitude compared with the values before assimilation. The results indicate that the model performs great effect on estimating the OBCs with the optimal IPSs.Table 2
The MAEs in OBCs before assimilation (BA) and after assimilation (AA). The values in parentheses are the optimal numbers of independent points.
Distributions
1 (2)
2 (4)
3 (5)
4 (9)
5 (13)
6 (17)
BA
0.550
0.550
0.567
0.305
0.326
0.305
AA
0.017
0.033
0.029
0.020
0.023
0.058The cost function and theL2 norms of its gradient with respect to OBCs are two important criteria to evaluate the assimilation efficiency. The iteration histories of the cost functions and their gradient norms (all values are normalized by their own initial values) for the 6 kinds of distributions of OBCs are shown in Figures 4 and 5 (only the histories of the previous 50 iterations are plotted), respectively.Figure 4
The iteration histories of the cost functions for group 1.Figure 5
The iteration histories of theL2 norms of gradients of the cost functions with respect to OBCs for group 1.The cost functions are declining dramatically during the assimilation process. After 100 iterations, the cost functions and their gradient norms are reduced by about 2 orders of magnitude and by about 1 order of magnitude compared with their initial values, respectively. The declining cost functions and gradient norms indicate that the OBCs have been indeed optimized and are getting close to their prescribed values.Statistics of the vector differences in horizontal tidal current before and after assimilation are shown in Table3. From Table 3 we can see that the vector differences in horizontal tidal current are reduced by more than 23.53% compared with the values before assimilation. It indicates that all the 6 kinds of spatial distributions of OBCs have been inverted successfully.Table 3
Statistics of the errors before assimilation (BA) and after assimilation (AA). The values in parentheses in the first column are the optimal numbers of independent points and the unit of the errors is cms−1.
Distributions
u-vector difference
v-vector difference
BA
AA
1-AA/BA (%)
BA
AA
1-AA/BA (%)
1 (2)
2.65
0.19
92.83
0.92
0.17
81.52
2 (4)
2.72
0.40
85.29
0.85
0.60
29.41
3 (5)
1.45
0.07
95.17
0.32
0.12
62.50
4 (9)
0.43
0.19
55.81
0.41
0.11
83.17
5 (13)
0.37
0.05
86.49
0.18
0.06
66.67
6 (17)
0.27
0.20
25.93
0.17
0.13
23.53Different IPSs are employed to invert the 6 kinds of prescribed distributions of OBCs in group 1. After 100 iterations, the cost functions and their gradient norms are reduced by about 2 orders of magnitude, and by about 1 order of magnitude, respectively, compared with their initial values. Meanwhile, the MAEs in OBCs are reduced by about 1 order of magnitude, and the vector differences in horizontal tidal current are reduced by more than 23.53% compared with their values before assimilation. The results demonstrate that the model is effective and it has a good ability in estimating the spatially varying OBCs. Moreover, utilization of fewer independent points is feasible when the distribution of OBCs is simpler, and appropriately more independent points are required when the distribution of OBCs to be inverted is more complex. For each prescribed distribution of OBCs, the IPSs leading to the best result is taken as the optimal one. According to the results of group 1, the optimal numbers of independent points are 2, 4, 5, 9, 13, and 17 corresponding to the 6 kinds of prescribed OBC distributions, respectively.


### 3.2. Results and Discussions of Group 2

In this section, six experiments are carried out on a real topography in the Luzon Strait area. Corresponding to these 6 experiments, six kinds of prescribed distributions of OBCs are inverted, respectively, with the optimal IPSs derived in group 1. The results are shown in Figure6, which show that the inverted distributions are in good consistency with the prescribed distributions. And the results demonstrate that the 6 kinds of prescribed distributions are all inverted successfully.Figure 6
Six kinds of prescribed distributions of OBCs (the black solid lines) and their inverted distributions.The iteration histories of the cost functions and their gradient norms (all values are normalized by their own initial values) are shown in Figure7 (only the histories of previous 50 iterations are plotted). After 100 iterations, the cost functions and their gradient norms are both reduced by about 1 order of magnitude compared with their initial values. Table 4 shows the MAEs in OBCs before and after assimilation. The MAEs in OBCs are reduced by 1 order of magnitude compared with the values before assimilation. The experimental results indicate that the model can invert the 6 kinds of prescribed distributions of OBCs successfully on a real topography in the Luzon Strait area.Table 4
The MAEs in OBCs before assimilation (BA) and after assimilation (AA).
Distributions
1
2
3
4
5
6
BA
0.550
0.550
0.567
0.305
0.326
0.305
AA
0.029
0.031
0.043
0.025
0.018
0.061Figure 7
The iteration histories of the cost functions and their gradient norms with respect to OBCs for group 2.


## 3.1. Results and Discussions of Group 1

In this section, six kinds of spatial distributions of OBCs are inverted on an ideal topography. The results are shown in Figure3. The prescribed distributions are all inverted successfully. The inverted distributions obtained without IPS have some obvious fluctuations, while those obtained with IPSs show good consistency with the prescribed distributions. The results demonstrate that the model has a good ability in estimating the spatially varying OBCs. Moreover, the IPSs can further improve the estimation precision.The mean absolute errors (MAEs) between the prescribed OBCs and the values show that utilization of fewer independent points is feasible when the distribution of OBCs is simpler, and appropriately more independent points can lead to higher precision when the distribution of OBCs is more complex. The result that the use of IPSs can indeed effectively improve the precision of parameter estimation is consistent with the results in [22–24]. The MAEs between the prescribed OBCs and the values obtained with the optimal IPSs corresponding to the 6 kinds of spatial distributions are shown in Table 2. The MAEs in OBCs are reduced by 1 order of magnitude compared with the values before assimilation. The results indicate that the model performs great effect on estimating the OBCs with the optimal IPSs.Table 2
The MAEs in OBCs before assimilation (BA) and after assimilation (AA). The values in parentheses are the optimal numbers of independent points.
Distributions
1 (2)
2 (4)
3 (5)
4 (9)
5 (13)
6 (17)
BA
0.550
0.550
0.567
0.305
0.326
0.305
AA
0.017
0.033
0.029
0.020
0.023
0.058The cost function and theL2 norms of its gradient with respect to OBCs are two important criteria to evaluate the assimilation efficiency. The iteration histories of the cost functions and their gradient norms (all values are normalized by their own initial values) for the 6 kinds of distributions of OBCs are shown in Figures 4 and 5 (only the histories of the previous 50 iterations are plotted), respectively.Figure 4
The iteration histories of the cost functions for group 1.Figure 5
The iteration histories of theL2 norms of gradients of the cost functions with respect to OBCs for group 1.The cost functions are declining dramatically during the assimilation process. After 100 iterations, the cost functions and their gradient norms are reduced by about 2 orders of magnitude and by about 1 order of magnitude compared with their initial values, respectively. The declining cost functions and gradient norms indicate that the OBCs have been indeed optimized and are getting close to their prescribed values.Statistics of the vector differences in horizontal tidal current before and after assimilation are shown in Table3. From Table 3 we can see that the vector differences in horizontal tidal current are reduced by more than 23.53% compared with the values before assimilation. It indicates that all the 6 kinds of spatial distributions of OBCs have been inverted successfully.Table 3
Statistics of the errors before assimilation (BA) and after assimilation (AA). The values in parentheses in the first column are the optimal numbers of independent points and the unit of the errors is cms−1.
Distributions
u-vector difference
v-vector difference
BA
AA
1-AA/BA (%)
BA
AA
1-AA/BA (%)
1 (2)
2.65
0.19
92.83
0.92
0.17
81.52
2 (4)
2.72
0.40
85.29
0.85
0.60
29.41
3 (5)
1.45
0.07
95.17
0.32
0.12
62.50
4 (9)
0.43
0.19
55.81
0.41
0.11
83.17
5 (13)
0.37
0.05
86.49
0.18
0.06
66.67
6 (17)
0.27
0.20
25.93
0.17
0.13
23.53Different IPSs are employed to invert the 6 kinds of prescribed distributions of OBCs in group 1. After 100 iterations, the cost functions and their gradient norms are reduced by about 2 orders of magnitude, and by about 1 order of magnitude, respectively, compared with their initial values. Meanwhile, the MAEs in OBCs are reduced by about 1 order of magnitude, and the vector differences in horizontal tidal current are reduced by more than 23.53% compared with their values before assimilation. The results demonstrate that the model is effective and it has a good ability in estimating the spatially varying OBCs. Moreover, utilization of fewer independent points is feasible when the distribution of OBCs is simpler, and appropriately more independent points are required when the distribution of OBCs to be inverted is more complex. For each prescribed distribution of OBCs, the IPSs leading to the best result is taken as the optimal one. According to the results of group 1, the optimal numbers of independent points are 2, 4, 5, 9, 13, and 17 corresponding to the 6 kinds of prescribed OBC distributions, respectively.


## 3.2. Results and Discussions of Group 2

In this section, six experiments are carried out on a real topography in the Luzon Strait area. Corresponding to these 6 experiments, six kinds of prescribed distributions of OBCs are inverted, respectively, with the optimal IPSs derived in group 1. The results are shown in Figure6, which show that the inverted distributions are in good consistency with the prescribed distributions. And the results demonstrate that the 6 kinds of prescribed distributions are all inverted successfully.Figure 6
Six kinds of prescribed distributions of OBCs (the black solid lines) and their inverted distributions.The iteration histories of the cost functions and their gradient norms (all values are normalized by their own initial values) are shown in Figure7 (only the histories of previous 50 iterations are plotted). After 100 iterations, the cost functions and their gradient norms are both reduced by about 1 order of magnitude compared with their initial values. Table 4 shows the MAEs in OBCs before and after assimilation. The MAEs in OBCs are reduced by 1 order of magnitude compared with the values before assimilation. The experimental results indicate that the model can invert the 6 kinds of prescribed distributions of OBCs successfully on a real topography in the Luzon Strait area.Table 4
The MAEs in OBCs before assimilation (BA) and after assimilation (AA).
Distributions
1
2
3
4
5
6
BA
0.550
0.550
0.567
0.305
0.326
0.305
AA
0.029
0.031
0.043
0.025
0.018
0.061Figure 7
The iteration histories of the cost functions and their gradient norms with respect to OBCs for group 2.


## 4. Summary and Conclusions

The isopycnic-coordinate numerical model for simulation of internal tides with adjoint assimilation method constructed by Chen et al. [17] is developed into a three-layer version. Two groups of ideal experiments are carried out in order to investigate the estimation of spatially varying OBCs.In group 1, different IPSs are employed to invert 6 kinds of prescribed spatial distributions of OBCs on an ideal topography. After 100 iterations, the cost functions and their gradient norms are reduced by about 2 orders of magnitude and by about 1 order of magnitude, respectively, compared with their initial values. Meanwhile, compared with the values before assimilation, the MAEs in OBCs after assimilation are reduced by 1 order of magnitude, and the vector differences in horizontal tidal current after assimilation are reduced by more than 23.53%, respectively. The results of group 1 demonstrate that the inversion test is successful, and this model has a good ability in estimating the spatially varying OBCs. Meanwhile, the use of IPSs can indeed effectively improve the estimation precision of OBCs. Moreover, utilization of fewer independent points is feasible when the distribution of OBCs is simpler, and appropriately more independent points are required when the distribution of OBCs is more complex.In group 2, the optimal IPSs selected from group 1 are used to invert the 6 kinds of prescribed distributions of OBCs, and the real topography in the Luzon Strait area is installed. The cost functions and their gradient norms after 100 iterations are both reduced by about 1 order of magnitude compared with their initial values. The MAEs in OBCs are reduced by 1 order of magnitude compared with the values before assimilation. The results of group 2 demonstrate that the model also has a good ability in estimating the spatially varying OBCs on a real topography, which lays a foundation for the coming real experiments.In conclusion, the work in this paper shows us that the IPSs can improve the estimation of OBCs which play an important role in this model. Thus, the work is significant to the simulation of internal tides in the South China Sea. In addition, some particular choices, such as the nonuniformly distributed independent points and the weighted coefficients of the interpolation, can also be employed by other IPSs in the work of parameter estimation. However, the performance of IPS is likely improved in a practical application, further highlighting its advantage.



---

*Source: 321387-2013-06-12.xml*