Assimilating the along-track sea level anomaly into the regional ocean modeling system using the ensemble optimal interpolation

The ensemble optimal interpolation （EnOI） is applied to the regional ocean modeling system （ROMS） with the ability to assimilate the along-track sea level anomaly （TSLA）. This system is tested with an eddy-resolving system of the South China Sea （SCS）. Background errors are derived from a running seasonal ensemble to account for the seasonal variability within the SCS. A fifth-order localization function with a 250 km localization radius is chosen to reduce the negative effects of sampling errors. The data assimilation system is tested from January 2004 to December 2006. The results show that the root mean square deviation （RMSD） of the sea level anomaly decreased from 10.57 to 6.70 cm, which represents a 36.6% reduction of error. The data assimilation reduces error for temperature within the upper 800 m and for salinity within the upper 200 m, although error degrades slightly at deeper depths. Surface currents are in better agreement with trajectories of surface drifters after data assimilation. The variance of sea level improves significantly in terms of both the amplitude and position of the strong and weak variance regions after assimilating TSLA. Results with AGE error （AGE） perform better than no AGE error （NoAGE） when considering the improvements of the temperature and the salinity. Furthermore, reasons for the extremely strong variability in the northern SCS in high resolution models are investigated. The results demonstrate that the strong variability of sea level in the high resolution model is caused by an extremely strong Kuroshio intrusion. Therefore, it is demonstrated that it is necessary to assimilate the TSLA in order to better simulate the SCS with high resolution models.

Down-scaled regional ocean modeling system (ROMS) for high-resolution coastal hydrodynamics in Korea

A down-scaled operational oceanographic system is developed for the coastal waters of Korea using a re- gional ocean modeling system （ROMS）. The operational oceanographic modeling system consists of at- mospheric and hydrodynamic models. The hydrodynamic model, ROMS, is coupled with wave, sediment transport, and water quality modules. The system forecasts the predicted results twice a day on a 72 h basis, including sea surface elevation, currents, temperature, salinity, storm surge height, and wave information for the coastal waters of Korea. The predicted results are exported to the web-GIS-based coastal informa- tion system for real-time dissemination to the public and validation with real-time monitoring data using visualization technologies. The ROMS is two-way coupled with a simulating waves nearshore model, SWAN, for the hydrodynamics and waves, nested with the meteorological model, WRE for the atmospheric surface forcing, and externally nested with the eutrophication model, CE-QUAL-ICM, for the water quality. The op- erational model, ROMS, was calibrated with the tidal surface observed with a tide-gage and verified with current data observed by bottom-mounted ADCP or AWAC near the coastal waters of Korea. To validate the predicted results, we used real-time monitoring data derived from remote buoy system, HF-radar, and geostationary ocean color imager （GOCI）. This down-scaled operational coastal forecasting system will be used as a part of the Korea operational oceanographic system （KOOS） with other operational oceanographic systems.

Overview on high-resolution ocean modeling in JAMSTEC

In view of the importance of ocean component for representing climate change,efforts are underway to implement a high-resolution nesting model system in Model for Interdisciplinary Research on Climate(MIROC) for the North Pacific using the same ocean model as used in the coupled model MIROC5. By comparing double(10 km for the northwestern Pacific,50 km for the rest of the Pacific) and triple(double nesting plus 2 km resolution near Japan) nesting,it turns out that relative vorticity is drastically enhanced near Japan with 2 km resolution. It is hoped that such an elaborated nesting system will reveal detailed processes for the ocean heat uptake by,e.g.,intermediate water and mode water formation for which the"perturbed region"near Japan is the key region.

The performance of a z-level ocean model in modeling the global tide

The performance of a z-level ocean model, the Modular Ocean Model Version 4（MOM4）, is evaluated in terms of simulating the global tide with different horizontal resolutions commonly used by climate models. The performance using various sets of model topography is evaluated. The results show that the optimum filter radius can improve the simulated co-tidal phase and that better topography quality can lead to smaller rootmean square（RMS） error in simulated tides. Sensitivity experiments are conducted to test the impact of spatial resolutions. It is shown that the model results are sensitive to horizontal resolutions. The calculated absolute mean errors of the co-tidal phase show that simulations with horizontal resolutions of 0.5° and 0.25° have about 35.5% higher performance compared that with 1° model resolution. An internal tide drag parameterization is adopted to reduce large system errors in the tidal amplitude. The RMS error of the best tuned 0.25° model compared with the satellite-altimetry-constrained model TPXO7.2 is 8.5 cm for M_2. The tidal energy fluxes of M_2 and K_1 are calculated and their patterns are in good agreement with those from the TPXO7.2. The correlation coefficients of the tidal energy fluxes can be used as an important index to evaluate a model skill.

