Adjoint Data Assimilation of the Flow on the Southern Flank of Georges Bank: March–June 1999

: An adjoint data assimilation method was incorporated into the ECOM-si coastal ocean circulation model and applied to assimilate the ﬂow ﬁeld on the southern ﬂank of Georges Bank from March to June 1999. The model was driven by tidal forcing consisting of ten tidal constituents at the open boundary, observed winds, and surface heat ﬂuxes. Numerical experiments were conducted following a strategy to understand the critical issues affecting the efﬁciency and accuracy of the assimilated ﬂow ﬁeld. The adjoint data assimilation method signiﬁcantly improved the computational accuracy of subtidal currents, especially for the along-isobath velocity. The integration window length and iteration number were two parameters affecting the assimilation convergence rate toward the observations. In such a nonlinear dynamical system, using a window length close to the M 2 tidal period could make the adjoint model difﬁcult to converge, no matter how many iterations were made. Reducing the time interval for the window length speeded up the convergence rate, but it was paid out with the sacriﬁce of statistical conﬁdence. This assimilation experiment used a 6-h window length, which led to a faster convergence rate in the ﬁrst ten iterations. The best-assimilated ﬁelds that satisﬁed the error criteria were obtained as the iteration number increased.


Introduction
One of the primary goals of developing a numerical model is to simulate and predict the temporospatial variability of temperature, salinity, and currents in the ocean.The four-dimensional (4-D) data assimilation is aimed at integrating observational data into a model for improving the simulation systematically.A data assimilation system, such as nudging, optimal interpolation, an adjoint model, and Kalman filters, consists of three components: an observational dataset, a hydrodynamical model, and an assimilation scheme.Since the data have measurement uncertainty errors and the model dynamics and parameters are imperfect, a well-constructed melding algorithm should enable to improve the match between observed and simulated variables within an error range of observations and modeling [1] Therefore, once a hydrodynamical model is developed, an appropriate assimilation algorithm becomes more critical in ensuring the success of its application to the ocean.
Optimal interpolation (IO), ensemble Kalman filtering (EnKF), and variance or adjoint methods are the most popular algorithms used in data assimilation.The IO and EnKF are based on the estimation theory [2].They involve a sequence of estimation processes based on the minimization of the expected estimation errors in terms of the statistics for both the model and observations.The IO builds a linear combination of the observed and modelpredicted variables by calculating their error covariances [3].It is the most straightforward assimilation algorithm requiring a priori statistical assumption about the model noises and measurement uncertainties in observational data.The EnKF is the most sophisticated method, but ensemble simulations are computationally intensive to run [4][5][6][7][8].The variance

The Adjoint Data Assimilation Model
ECOM-si used in this study is a three-dimensional, free-surface, σ-coordinated vertical transformation coastal ocean circulation model [23].In ECOM-si, the Mellor and Yamada level 2.5 turbulent closure scheme was used for vertical viscosity parameterization [24].Unlike the time-splitting solver in POM, ECOM-si is integrated using a semiimplicit scheme, which makes the Fortran code simpler by removing explicit time smoothing [19].

