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.

Identification of thermal front dynamics in the northern Malacca Strait using ROMS 3D-model

The thermal front in the oceanic system is believed to have a significant effect on biological activity.During an era of climate change,changes in heat regulation between the atmosphere and oceanic interior can alter the characteristics of this important feature.Using the simulation results of the 3D Regional Ocean Modelling System(ROMS),we identified the location of thermal fronts and determined their dynamic variability in the area between the southern Andaman Sea and northern Malacca Strait.The Single Image Edge Detection(SIED)algorithm was used to detect the thermal front from model-derived temperature.Results show that a thermal front occurred every year from 2002 to 2012 with the temperature gradient at the location of the front was 0.3°C/km.Compared to the years affected by El Ni?o and negative Indian Ocean Dipole(IOD),the normal years(e.g.,May 2003)show the presence of the thermal front at every selected depth(10,25,50,and 75 m),whereas El Ni?o and negative IOD during 2010 show the presence of the thermal front only at depth of 75 m due to greater warming,leading to the thermocline deepening and enhanced stratification.During May 2003,the thermal front was separated by cooler SST in the southern Andaman Sea and warmer SST in the northern Malacca Strait.The higher SST in the northern Malacca Strait was believed due to the besieged Malacca Strait,which trapped the heat and make it difficult to release while higher chlorophyll a in Malacca Strait is due to the freshwater conduit from nearby rivers(Klang,Langat,Perak,and Selangor).Furthermore,compared to the southern Andaman Sea,the chlorophyll a in the northern Malacca Strait is easier to reach the surface area due to the shallower thermocline,which allows nutrients in the area to reach the surface faster.

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.

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.

Climate variability and predictability associated with the Indo-Pacific Oceanic Channel Dynamics in the CCSM4 Coupled System Model

An experiment using the Community Climate System Model(CCSM4), a participant of the Coupled Model Intercomparison Project phase-5(CMIP5), is analyzed to assess the skills of this model in simulating and predicting the climate variabilities associated with the oceanic channel dynamics across the Indo-Pacific Oceans. The results of these analyses suggest that the model is able to reproduce the observed lag correlation between the oceanic anomalies in the southeastern tropical Indian Ocean and those in the cold tongue in the eastern equatorial Pacific Ocean at a time lag of 1 year. This success may be largely attributed to the successful simulation of the interannual variations of the Indonesian Throughflow, which carries the anomalies of the Indian Ocean Dipole(IOD) into the western equatorial Pacific Ocean to produce subsurface temperature anomalies, which in turn propagate to the eastern equatorial Pacific to generate ENSO. This connection is termed the "oceanic channel dynamics" and is shown to be consistent with the observational analyses. However, the model simulates a weaker connection between the IOD and the interannual variability of the Indonesian Throughflow transport than found in the observations. In addition, the model overestimates the westerly wind anomalies in the western-central equatorial Pacific in the year following the IOD, which forces unrealistic upwelling Rossby waves in the western equatorial Pacific and downwelling Kelvin waves in the east. This assessment suggests that the CCSM4 coupled climate system has underestimated the oceanic channel dynamics and overestimated the atmospheric bridge processes.

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.

Spatiotemporal characteristics of water exchange between the Andaman Sea and the Bay of Bengal

A high-resolution customized numerical model is used to analyze the water transport in the three major water passages between the Andaman Sea(AS)and the Bay of Bengal,i.e.,the Preparis Channel(PC),the Ten Degree Channel(TDC),and the Great Channel(GC),based on the daily averaged simulation results ranging from 2010 to 2019.Spectral analysis and Empirical Orthogonal Function(EOF)methods are employed to investigate the spatiotemporal variability of the water exchange and controlling mechanisms.The results of model simulation indicate that the net average transports of the PC and GC,as well as their linear trend,are opposite to that of the TDC.This indicates that the PC and the GC are the main inflow channels of the AS,while the TDC is the main outflow channel of the AS.The transport variability is most pronounced at surface levels and between 100 m and 200 m depth,likely affected by monsoons and circulation.A 182.4-d semiannual variability is consistently seen in all three channels,which is also evident in their second principal components.Based on sea level anomalies and EOF analysis results,this is primarily due to equatorial winds during the monsoon transition period,causing eastward movement of Kelvin waves along the AS coast,thereby affecting the spatiotemporal characteristics of the flow in the AS.The first EOF of the PC flow field section shows a split at 100 m deep,likely due to topography.The first EOF of the TDC flow field section is steady but has potent seasonal oscillations in its time series.Meanwhile,the first EOF of the GC flow field section indicates a stable surface inflow,probably influenced by the equatorial Indian Ocean’s eastward current.

Estimation of free core nutation parameters and availability of computing options