Fourth-Order Conservative Transport on Overset Grids Using Multi-Moment Constrained Finite Volume Scheme for Oceanic Modeling

With an increase in model resolution,compact high-order numerical advection scheme can improve its effectiveness and competitiveness in oceanic modeling due to its high accuracy and scalability on massive-processor computers.To provide high-quality numerical ocean simulation on overset grids,we tried a novel formulation of the fourth-order multi-moment constrained finite volume scheme to simulate continuous and discontinuous problems in the Cartesian coordinate.Utilizing some degrees of freedom over each cell and derivatives at the cell center,we obtained a two-dimensional(2D)cubic polynomial from which point values on the extended overlap can achieve fourth-order accuracy.However,this interpolation causes a lack of conservation because the flux between the regions are no longer equal;thus,a flux correction is implemented to ensure conservation.A couple of numerical experiments are presented to evaluate the numerical scheme,which confirms its approximately fourth-order accuracy in conservative transportation on overset grid.The test cases reveal that the scheme is effective to suppress numerical oscillation in discontinuous problems,which may be powerful for salinity advection computing with a sharp gradient.

Effect of Spatial and Temporal Scales on Habitat Suitability Modeling:A Case Study of Ommastrephes bartramii in the Northwest Pacific Ocean

Temporal and spatial scales play important roles in fishery ecology,and an inappropriate spatio-temporal scale may result in large errors in modeling fish distribution.The objective of this study is to evaluate the roles of spatio-temporal scales in habitat suitability modeling,with the western stock of winter-spring cohort of neon flying squid （Ornmastrephes bartramii） in the northwest Pacific Ocean as an example.In this study,the fishery-dependent data from the Chinese Mainland Squid Jigging Technical Group and sea surface temperature （SST） from remote sensing during August to October of 2003-2008 were used.We evaluated the differences in a habitat suitability index model resulting from aggregating data with 36 different spatial scales with a combination of three latitude scales （0.5°,1 ° and 2°）,four longitude scales （0.5°,1°,2° and 4°）,and three temporal scales （week,fortnight,and month）.The coefficients of variation （CV） of the weekly,biweekly and monthly suitability index （SI） were compared to determine which temporal and spatial scales of SI model are more precise.This study shows that the optimal temporal and spatial scales with the lowest CV are month,and 0.5° latitude and 0.5° longitude for O.bartramii in the northwest Pacific Ocean.This suitability index model developed with an optimal scale can be cost-effective in improving forecasting fishing ground and requires no excessive sampling efforts.We suggest that the uncertainty associated with spatial and temporal scales used in data aggregations needs to be considered in habitat suitability modeling.

Quantifying the Role of the Eddy Transfer Coefficient in Simulating the Response of the Southern Ocean Meridional Overturning Circulation to Enhanced Westerlies in a Coarse-resolution Model

This study assesses the capability of a coarse-resolution ocean model to replicate the response of the Southern Ocean Meridional Overturning Circulation(MOC) to intensified westerlies,focusing on the role of the eddy transfer coefficient(κ).κ is a parameter commonly used to represent the velocities induced by unresolved eddies.Our findings reveal that a stratification-dependent κ,incorporating spatiotemporal variability,leads to the most robust eddy-induced MOC response,capturing 82% of the reference eddy-resolving simulation.Decomposing the eddy-induced velocity into its vertical variation(VV) and spatial structure(SS) components unveils that the enhanced eddy compensation response primarily stems from an augmented SS term,while the introduced VV term weakens the response.Furthermore,the temporal variability of the stratification-dependent κ emerges as a key factor in enhancing the eddy compensation response to intensified westerlies.The experiment with stratification-dependent κ exhibits a more potent eddy compensation response compared to the constant κ,attributed to the structure of κ and the vertical variation of the density slope.These results underscore the critical role of accurately representing κ in capturing the response of the Southern Ocean MOC and emphasize the significance of the isopycnal slope in modulating the eddy compensation mechanism.