The Adjoint Data Assimilation Model
ECOM-si used in this study is a three-dimensional, free-surface, σ-coordinated vertical transformation coastal ocean circulation model [23].In ECOM-si, the Mellor and Yamada level 2.5 turbulent closure scheme was used for vertical viscosity parameterization [24].Unlike the time-splitting solver in POM, ECOM-si is integrated using a semi-implicit scheme, which makes the Fortran code simpler by removing explicit time smoothing [19].
The adjoint ECOM-si model was conducted using a variational method.The primitive equations ECOM-si can be written in the vector form given as where x is a matrix array consisting of dependent variables such as the free-surface water level (ζ), water temperature (θ), water salinity (s) and the three-dimensional (3D) velocity components (u, v, and w).F is a nonlinear operator representing the ocean-governing equations, including the advective, Coriolis forcing, barotropic and baroclinic pressure gradient forcings, and diffusive terms in the momentum equations, as well as advective and diffusive terms for temperature and salinity equations.c is a matrix array containing the model parameters such as surface and bottom drag coefficients, light attenuation lengths, and open boundary and initial conditions.The model uses the Cartesian coordinates, in which x, y, and z are the cross-bank (northward: positive), along-bank (eastward: positive), and vertical axes of the orthogonal coordinate).The cost function in the adjoint model is defined as where x o and c o are the observed vectors for dependent variables and model parameters, respectively; and K and K c are validity coefficients, which were chosen as K = 0.5 for u and v, and K c was taken as 1 in our numerical experiments.Ω is the numerical computational domain.t represents time, and (t 1 ,t 2 ) is the time assimilation window.The cost function measures a "distance" (error) between observed and simulated quantities.A Lagrange function was defined as where λ is a Lagrange multiplier's matrix array for x (also denoted as adjoint variables), and the superscript "T" is the sign for the transposed matrix.The variational method used in the adjoint assimilation model was designed to determine x with minimum L with respect to λ, x, and c in terms of least-square fitting, i.e., This condition was equivalent to solving the Lagrange-Euler equations that satisfied the constraints in the form of Three equations can be derived from these three conditions of restraint.They are the Lagrange-Euler equations from the original (forward) model equations of adjoint equations of where Σ is the boundary of the computational domain.The gradient of the cost function can be derived by ∂L/∂c = 0.That is, Equation ( 8) can be rewritten as an iterative form as where c n is the previous nth iteration value of c, which is also equivalent to c o .n indicates the iteration number.
The forward model of ECOM-si was converted to the adjoint model using the tangent linear adjoint model compiler (OpenAD) developed by [22].The OpenAD is a numerical computer program that could automatically convert the differentiation equations to the adjoint code (https://www.mcs.anl.gov/OpenAD/,accessed on 24 November 2023).We used this software to automatically convert all subroutines in ECOM-si into adjoint codes.Then, the adjoint model was built by assembling all these subroutines based on the mathematical approach shown in Equations ( 6)- (9).A detailed description of OpenAD can be found at the website listed above.
A schematic of the adjoint assimilation approach is shown in Figure 2. First, the forward model of Equation ( 6) was run with measured wind forcing and surface heat flux in a tidal environment.Second, the model-simulated u and v within each integration window were output and compared with observed u and v at all available measurement positions.Then, the model-predicted and observed u and v differences were input into the adjoint model of Equation (7).Third, Equation ( 7) was run backward within the assimilation window, which provided a new set of λ.A new set of c (initial condition of u and v and the bottom drag coefficient) at n + 1 was then determined by solving Equation ( 9) with the input of λ from the adjoint model.A criterion based on the difference (ε) between c n+1 and c n was used to judge if the model provided an optimal flow field of u and v.If not, the simulation returned to the forward model with a new set of c as its initial condition and all three procedures were repeated.The iteration was conducted until the criterion was satisfied.This approach assumed that the inaccurate simulation results were due to the improper model setup in parameterization [25].In this study, the model parameters only considered the bottom stress coefficient since it was specified initially from an empirical formulation.However, in this case, a test showed that it required at least 40 iterations to meet the criterion of ε < 2.5 cm/s.For an assimilation window of 6 h, this iteration number produced the most rapid convergence of the assimilated currents to the observations.A detailed discussion on the accuracy of the assimilation regarding iteration number and assimilation window is given in Section 3.
where Σ is the boundary of the computational domain.The gradient of the cost function can be derived by   ⁄ = 0.That is, Equation ( 8) can be rewritten as an iterative form as where n c is the previous nth iteration value of c, which is also equivalent to o c .n indi- cates the iteration number.
The forward model of ECOM-si was converted to the adjoint model using the tangent linear adjoint model compiler (OpenAD) developed by [22].The OpenAD is a numerical computer program that could automatically convert the differentiation equations to the adjoint code (https://www.mcs.anl.gov/OpenAD/,accessed on 24 November 2023).We used this software to automatically convert all subroutines in ECOM-si into adjoint codes.Then, the adjoint model was built by assembling all these subroutines based on the mathematical approach shown in Equations ( 6)- (9).A detailed description of OpenAD can be found at the website listed above.
A schematic of the adjoint assimilation approach is shown in Figure 2. First, the forward model of Equation ( 6) was run with measured wind forcing and surface heat flux in a tidal environment.Second, the model-simulated u and v within each integration window were output and compared with observed u and v at all available measurement positions.Then, the model-predicted and observed u and v differences were input into the adjoint model of Equation (7).Third, Equation (7) was run backward within the assimilation window, which provided a new set of  .A new set of c (initial condition of u and v and the bottom drag coefficient) at n + 1 was then determined by solving Equation ( 9) with the input of  from the adjoint model.A criterion based on the difference (  ) between If not, the simulation returned to the forward model with a new set of c as its initial condition and all three procedures were repeated.The iteration was conducted until the criterion was satisfied.This approach assumed that the inaccurate simulation results were due to the improper model setup in parameterization [25].In this study, the model parameters only considered the bottom stress coefficient since it was specified initially from an empirical formulation.However, in this case, a test showed that it required at least 40 iterations to meet the criterion of  < 2.5 cm/s.For an assimilation window of 6 h, this iteration number produced the most rapid convergence of the assimilated currents to the observations.A detailed discussion on the accuracy of the assimilation regarding iteration number and assimilation window is given in Section 3. Adjoint assimilation modeling was conducted for a two-dimensional (2-D) numerical experiment on the assumption that the change of all variables in the along-isobath direction was negligible.The experiments used the 2-D version of ECOM-si developed by [26].This model was configured for a cross-bank section through the moored ADCP array, where the water depth varied from 100 to 35 m on the southern flank (Figure 1).[26] examined the dynamics attributed to tidal fronts over GB.They found that the tidal mixing front formation and along-frontal current were predominated by 2-D dynamical features, where the along-frontal barotropic and baroclinic pressure gradient could be ignored.The 2-D model, assuming no along-bank variations, was robust enough to simulate the stratified frontal currents observed on the south-north transect across the central bank (Figure 1).The tripod layout was on the cross-bank transect used in the 2-D experiments.This 2-D simplification made it easy to incorporate the observed data from the current meter array.It also helped avoid internal wave reflection at the northern open boundary.

Adjoint Data Assimilation Approach
The model used the σ-coordinate transformation defined by σ = (z − ζ)/(H + ζ) where σ varied from 0 at the surface to −1 at the bottom, and H was the mean sea level depth.A total of 50 uniform σ-layers were specified in the vertical, with a resolution of 2 m at the 100 m isobath and 0.7 m at the 33.4 m isobath.The cross-bank domain was configured with a nonuniform horizontal grid.The horizontal resolution (∆x) was 471 m near and across the bank area and increased linearly to 10 km over 33 grid points outside the domain of interest.Twenty additional grid points with a resolution of 10 km were added at the northern end to extend the northern domain 200 km longer.A gravity wave radiation condition was specified for the surface elevation and currents at the northern open boundary.The internal model time step was 110.4 s, 405 time steps over an M 2 (12.42 h.) tidal period.
Under a stratified condition, a free-surface model can generate surface and internal gravity waves.The wavelengths of internal waves were much shorter than surface gravity waves.Since the radiation condition could not transfer all wave energies from the boundary, reducing the horizontal resolution near the northern open boundary could help filter shortwavelength interval waves [26].
The model was driven by the tidal forcing at the southern open boundary and calculated wind stress plus heat flux at the sea surface.The tidal forcing consisted of ten tidal elevation constituents, including M 2 (12.42 h), S 2 (12 h), N 2 (12.66 h), L 2 (12.19 h), K 2 (11.97 h), MU 2 (12.87 h), K 1 (23.93 h), P 1 (24.07 h), O 1 (25.82 h), and MSF (14.77 days).The amplitudes and phases of these tidal constituents were specified by the output from the global 0.5 • × 0.5 • tidal assimilation model developed by [27].These values were adjusted to produce an optimal fit to tidal currents observed from the ADCP measurements.
At first, the model ran with tidal forcing only for three months after it was fully ramped up from the initial.After tidal currents in the model reached an equilibrium state, the calculated wind stress and surface heat flux were added at the correct tidal phase when ADCP data were recorded.The March-April climatologically hydrographic fields were used to set up the initial conditions of temperature and salinity (Figure 3).The climatological dataset was downloaded from the Dartmouth College Ocean Modeling Group website (http://www-nml.dartmouth.edu,accessed on 24 November 2023) [28,29].The data were processed by investigators at the Bedford Institute of Oceanography.The wind stress was calculated using the meteorological data recorded on the NOAA environmental buoy #44011 east of GB (Figures 1 and 4).Surface heat flux was calculated using an empirical formula based on the air-sea temperature difference recorded on buoy #44011.Relative humidity and cloud coverage were estimated using the data from the NOAA Climate Diagnostic Center (http://www.cdc.noaa.gov,accessed on 24 November 2023) [30,31].The amplitudes and phases of the observed and simulated tidal currents for each tidal constituent were calculated using a MATLAB R2014 harmonic analysis program developed by [32].To investigate processes in the coastal ocean using the adjoint method, the model was first run in the forward mode only for three months (from 19 March to 16 June 1999), then rerun with data assimilation.This strategy helped us first determine the level of accuracy for a 2-D model in predicting tidal and subtidal currents on GB and then evaluate if adding  To investigate processes in the coastal ocean using the adjoint method, the model was first run in the forward mode only for three months (from 19 March to 16 June 1999), then rerun with data assimilation.This strategy helped us first determine the level of accuracy for a 2-D model in predicting tidal and subtidal currents on GB and then evaluate if adding To investigate processes in the coastal ocean using the adjoint method, the model was first run in the forward mode only for three months (from 19 March to 16 June 1999), then rerun with data assimilation.This strategy helped us first determine the level of accuracy for a 2-D model in predicting tidal and subtidal currents on GB and then evaluate if adding the adjoint data assimilation method could improve the model accuracy.It also helped identify some key issues in conducting a data assimilation experiment in a coastal ocean environment with strong tidal currents.

Tidal Simulation
The time series of model-simulated currents at each grid point were harmonically reconstructed into the ten tidal constituents: M 2 , S 2 , N 2 , K 1 , P 1 , O 1 , L 2 , K 2 , MU 2 , and MSF.The same method was used for the observed currents recorded at the seven ADCP moorings.The model data comparisons were made for each tidal constituent and the time series of total tidal currents (the sum of ten tidal constituents).For example, station 567 was at the 62 m isobath.At this site, the vertical profile of the simulated M 2 tidal current was in good agreement with the observation (Figure 5), except for an underestimation in the minor axis in a ~25 m water column above the bottom and a slight overestimation in the major axis near a depth of ~35 m.Similar results were found in S 2 , N 2 , K 1 , P 1 , O 1 , L 2 , K 2 , MU 2 , and MSF, but their amplitudes were one order of magnitude smaller.The vertically averaged absolute values of the differences between model-simulated and observed currents were 0.68 and 1.05 cm/s in major and minor axes, 1.0 • and 1.17 • in orientation and phase for M 2 .Statistics for the other nine tidal constituents are given in Table 1.These differences were within the estimated uncertainty of the ADCP measurements at this site.However, the overall simulation errors in the along-and cross-bank tidal velocity components over the seven ADCP sites were 3.81 cm/s and 2.33 cm/s, respectively.These errors are about 1.0-1.5 cm/s over the estimated overall uncertainty of the ADCP current measurements.the adjoint data assimilation method could improve the model accuracy.It also helped identify some key issues in conducting a data assimilation experiment in a coastal ocean environment with strong tidal currents.

Tidal Simulation
The time series of model-simulated currents at each grid point were harmonically reconstructed into the ten tidal constituents: M2, S2, N2, K1, P1, O1, L2, K2, MU2, and MSF.The same method was used for the observed currents recorded at the seven ADCP moorings.The model data comparisons were made for each tidal constituent and the time series of total tidal currents (the sum of ten tidal constituents).For example, station 567 was at the 62 m isobath.At this site, the vertical profile of the simulated M2 tidal current was in good agreement with the observation (Figure 5), except for an underestimation in the minor axis in a ~25 m water column above the bottom and a slight overestimation in the major axis near a depth of ~35 m.Similar results were found in S2, N2, K1, P1, O1, L2, K2, MU2, and MSF, but their amplitudes were one order of magnitude smaller.The vertically averaged absolute values of the differences between model-simulated and observed currents were 0.68 and 1.05 cm/s in major and minor axes, 1.0° and 1.17° in orientation and phase for M2.Statistics for the other nine tidal constituents are given in Table 1.These differences were within the estimated uncertainty of the ADCP measurements at this site.However, the overall simulation errors in the along-and cross-bank tidal velocity components over the seven ADCP sites were 3.81 cm/s and 2.33 cm/s, respectively.These errors are about 1.0-1.5 cm/s over the estimated overall uncertainty of the ADCP current measurements.Table 1.Difference between model-predicted and observed tidal current ellipse parameters.

Depth(m)
∆U maj (cm/s) cross-bank total velocities approximately agreed with each other on both semidiurnal (shortterm) and weekly-to-monthly (long-term) time scales.However, the model overestimated the flow magnitude, especially in the upper 40 m, where wind forcing was more critical.The simulated and observed velocity difference was much larger than the tidal-induced errors, implying that such a difference was partially due to the inaccurate wind stress and heat flux forcing conditions.It was supported by a comparison between simulated and observed subtidal currents shown in Figure 7.The model significantly overestimated the low-frequency variation of the along-bank velocity, especially during the strong southward wind period in April.In addition, the model also underestimated the wind-induced, lowfrequency variation of the cross-bank velocity near the bottom.The observed velocity showed the time-coherent variation near the bottom and surface, but the model did not capture these features.
should be robust enough to reproduce the significant tidal and subtidal features observed in the natural ocean.For this reason, in addition to tidal forcing, the wind stress and heat flux were also considered in numerical experiments over 19 March-16 June 1999.The simulated total (tidal plus subtidal) and subtidal flows were compared with observed currents at measurement depths of seven moorings.
The forward model reasonably simulated the total velocity measured at the ADCP sites.An example is given in Figure 6, in which the model-predicted and observed alongand cross-bank total velocities approximately agreed with each other on both semidiurnal (short-term) and weekly-to-monthly (long-term) time scales.However, the model overestimated the flow magnitude, especially in the upper 40 m, where wind forcing was more critical.The simulated and observed velocity difference was much larger than the tidalinduced errors, implying that such a difference was partially due to the inaccurate wind stress and heat flux forcing conditions.It was supported by a comparison between simulated and observed subtidal currents shown in Figure 7.The model significantly overestimated the low-frequency variation of the along-bank velocity, especially during the strong southward wind period in April.In addition, the model also underestimated the windinduced, low-frequency variation of the cross-bank velocity near the bottom.The observed velocity showed the time-coherent variation near the bottom and surface, but the model did not capture these features.Despite these deficiencies, the forward model did capture the primary pattern of the low-frequency currents in response to the varying wind stress at ADCP measurement sites.It suggested that the model was robust enough to resolve the lowest order of the correct physics (tidal plus surface forcing) controlling current variability on the southern flank of GB during the simulation period.Meanwhile, the 2-D model did not resolve the basin-scale pressure field set up by the larger-scale wind stress [34].Although errors due to the basin-scale pressure gradient forcing and the wind stress spatial variability could cause the subtidal model-data differences, they were secondary to surface and boundary forcing conditions on the cross-bank moored ADCP transect.Despite these deficiencies, the forward model did capture the primary pattern of the low-frequency currents in response to the varying wind stress at ADCP measurement sites.It suggested that the model was robust enough to resolve the lowest order of the correct physics (tidal plus surface forcing) controlling current variability on the southern of GB during the simulation period.Meanwhile, the 2-D model did not resolve the basin-scale pressure field set up by the larger-scale wind stress [34].Although errors due to the basin-scale pressure gradient forcing and the wind stress spatial variability could cause the subtidal model-data differences, they were secondary to surface and boundary forcing conditions on the cross-bank moored ADCP transect.

Assimilation Results for Tidal and Subtidal Currents
Turning on adjoint assimilation significantly improved the simulation accuracy of the total and subtidal currents at the ADCP measurement sites.It was evident by comparing Figure 8 with Figure 6.This comparison showed that the assimilated along-bank velocity quickly converged to the observed velocity after 12 iterations within a 6 h integration window, even though no significant improvement was found for the cross-bank velocity.Without data assimilation, the overall model errors in total along-and cross-bank velocities were 7.56 and 9.91 cm/s at the ADCP measurement points.After the adjoint assimilation was applied, the overall uncertainties were reduced to 7.10 and 6.63 cm/s (Table 2).A significant improvement was also evident for the subtidal current, shown in Figure 9.In

Assimilation Results for Tidal and Subtidal Currents
Turning on adjoint assimilation significantly improved the simulation accuracy of the total and subtidal currents at the ADCP measurement sites.It was evident by comparing Figure 8 with Figure 6.This comparison showed that the assimilated along-bank velocity quickly converged to the observed velocity after 12 iterations within a 6 h integration window, even though no significant improvement was found for the cross-bank velocity.Without data assimilation, the overall model errors in total along-and cross-bank velocities were 7.56 and 9.91 cm/s at the ADCP measurement points.After the adjoint assimilation was applied, the overall uncertainties were reduced to 7.10 and 6.63 cm/s (Table 2).A significant improvement was also evident for the subtidal current, shown in Figure 9.In the case without assimilations, overall simulation errors in along-and cross-bank subtidal velocities were 3.75 and 7.58 cm/s, respectively, at the ADCP measurement points.These errors were reduced to 2.77 and 4.16 cm/s in the assimilation case.The overall improvement at these comparison sites was 26 and 45% in the cross-and along-bank directions.Note: ∆u and ∆v are the standard deviation of the along-and cross-bank current differences, respectively.
the case without assimilations, overall simulation errors in along-and cross-bank subtidal velocities were 3.75 and 7.58 cm/s, respectively, at the ADCP measurement points.These errors were reduced to 2.77 and 4.16 cm/s in the assimilation case.The overall improvement at these comparison sites was 26 and 45% in the cross-and along-bank directions.the case without assimilations, overall simulation errors in along-and cross-bank subtidal velocities were 3.75 and 7.58 cm/s, respectively, at the ADCP measurement points.These errors were reduced to 2.77 and 4.16 cm/s in the assimilation case.The overall improvement at these comparison sites was 26 and 45% in the cross-and along-bank directions.The assimilated subtidal velocity exhibited the same spatial distribution patterns as observations at seven ADCP sites.In Figure 10, for example, the data showed a subsurface along-bank velocity maximum at a depth of 20-30 m on 11 June 1999.This feature was reproduced in the model solution with assimilation.A relatively strong vertical shear of the cross-bank subtidal velocity was found in the upper 20 m on the same day, which also appeared in the vertical distribution of the assimilated subtidal flow.u  and v  are the standard deviation of the along-and cross-bank current differences, respectively.
The assimilated subtidal velocity exhibited the same spatial distribution patterns as observations at seven ADCP sites.In Figure 10, for example, the data showed a subsurface along-bank velocity maximum at a depth of 20-30 m on 11 June 1999.This feature was reproduced in the model solution with assimilation.A relatively strong vertical shear of the cross-bank subtidal velocity was found in the upper 20 m on the same day, which also appeared in the vertical distribution of the assimilated subtidal flow.

Discussion
Our simple model experiments show that the adjoint assimilation method can minimize model data differences in a least-square sense for current assimilation.However, two critical issues arise regarding this method: (1) computational efficiency and (2) convergence.We had difficulty using the ε-criterion to control the number of iterations required for a longterm (three-month) assimilation since the time required to reach a specific criterion was too long.Therefore, based on the time required to run an ocean forecast model, the adjoint data assimilation approach used in these experiments seems inefficient and practical enough to be used for the forecast application in the coastal ocean without a high-performance MPI cluster computer.
Convergence is sensitive to the assimilation window and iterations.Unlike the open ocean, many coastal regions are characterized by dominant semidiurnal tidal motion.Previous current measurements on GB revealed that the independent time scale of the tidal flow is about one hour [35].This means that a semidiurnal tidal cycle only contains about 12 independent data records.Using a short time integration window in an adjoint model could reach a faster convergence rate toward the observation, but it is accompanied by a sacrifice of statistical confidence level.Determinations of a proper integration window with a rapid convergence rate and sufficient statistical confidence are the critical issues that should be considered when the adjoint data assimilation method is applied to the tidal-dominant coastal oceans.
Our numerical experiments showed that the error convergence rate in the assimilated velocity was sensitive to the forward/backward integration windows.When the window length was close to the M 2 tidal period, the adjoint model failed to reach a convergence solution.When a 12-h window was chosen, the standard deviation error between assimilated and observed along-bank velocities quickly decreased in the first three iterations, but gradually increased after the 6th iteration (Figure 11).As the window length was reduced to 6 h (about half a M 2 tidal period), the convergence was significantly improved.The standard deviation error exhibited a rapid decrease rate in the first ten iterations and gradually converged to the observation as iterations increased (Figure 11).The deviation error was reduced to 2.0 cm/s after 40 iterations.
Convergence is sensitive to the assimilation window and iterations.Unlike the open ocean, many coastal regions are characterized by dominant semidiurnal tidal motion.Previous current measurements on GB revealed that the independent time scale of the tidal flow is about one hour [35].This means that a semidiurnal tidal cycle only contains about 12 independent data records.Using a short time integration window in an adjoint model could reach a faster convergence rate toward the observation, but it is accompanied by a sacrifice of statistical confidence level.Determinations of a proper integration window with a rapid convergence rate and sufficient statistical confidence are the critical issues that should be considered when the adjoint data assimilation method is applied to the tidal-dominant coastal oceans.
Our numerical experiments showed that the error convergence rate in the assimilated velocity was sensitive to the forward/backward integration windows.When the window length was close to the M2 tidal period, the adjoint model failed to reach a convergence solution.When a 12-h window was chosen, the standard deviation error between assimilated and observed along-bank velocities quickly decreased in the first three iterations, but gradually increased after the 6th iteration (Figure 11).As the window length was reduced to 6 h (about half a M2 tidal period), the convergence was significantly improved.The standard deviation error exhibited a rapid decrease rate in the first ten iterations and gradually converged to the observation as iterations increased (Figure 11).The deviation error was reduced to 2.0 cm/s after 40 iterations.In our experiments, no attempts have been made to assimilate the temperature and salinity because sufficient hydrographic data were unavailable across the entire study area during the ADCP measurement period.It is still being determined whether the convergence rate can change as more control parameters are added to this M 2 tidal-dominant coastal region.

Summary
The ECOM-si adjoint model was developed and applied to assimilate the March-June 1999 flow field across the southern flank of GB.The current measurement data in the assimilation were from seven ADCP tripod sites centered at the 60 m isobath.Numerical experiments were carried out to understand critical issues affecting the assimilation efficiency and accuracy.Without assimilation, the model was robust enough to reproduce the temporospatial variability of tidal and subtidal flows observed at seven ADCP measurement sites.However, the tidal currents were underestimated near the bottom and overestimated near the surface, with an overall error of 2-4 cm/s.Adding wind stress and surface heat flux in the model captured the variability of the low-frequency flow, but the subtidal flow was significantly overpredicted during relatively strong wind periods.The maximum simulation error was up to ~20 cm/s.The adjoint assimilation significantly improved tidal and subtidal current predictions, especially during the storm periods in April and May.The assimilation accounted for a 45% improvement for the along-bank current and 26% for the cross-bank current.
The convergence rate of the assimilated currents toward observed currents depends on the assimilation window and iteration number.When the window length was close to the M 2 tidal period, the assimilated flow could not converge to the observed flow as the iteration number increased.Shortening the window length provided a faster convergence rate, but the statistical confidence level was reduced.As a practical compromise, a 6-+h window was recommended for the semidiurnal tidal dominant region.In our experiments, this window length led to a rapid convergence rate in the first ten iterations.After 40 iterations, the model reached a numerical solution that satisfied a given ε-criterion.
Our assimilation experiments focused only on using observed currents.Since no hydrographic measurements were made near the ADCP array during the study period, no efforts were made to assimilate temperature and salinity data.

Figure 1 .
Figure 1.The ADCP mooring locations and numerical computational domain (a) and numerical grids used in the 2-D model (b).In (a), filled triangles and filled dots are the locations of the NOAA environmental buoy 44011 and ADCP measurement sites.In (b), the upper panel grids were plotted every 4 points in the vertical and 10 points in the horizontal.

Figure 1 .
Figure 1.The ADCP mooring locations and numerical computational domain (a) and numerical grids used in the 2-D model (b).In (a), filled triangles and filled dots are the locations of the NOAA environmental buoy 44011 and ADCP measurement sites.In (b), the upper panel grids were plotted every 4 points in the vertical and 10 points in the horizontal.
to judge if the model provided an optimal flow field of u and v.

Figure 2 .
Figure 2. Schematic chart of the numerical operations for adjoint assimilation.

Figure 3 .
Figure 3.The temperature and salinity initial conditions used in the 2-D experiments.The data used to construct these conditions were from the March-April climatologically averaged dataset [33].

Figure 4 .
Figure 4.The time series of surface wind stress calculated using the meteorological data from 20 March to 17 June 1999, on the NOAA environmental buoy 44011 (shown in Figure 1).

Figure 3 . 18 Figure 3 .
Figure 3.The temperature and salinity initial conditions used in the 2-D experiments.The data used to construct these conditions were from the March-April climatologically averaged dataset [33].

Figure 4 .
Figure 4.The time series of surface wind stress calculated using the meteorological data from 20 March to 17 June 1999, on the NOAA environmental buoy 44011 (shown in Figure 1).

Figure 4 .
Figure 4.The time series of surface wind stress calculated using the meteorological data from 20 March to 17 June 1999, on the NOAA environmental buoy 44011 (shown in Figure 1).

Figure 5 .
Figure 5. Comparisons between simulated (dashed line) and observed tidal (dots) currents (major and minor axes) at the ADCP site 567.The horizontal bar line: the measurement uncertainty.

Figure 5 .
Figure 5. Comparisons between simulated (dashed line) and observed tidal (dots) currents (major and minor axes) at the ADCP site 567.The horizontal bar line: the measurement uncertainty.

Figure 6 .
Figure 6.Comparison between simulated (dashed line) and observed (solid line) currents at selected depths of 5, 20, and 50 m at the ADCP site 567 during April 1999.Here, u is the cross-bank velocity and v is the along-bank velocity.

Figure 6 .
Figure 6.Comparison between simulated (dashed line) and observed (solid line) currents at selected depths of 5, 20, and 50 m at the ADCP site 567 during April 1999.Here, u is the cross-bank velocity and v is the along-bank velocity.

Figure 7 .
Figure 7.Comparison between simulated (dashed line) and observed (solid line) subtidal currents at selected depths of 5, 20, and 50 m at the ADCP site 567 over 20 March-15 June 1999.In the figure, u is the cross-bank velocity and v is the along-bank velocity.

Figure 7 .
Figure 7.Comparison between simulated (dashed line) and observed (solid line) subtidal currents at selected depths of 5, 20, and 50 m at the ADCP site 567 over 20 March-15 June 1999.In the figure, u is the cross-bank velocity and v is the along-bank velocity.

Figure 8 .
Figure 8.Comparison between assimilated (dashed line) and observed (solid line) currents at depths of 5, 20, and 50 m at the ADCP site 567 during April 1999.u and v are the cross-and along-bank components of the velocity.

Figure 8 .
Figure 8.Comparison between assimilated (dashed line) and observed (solid line) currents at depths of 5, 20, and 50 m at the ADCP site 567 during April 1999.u and v are the cross-and along-bank components of the velocity.

Figure 8 .
Figure 8.Comparison between assimilated (dashed line) and observed (solid line) currents at depths of 5, 20, and 50 m at the ADCP site 567 during April 1999.u and v are the cross-and along-bank components of the velocity.

Figure 9 .
Figure 9.Comparison between assimilated (dashed line) and observed (solid line) subtidal currents at depths of 5, 20, and 50 m at the ADCP site 567 over 20 March-15 June 1999.u and v are the cross-and along-bank components of the velocity.

Figure 11 .Figure 11 .
Figure 11.Upper panel: comparisons between assimilated and observed along-bank velocities over an M2 tidal cycle on 24-25 March 1999 for the cases with iteration numbers of 1, 10, and 39.Lower

Table 1 .
Difference between model-predicted and observed tidal current ellipse parameters.

Table 2 .
The standard deviation of the difference between simulated/assimilated and observed along-bank and cross-bank currents.