The Earth’s Free Core Nutation(FCN) causes Earth tides and forced nutation with frequencies close to the FCN that exhibit resonance effects.High-precision superconducting gravimeter(SG) and very long baseline interferometry(VLBI) provide good observation techniques for detecting the FCN parameters.However,some choices in data processing and solution procedures increase the uncertainty of the FCN parameters.In this study,we analyzed the differences and the effectiveness of weight function and ocean tide corrections in the FCN parameter detection using synthetic data,SG data from thirty-one stations,and the 10 celestial pole offset(CPO) series.The results show that significant discrepancies are caused by different computing options for a single SG station.The stacking method,which results in a variation of0.24-5 sidereal days(SDs) in the FCN period(T) and 10^(3)-10^(4) in the quality factor(Q) due to the selection of the weighting function and the ocean tide model(OTM),can effectively suppress this influence.The statistical analysis results of synthetic data shows that although different weight choices,while adjusting the proportion of diurnal tidal waves involved,do not significantly improve the accuracy of fitted FCN parameters from gravity observations.The study evaluated a series of OTMs using the loading correction efficiency.The fitting of FCN parameters can be improved by selecting the mean of appropriate OTMs based on the evaluation results.Through the estimation of the FCN parameters based on the forced nutation,it was found that the weight function P_(1) is more suitable than others,and different CPO series(after 2009) resulted in a difference of 0.4 SDs in the T and of 103 in the Q.We estimated the FCN parameters for SG(T=430.4±1.5 SDs and Q=1.52×10^(4)±2.5×10^(3)) and for VLBI(T=429.8±0.7 SDs,Q=1.88×10^(4)±2.1×10^(3)).

Numerical simulations of Atlantic meridional overturning circulation(AMOC)from OMIP experiments and its sensitivity to surface forcing

Atlantic meridional overturning circulation(AMOC)plays an important role in transporting heat meridionally in the Earth’s climate system and is also a key metrical tool to verify oceanic general circulation models.Two OMIP(Ocean Model Intercomparison Project phase 1 and 2)simulations with LICOM3(version 3 of the LASG/IAP Climate System Ocean Model)developed at the State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics(LASG),Institute of Atmospheric Physics(IAP),are compared in this study.Both simulations well reproduce the fundamental characteristics of the AMOC,but the OMIP1 simulation shows a significantly stronger AMOC than the OMIP2 simulation.Because the LICOM3 configurations are identical between these two experiments,any differences between them must be attributed to the surface forcing data.Further analysis suggests that sea surface salinity(SSS)differences should be mainly responsible for the enhanced AMOC in the OMIP1 simulation,but sea surface temperature(SST)also play an unignorable role in modulating AMOC.In the North Atlantic,where deep convection occurs,the SSS in OMIP1 is more saline than that in OMIP1.We find that in the major region of deep convection,the change of SSS has more significant effect on density than the change of SST.As a result,the SSS was more saline than that in OMIP2,leading to stronger deep convection and subsequently intensify the AMOC.We conduct a series of numerical experiments with LICOM3,and the results confirmed that the changes in SSS have more significant effect on the strength of AMOC than the changes in SST.

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.

LICOM Model Datasets for the CMIP6 Ocean Model Intercomparison Project

The datasets of two Ocean Model Intercomparison Project(OMIP)simulation experiments from the LASG/IAP Climate Ocean Model,version 3(LICOM3),forced by two different sets of atmospheric surface data,are described in this paper.The experiment forced by CORE-II(Co-ordinated Ocean–Ice Reference Experiments,Phase II)data(1948–2009)is called OMIP1,and that forced by JRA55-do(surface dataset for driving ocean–sea-ice models based on Japanese 55-year atmospheric reanalysis)data(1958–2018)is called OMIP2.First,the improvement of LICOM from CMIP5 to CMIP6 and the configurations of the two experiments are described.Second,the basic performances of the two experiments are validated using the climatological-mean and interannual time scales from observation.We find that the mean states,interannual variabilities,and long-term linear trends can be reproduced well by the two experiments.The differences between the two datasets are also discussed.Finally,the usage of these data is described.These datasets are helpful toward understanding the origin system bias of the fully coupled model.

A Numerical World Ocean General Circulation Model

This paper describes a numerical model of the world ocean based on the fully primitive equations. A 'Standard' ocean state is introduced into the equations of the model and the perturbed thermodynamic variables are used in the modlc's calculations. Both a free upper surface and a bottom topography are included in the model and a sigma coordinate is used to normalize the model's vertical component. The model has four unevenly-spaced layers and 4 × 5 horizontal resolution based on C-grid system. The finite-difference scheme of the model is designed to conserve the gross available energy in order to avoid fictitious energy generation or decay.The model has been tested in response to the annual mean surface wind stress, sea level air pressure and sea level air temperature as a preliminary step to its further improvement and its coupling with a global atmospheric general circulation model. Some of results, including currents, temperature and sea surface elevation simulated by the mode! arc presented.

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.

An extended variable-grid global ocean circulation model and its preliminary results of the equatorial Pacific circulation

To investigate the interaction between the tropical Pacific and China seas a variable-grid global ocean circulation model with fine grid covering the area from 20°S to 50°N and from 99° to 150°E is developed. Numerical computation of the annually cyclic circulation fields is performed. The results of the annual mean zonal currents and deep to abyssal western boundary currents in the equatorial Pacific Ocean are reported. The North Equatorial Current,the North Equatorial Countercurrent, the South Equatorial Current and the Equatorial Undercurrent are fairly well simulated. The model well reproduces the northward flowing abyssal western boundary current.From the model results a lower deep western boundary current east of the Bismarck-Solomon-New Hebrides Island chain at depths around 2 000 m has been found. The model results also show that the currents in the equatorial Pacific Ocean have multi-layer structures both in zonal currents and western boundary currents, indicating that the global ocean overturning thermohaline circulation appears of multi-layer pattern.

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.

Asymptotic solution of a weak nonlinear model for the mid-latitude stationary wind field of a two-layer barotropic ocean

A weak nonlinear model of a two-layer barotropic ocean with Rayleigh dissipation is built.The analytic asymptotic solution is derived in the mid-latitude stationary wind field,and the physical meaning of the corresponding problem is discussed.