Tidal modeling based on satellite altimetry observations of TOPEX/ Poseidon, Jason1, Jason2, and Jason3 with high prediction capability: A case study of the Baltic Sea

This research aims to optimize the utilization of long-term sea level data from the TOPEX/Poseidon,Jason1,Jason2,and Jason3 altimetry missions for tidal modeling.We generate a time series of along-track observations and apply a developed method to produce tidal models with specific tidal constituents for each location.Our tidal modeling methodology follows an iterative process:partitioning sea surface height(SSH)observations into analysis/training and prediction/validation parts and ultimately identi-fying the set of tidal constituents that provide the best predictions at each time series location.The study focuses on developing 1256 time series along the altimetry tracks over the Baltic Sea,each with its own set of tidal constituents.Verification of the developed tidal models against the sSH observations within the prediction/validation part reveals mean absolute error(MAE)values ranging from 0.0334 m to 0.1349 m,with an average MAE of 0.089 m.The same validation process is conducted on the FES2014 and EOT20 global tidal models,demonstrating that our tidal model,referred to as BT23(short for Baltic Tide 2023),outperforms both models with an average MAE improvement of 0.0417 m and 0.0346 m,respectively.In addition to providing details on the development of the time series and the tidal modeling procedure,we offer the 1256 along-track time series and their associated tidal models as supplementary materials.We encourage the satellite altimetry community to utilize these resources for further research and applications.

Evaluation of Nonbreaking Wave-Induced Mixing Parameterization Schemes Based on a One-Dimensional Ocean Model

Surface waves have a considerable effect on vertical mixing in the upper ocean.In the past two decades,the vertical mixing induced through nonbreaking surface waves has been used in ocean and climate models to improve the simulation of the upper ocean.Thus far,several nonbreaking wave-induced mixing parameterization schemes have been proposed;however,no quantitative comparison has been performed among them.In this paper,a one-dimensional ocean model was used to compare the performances of five schemes,including those of Qiao et al.(Q),Hu and Wang(HW),Huang and Qiao(HQ),Pleskachevsky et al.(P),and Ghantous and Babanin(GB).Similar to previous studies,all of these schemes can decrease the simulated sea surface temperature(SST),increase the subsurface temperature,and deepen the mixed layer,thereby alleviating the common thermal deviation problem of the ocean model for upper ocean simulation.Among these schemes,the HQ scheme exhibited the weakest wave-induced mixing effect,and the HW scheme exhibited the strongest effect;the other three schemes exhibited roughly the same effect.In particular,the Q and P schemes exhibited nearly the same effect.In the simulation based on observations from the Ocean Weather Station Papa,the HQ scheme exhibited the best performance,followed by the Q scheme.In the experiment with the HQ scheme,the root-mean-square deviation of the simulated SST from the observations was 0.43℃,and the mixed layer depth(MLD)was 2.0 m.As a contrast,the deviations of the SST and MLD reached 1.25℃ and 8.4 m,respectively,in the experiment without wave-induced mixing.

Tide simulation in a global eddy-resolving ocean model

The tide plays a pivotal role in the ocean,affecting the global ocean circulation and supplying the bulk of the energy for the global meridional overturning circulation.To further investigate internal tides and their impacts on circulation,it is imperative to incorporate tidal forcing into the eddy-resolving global ocean circulation model.In this study,we successfully incorporated explicit tides(eight major constituents)into a global eddy-resolving general ocean circulation model and evaluated its tidal simulation ability.We obtained harmonic constants by analyzing sea surface height through tidal harmonic analysis and compared them with the analysis data Topex Poseidon Cross-Overs v9(TPXO9),the open ocean tide dataset from 102 open-ocean tide observations,and tide gauge stations from World Ocean Circulation Experiment.The results demonstrated that the State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics/Institute of Atmospheric Physics(LASG/IAP)Climate System Ocean Model 3.0(LICOM3.0)effectively simulated tides,with errors predominantly occurring in nearshore regions.The tidal amplitude simulated in LICOM3.0 was greater than that of TPXO9,and these high-amplitude areas exhibited greater errors.The amplitude error of the M_(2) constituent was larger,while the phase error of the K_(1) constituent was more significant.Furthermore,we further compared our results with those from other models.

Sea Surface Temperature Simulation of Tropical and North Pacific Basins Using a Hybrid Coordinate Ocean Model (HYCOM)

A Hybrid Coordinate Ocean Model （HYCOM） is used to simulate the sea surface temperature of the Tropical and North Pacific. Based on the different combinations of two air-Sea flux data sets （COADS and ECMWF） and two bulk parameter formulas （non-constant and constant）, four numerical experiments are carried out. The following conclusions can be deduced from the numerical results. （1） The numerical results using non-constant bulk parameter formula are much better than those using constant one. In the Pacific area from 40°N to 20°S, the annual average SST obtained from the experiment using non-constant bulk parameter formula is 0.21 ℃ higher than that from the satellite-based SST climatology （the pathfinder data）. However, the difference is 0.63 ℃ for the experiment when the using constant one. （2） HYCOM successfully simulates the monthly variation of climatological SST in tropical and north Pacific basins and monthly spatial variation of Western Pacific Warm Pool. Especially in the Pacific area from 40°N to 20°S, the difference of the seasonal averaged SST between pathfinder data and the result of experiment 2 （using COADS data set and non-constant bulk parameter formula） is only about 0.02 ℃. （3）The simulation results using different Air-Sea flux data are different and the difference is very large in some regions. In the northwest of the model region, the annual average SST obtained from experiment 2 （using COADS data set） is 1℃ higher than that obtained from experiment 4 （using ECMWF data set）. Contrarily, the result of experiment 4 is 1 ℃ larger than that of experiment 2 in the southeast of the model region. The largest difference is about 4 ℃ occurred near the area of 58°N, 140°E and the Bohai sea.

A Coastal Ocean Model of Semi-Implicit Finite Volume Unstructured Grid

A two-dimensional coastal ocean model based on unstructured C-grid is built, in which the momentum equation is discretized on the faces of each cell, and the continuity equation is discretized on the cell. The model is discretized by semi-implicit finite volume method, in that the free surface is semi-implicit and the bottom friction is implicit, thereby removing stability limitations associated with the surface gravity wave and friction. The remaining terms in the momentum equations are discretized explicitly by integral finite volume method and second-order Adams-Bashforth method. Tidal flow in the polar quadrant with known analytic solution is employed to test the proposed model. Finally, the performance of the present model to simulate tidal flow in a geometrically complex domain is examined by simulation of tidal currents in the Pearl River Estuary.

Remote Efects of Tropical Cyclone Wind Forcing over the Western Pacific on the Eastern Equatorial Ocean

An ocean general circulation model （OGCM） is used to demonstrate remote effects of tropical cyclone wind （TCW） forcing in the tropical Pacific. The signature of TCW forcing is explicitly extracted using a locally weighted quadratic least=squares regression （called as LOESS） method from six-hour satellite surface wind data; the extracted TCW component can then be additionally taken into account or not in ocean modeling, allowing isolation of its effects on the ocean in a clean and clear way. In this paper, seasonally varying TCW fields in year 2008 are extracted from satellite data which are prescribed as a repeated annual cycle over the western Pacific regions off the equator （poleward of 10°N/S）; two long-term OGCM experiments are performed and compared, one with the TCW forcing part included additionally and the other not. Large, persistent thermal perturbations （cooling in the mixed layer （ML） and warming in the thermocline） are induced locally in the western tropical Pacific, which are seen to spread with the mean ocean circulation pathways around the tropical basin. In particular, a remote ocean response emerges in the eastern equatorial Pacific to the prescribed off-equatorial TCW forcing, characterized by a cooling in the mixed layer and a warming in the thermocline. Heat budget analyses indicate that the vertical mixing is a dominant process responsible for the SST cooling in the eastern equatorial Pacific. Further studies are clearly needed to demonstrate the significance of these results in a coupled ocean-atmosphere modeling context.

Modeling the circulation and sediment transport in the Beibu Gulf

Water circulation and sediment transport in the Beibu Gulf are important for its environmental protection and resource exploitation.By employing the Regional Ocean Modeling System（ROMS）,we studied the seasonal variation of circulation,sediment transport and long-term morphological evolution in the Beibu Gulf.The simulation results show that the circulation induced by tide and wind is cyclonic both in winter and summer in the gulf and that the wind-driven circulation is stronger in winter than that in summer.The sediment concentration is higher in the Qiongzhou Strait,west of the Hainan Island and the coast of Vietnam and the Leizhou Peninsula.The sediment is transported westwards in winter and eastwards in summer in the Qiongzhou Strait.The west entrance of the Qiongzhou Strait is dominated by westward transport all the year round.The sediment discharged by rivers is deposited near the river mouths.The simulated result demonstrates that the sediment transport is mainly controlled by tidal induced bottom resuspension in the Beibu Gulf.Four characteristics are summarized for the distribution patterns of erosion and deposition.（1） The erosion and deposition are insignificant in most area of the gulf.（2） Sediment deposition is more significant in the mouths of Qiongzhou Strait.（3） The erosion is observed in the seabed of Qiongzhou Strait.（4） Erosion and deposition occur alternatively in the west of Hainan Island.

Model Evidence for Interdecadal Pathway Changes in the Subtropics and Tropics of the South Pacific Ocean

Numerical simulations using a version of the GFDL/NOAA Modular Ocean Model （MOM 3） are analyzed to demonstrate interdecadal pathway changes from the subtropics to the tropics in the South Pacific Ocean. After the 1976 -77 climate shift, the subtropical gyre of the South Pacific underwent significant changes, characterized by a slowing down in its circulation and a southward displacement of its center by about 5°- 10° latitude on the western side. The associated circulation altered its flow path in the northwestern part of the subtropical gyre, changing from a direct pathway connecting the subtropics to the tropics before the shift to a more zonal one after. This effectively prevented some subtropical waters from directly entering into the western equatorial Pacific. Since waters transported onto the equator around the subtropical gyre are saline and warm, such changes in the direct pathway and the associated reduction in equatorward exchange from the subtropics to the tropics affected water mass properties downstream in the western equatorial Pacific, causing persisted freshening and cooling of subsurface water as observed after the late 1970s. Previously, changes in gyre strength and advection of temperature anomalies have been invoked as mechanisms for linking the subtropics and tropics on interdecadal time scales. Here we present an additional hypothesis in which geographic shifts in the gyre structure and location （a pathway change） could play a similar role.

On the Importance of High-Resolution in Large-Scale Ocean Models

Eddying global ocean models are now routinely used for ocean prediction,and the value-added of a better representation of the observed ocean variability and western boundary currents at that resolution is currently being evaluated in climate models.This overview article begins with a brief summary of the impact on ocean model biases of resolving eddies in several global ocean-sea ice numerical simulations.Then,a series of North and Equatorial Atlantic configurations are used to show that an increase of the horizontal resolution from eddy-resolving to submesoscale-enabled together with the inclusion of high-resolution bathymetry and tides significantly improve the models’abilities to represent the observed ocean variability and western boundary currents.However,the computational cost of these simulations is extremely large,and for these simulations to become routine,close collaborations with computer scientists are essential to ensure that numerical codes can take full advantage of the latest computing architecture.

Numerical study of the seasonal salinity budget of the upper ocean in the Bay of Bengal in 2014

Impact factors on the salinity budget, especially the eddy salt fluxes and smaller-scale diffusive salt fluxes for the upper 50 m of the Bay of Bengal(BoB) in 2014 are investigated using a box model based on the Regional Ocean Modeling System(ROMS) daily outputs. The model results reproduce that the precipitation and river runoff s are the dominant factors modulating the sharp salinity decrease during the summer monsoon season. The analysis shows that the salinity increase after the summer monsoon is mostly due to the meridional advective and diffusive salt fluxes. The vertical advective salt flux, which is sensitive to the different signals of the wind stress curl, plays an important role in balancing the salinity change induced by the meridional advective salt flux during both the summer and winter monsoon seasons. Distinctive spatial mesoscale structures are presented in the eddy salt flux throughout the year, and their contributions are sizeable(over 30% in the meridional direction and about 10%–30% in the vertical direction). The meridional eddy salt flux is larger in the monsoon seasons than that in the inter-monsoon seasons, and in a positive pattern near the western boundary during the winter monsoon and autumn inter-monsoon. The vertical eddy salt flux makes an important contribution to the salinity budget, especially along the coastal area and around the Andaman and Nicobar Islands. The vertical eddy salt flux becomes large when a tropical cyclone passes the area.

Fundamental Framework and Experiments of the ThirdGeneration of IAP/LASG World OceanGeneral Circulation Model

A new generation of the IAP / LASG world ocean general circulation model is designed and presented based on the previous 20-layer model, with enhanced spatial resolutions and improved parameterizations. The model uses a triangular-truncated spectral horizontal grid system with its zonal wave number of 63 (T63) to match its atmospheric counterpart of a T63 spectral atmosphere general circulation model in a planned coupled ocean-atmosphere system. There are 30 layers in vertical direction, of which 20 layers are located above 1000 m for better depicting the permanent thermocline. As previous ocean models developed in IAP / LASG, a free surface (rather than “rigid-lid” approximation) is included in this model. Compared with the 20-layer model, some more detailed physical parameterizations are considered, including the along / cross isopycnal mixing scheme adapted from the Gent-MacWilliams scheme. The model is spun up from a motionless state. Initial conditions for temperature and salinity are taken from the three-dimensional distributions of Levitus’ annual mean observation. A preliminary analysis of the first 1000-year integration of a control experiment shows some encouraging improvements compared with the twenty-layer model, particularly in the simulations of permanent thermocline, thermohaline circulation, meridional heat transport, etc. resulted mainly from using the isopycnal mixing scheme. However, the use of isopycnal mixing scheme does not significantly improve the simulated equatorial thermocline. A series of numerical experiments show that the most important contribution to the improvement of equatorial thermocline and the associated equatorial under current comes from reducing horizontal viscosity in the equatorial regions. It is found that reducing the horizontal viscosity in the equatorial Atlantic Ocean may slightly weaken the overturning rate of North Atlantic Deep Water.

Progress in the Development and Application of Climate Ocean Models and Ocean-Atmosphere Coupled Models in China

A review is presented about the development and application of climate ocean models and oceanatmosphere coupled models developed in China as well as a review of climate variability and climate change studies performed with these models. While the history of model development is briefly reviewed, emphasis has been put on the achievements made in the last five years. Advances in model development are described along with a summary on scientific issues addressed by using these models. The focus of the review is the climate ocean models and the associated coupled models, including both global and regional models, developed at the Institute of Atmospheric Physics, Chinese Academy of Sciences. The progress of either coupled model development made by other institutions or climate modeling using internationally developed models also is reviewed.

An Eddy-Permitting Oceanic General Circulation Model and Its Preliminary Evaluation

An eddy-permitting, quasi-global oceanic general circulation model, LICOM (LASG/IAP (State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics, Institute of Atmospheric Physics) Climate System Ocean Model), with a uniform grid of 0.5? × 0.5? is established. Forced by wind stresses from Hellerman and Rosenstain (1983), a 40-yr integration is conducted with sea surface temperature and salinity being restored to the Levitus 94 datasets. The evaluation of the annual mean climatology of the LICOM control run shows that the large-scale circulation can be well reproduced. A comparison between the LICOM control run and a parallel integration of L30T63, which has the same framework but a coarse resolution, is also made to con?rm the impact of resolution on the model performance. On account of the reduction of horizontal viscosity with the enhancement of the horizontal resolution, LICOM improves the simulation with respect to not only the intensity of the large scale circulations, but also the magnitude and structure of the Equatorial Undercurrent and South Equatorial Current. Taking advantage of the ?ne grid size, the pathway of the Indonesian Through?ow (ITF) is better represented in LICOM than in L30T63. The transport of ITF in LICOM is more convergent in the upper layer. As a consequence, the Indian Ocean tends to get warmer in LICOM. The poleward heat transports for both the global and individual basins are also signi?cantly improved in LICOM. A decomposed analysis indicates that the transport due to the barotropic gyre, which primarily stands for the barotropic e?ect of the western boundary currents, plays a crucial role in making the di?erence.