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.

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.

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.

A “Dressed” Ensemble Kalman Filter Using the Hybrid Coordinate Ocean Model in the Pacific

The computational cost required by the Ensemble Kalman Filter （EnKF） is much larger than that of some simpler assimilation schemes, such as Optimal Interpolation （OI） or three-dimension variational as- similation （3DVAR）. Ensemble optimal interpolation （EnOI）, a crudely simplified implementation of EnKF, is sometimes used as a substitute in some oceanic applications and requires much less computational time than EnKF. In this paper, to compromise between computational cost and dynamic covariance, we use the idea of ＂dressing＂ a small size dynamical ensemble with a larger number of static ensembles in order to form an approximate dynamic covariance. The term ＂dressing＂ means that a dynamical ensemble seed from model runs is perturbed by adding the anomalies of some static ensembles. This dressing EnKF （DrEnKF for short） scheme is tested in assimilation of real altimetry data in the Pacific using the HYbrid Coordinate Ocean Model （HYCOM） over a four-year period. Ten dynamical ensemble seeds are each dressed by 10 static ensemble members selected from a 100-member static ensemble. Results are compared to two EnKF assimilation runs that use 10 and 100 dynamical ensemble members. Both temperature and salinity fields from the DrEnKF and the EnKF are compared to observations from Argo floats and an OI SST dataset. The results show that the DrEnKF and the 100-member EnKF yield similar root mean square errors （RMSE） at every model level. Error covariance matrices from the DrEnKF and the 100-member EnKF are also compared and show good agreement.

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.

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.

Simulation of Typhoon Muifa using a mesoscale coupled atmosphere–ocean model

A mesoscale coupled atmosphere–ocean model has been developed based on the GRAPES（Global and Regional Assimilation and Prediction System） regional typhoon model（GRAPES_TYM） and ECOM-si（estuary, coast and ocean model（semi-implicit））. Coupling between the typhoon and ocean models was conducted by exchanging wind stress, heat, moisture fluxes, and sea surface temperatures（SSTs） using the coupler OASIS3.0. Numerical prediction experiments were run with and without coupling for the case of Typhoon Muifa in the western North Pacific. To investigate the impact of using more accurate SST information on the simulation of the track and the intensity of Typhoon Muifa, experiments were also conducted using increased SST resolution in the initial condition field of the control test. The results indicate that increasing SST resolution in the initial condition field somewhat improved the intensity forecast, and use of the coupled model improved the intensity forecast significantly, with mean absolute errors in maximum wind speed within 48 and 72 h reduced by 32% and 20%, respectively. Use of the coupled model also resulted in less pronounced over-prediction of the intensity of Typhoon Muifa by the GRAPES_TYM. Moreover, the effects of using the coupled model on the intensity varied throughout the different stages of the development of Muifa owing to changes in the oceanic mixed layer depth. The coupled model had pronounced effects during the later stage of Muifa but had no obvious effects during the earlier stage. The SSTs predicted by the coupled model decreased by about 5–6℃ at most after the typhoon passed, in agreement with satellite data. Furthermore, based on analysis on the sea surface heat flux, wet static energy of the boundary layer, atmospheric temperature, and precipitation forecasted by the coupled model and the control test, the simulation results of this coupled atmosphere–ocean model can be considered to reasonably reflect the primary mechanisms underlying the interactions between tropical cyclones and oceans.

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

评论关于联合模型在中国以及研究与这些模型一起执行了的气候可变性和气候变化的评论开发了的气候海洋模型和海洋空气的发展和应用程序被介绍。当模型开发的历史简短被考察时，强调被放在在最后五年里做的成就上。在模型开发的进展在使用这些模型处理的科学问题上与一篇摘要一起被描述。评论的焦点是气候海洋包括全球、地区性的模型，模型和联系联合模型在大气的物理的研究所发展了，中国科学院。联合的任何一个的进步为另外的机构或气候用国际性地发达的模型当模特儿做的开发建模也被考察。

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.

Estimates of Anthropogenic CO_2 Uptake in a Global Ocean Model

A global ocean general circulation model （L30T63） is employed to study the uptake and distribution of anthropogenic CO2 in the ocean. A subgrid-scale mixing scheme called GM90 is used in the model. There are two main GM90 parameters including isopycnal diffusivity and skew （thickness） diffusivity. Sensitivities of the ocean circulation and the redistribution of dissolved anthropogenic CO2 to these two parameters are examined. Two runs estimate the global oceanic anthropogenic CO2 uptake to be 1.64 and 1.73 Pg C yr^-1 for the 1990s, and that the global ocean contained 86.8 and 92.7 Pg C of anthropogenic CO2 at the end of 1994, respectively. Both the total inventory and uptake from our model are smaller than the data-based estimates. In this presentation, the vertical distributions of anthropogenic CO2 at three meridional sections are discussed and compared with the available data-based estimates. The inventory in the individual basins is also calculated. Use of large isopycnal diffusivity can generally improve the simulated results, including the exchange flux, the vertical distribution patterns, inventory, storage, etc. In terms of comparison of the vertical distributions and column inventory, we find that the total inventory in the Pacific Ocean obtained from our model is in good agreement with the data-based estimate, but a large difference exists in the Atlantic Ocean, particularly in the South Atlantic. The main reasons are weak vertical mixing and that our model generates small exchange fluxes of anthropogenic CO2 in the Southern Ocean. Improvement in the simulation of the vertical transport and sea ice in the Southern Ocean is important in future work.

Meso-scale eddy in the South China Sea simulated by an eddy-resolving ocean model

Mesoscale eddies （MEs） in the South China Sea （SCS） simulated by a quasi-global eddy-resolving ocean general circulation model are evaluated against satellite data during 1993-2007. The modeled ocean data show more activity than shown by the satellite data and reproduces more eddies in the SCS. A total of 345 （428） cyclonic eddies （CEs） and 330 （371） anti-cyclonic eddies （AEs） generated for satellite （model） data are identified during the study period, showing increase of -24% and -12% for the model data, respectively. Compared with eddies in satellite, the simulated eddies tend to have smaller radii, larger amplitudes, a slightly longer lifetime, faster movement and rotation speed, a slightly larger nonlinear properties （U/c） in the model. However, the spatial distribution of generated eddies appears to be inhomogeneous, with more CEs in the northern part of SCS and fewer AEs in the southern part. This is attributed to the exaggerated Kuroshio intrusion in the model because the small islands in the Luzon Strait are still not well resolved although the horizontal resolution reaches （1/10）°. The seasonal variability in the number and the amplitude of eddies generated is also investigated.

ADI BAROTROPIC OCEAN MODEL FOR SIMULATON OF KUROSHIO INTRUSION INTO CHINA SOUTHEASTERN WATERS

In this numerical model for simulating the Kuroshio intrusion into the East and South China Seas,vertically averaged marine hydrodynamic equations governing ocean currents and long-period waves areapproximated by a set of two-time-level semi-implicit finite difference equations. The major terms in-cluding the local acceleration, sea-surface slope, Coriolis force and the bottom friction are approxi-mated with the Crank-Nicholson scheme, which is of second order accuracy. The advection terms are app-roximated with the Leith scheme. The difference equations are split into two sets of alternating directionimplicit quations, each of which has a tridiagonal matrix and can be easily solved. The model reproduces a major Kuroshio intrusion north of Luzon Island, one north of Taiwan Island, andone west of the Tokara Strait. The model shows a current system running from the Luzon Strait to the coastof Vietnam and Hainan Island, through the Taiwan Strait and then into the Tsushima Strait. The summerand winter

Numerical simulations and comparative analysis for two types of storm surges in the Bohai Sea using a coupled atmosphere-ocean model

The Bohai Sea is extremely susceptible to storm surges induced by extratropical storms and tropical cyclones in nearly every season. In order to relieve the impacts of storm surge disasters on structures and human lives in coastal regions, it is very important to understand the occurring of the severe storm surges. The previous research is mostly restricted to a single type of storm surge caused by extratropical storm or tropical cyclone. In present paper, a coupled atmosphere-ocean model is developed to study the storm surges induced by two types of extreme weather conditions. Two special cases happened in the Bohai Sea are simulated successively. The wind intensity and minimum sea-level pressure derived from the Weather Research and Forecasting (WRF) model agree well with the observed data. The computed time series of water level obtained from the Regional Ocean Modeling System (ROMS) also are in good agreement with the tide gauge observations. The structures of the wind fields and average currents for two types of storm surges are analyzed and compared. The results of coupled model are compared with those from the uncoupled model. The case studies indicate that the wind field and structure of the ocean surface current have great differences between extratropical storm surge and typhoon storm surge. The magnitude of storm surge in the Bohai Sea is shown mainly determined by the ocean surface driving force, but greatly affected by the coastal geometry and bathymetry.

An ocean-land-atmosphere coupled model for tropical cyclone landfall processes: The multi-layer ocean model and its verification

POM (Princeton ocean model) tentatively taken as the ocean part of an ocean-land atmosphere coupled modcl is verified for the ultimate purpose of studying the landfall process of tropical cyclone (TC) in the western North Pacific. The POM is tested with monthly mean wind stress in the summer and given lateral boundary conditions. The results indicate that the equilibrium state of the ocean is in accordance with the climate mean, with the error in sea surface temperature (salinity) less than 0.5 ℃ (0.5). The simulated occan currents are reasonable as well.Several numerical experiments are designed to verify the oceanic response to a stationary or moving TC. It is found that the results agree fairly well with the previous work, including both the drop magnitude and the distribution ofsca temperature. Compared with the simple two-layer ocean model used by some other studies, the response of the ocean to a TC is more logical here. The model is also verified in a real case with a TC passing the neighborhood of a buoy station. It is shown that the established ocean model can basically reproduce the sea surface temperature change as observed.

A harmonic analyzed parameterization of tide-induced mixing for ocean models

The tide-induced mixing plays an important role in the regulation of ocean circulation.Numerical simulation of continental shelf circulation is found to exhibit an unreasonable vertical thermohaline structure without consideration of tide effects.In this study,we establish a harmonic analyzed parameterization of tide-induced（HAT） mixing,by which means to derive time-depended function of mixing coefficient based on harmonic analysis of the vertical mixing coefficient.By employing HAT mixing parameterization scheme,a series of numerical experiments are conducted for the Yellow Sea.Numerical results show that an ocean circulation model with the HAT mixing involved is capable of reproducing the reasonable thermohaline structure of the Yellow Sea Cold Water Mass,similar to structures produced by explicit tidal forcing on the open boundary.The advantage of the HAT method is its faster computation time,compared with models that directly resolve explicit tidal motion.The HAT parameterization for the tide-induced mixing has potential to improve both the accuracy and efficiency of ocean circulation and climate models.

Water exchange through the Kerama Gap estimated with a 25-year Pacific HYbrid Coordinate Ocean Model

Variations in water exchange through the Kerama Gap(between Okinawa Island and Miyakojima Island) from 1979 to 2003 were estimated with the 0.08° Pacific HYbrid Coordinate Ocean Model(HYCOM). The model results show that the mean transport through the Kerama Gap(KGT) from the Pacific Ocean to the East China Sea(ECS) was 2.1 Sv, which agrees well with the observed mean KGT(2.0 Sv) for 2009–2010. Over the time period examined, the monthly KGT varied from-10.9 Sv to 15.8 Sv and had a standard deviation of ? 5.0 Sv. The water mainly enters the ECS via the subsurface layer(300–500 m) along the northeastern slope of the Kerama Gap and mainly flows out of the ECS into the southwest of the Kerama Gap. The seasonal and interannual variations of the KGT and the Kuroshio upstream transport were negatively correlated. The Kuroshio upstream transport was largest in summer and smallest in autumn while the KGT was smallest in summer(1.02 Sv) and largest in spring(2.94 Sv) and autumn(2.44 Sv). The seasonal and interannual variations in the Kuroshio downstream(across the PN-line) transport dif fered significantly from the Kuroshio upstream transport but corresponded well with the KGT and the sum of the transport through the Kerama Gap and the Kuroshio upstream, which indicates that information about variation in the KGT is important for determining variation in the Kuroshio transport along the PN-line.

A Numerical Study of a TOGA-COARE Squall-Line Using a Coupled Mesoscale Atmosphere-Ocean Model

An atmosphere-ocean coupled mesoscale modeling system is developed and used to investigate the interactions between a squall line and the upper ocean observed over the western Paci?c warm pool during the Tropical Ocean/Global Atmosphere Coupled Ocean and Atmosphere Response Experiment (TOGA-COARE). The modeling system is developed by coupling the Advanced Regional Prediction Sys- tem (ARPS) to the Princeton Ocean Model (POM) through precipitation and two-way exchanges of mo- mentum, heat, and moisture across the air-sea interface. The results indicate that the interaction between the squall-line and the upper ocean produced noticeable di?erences in the sensible and latent heat ?uxes, as compared to the uncoupled cases. Precipitation, which is often ignored in air-sea heat ?ux estimates, played a major role in the coupling between the mesoscale convective system and the ocean. Precipitation a?ected the air-sea interaction through both freshwater ?ux and sensible heat ?ux. The former led to the formation of a thin stable ocean layer underneath and behind the precipitating atmospheric convection. The presence of this stable layer resulted in a more signi?cant convection-induced sea surface temperature (SST) change in and behind the precipitation zone. However, convection-induced SST changes do not seem to play an important role in the intsensi?cation of the existing convective system that resulted in the SST change, as the convection quickly moved away from the region of original SST response.

A coupled regional Arctic sea ice-ocean model:configuration and application

A regional sea ice-ocean coupled model for the Arctic Ocean was developed, based on the MlTgcm ocean circulation model and classical Hibler79 type two category thermodynamics-dynamics sea ice model.The sea ice dynamics and thermodynamics were considered based on Viscous-Plastic（VP） and Winton three-layer models,respectively.A detailed configuration of coupled model has been introduced.Special attention has been paid to the model grid setup,subgrid paramerization,ice-ocean coupling and open boundary treatment.The coupled model was then applied and two test run examples were presented.The first model run was a climatology simulation with 10 years（1992—2002） averaged NCAR/NCEP reanalysis data as atmospheric forcing.The second model run was a seasonal simulation for the period of 1992—2007.The atmospheric forcing was daily NCAR/NCEP reanalysis.The climatology simulation captured the general pattern of the sea ice thickness distribution of the Arctic,i.e.,the thickest sea ice is situated around the Canada Archipelago and the north coast of the Greenland. For the second model run,the modeled September Sea ice extent anomaly from 1992—2007 was highly correlated with the observations,with a linear correlation coefficient of 0.88.The minimum of the Arctic sea ice area in the September of 2007 was unprecedented.The modeled sea ice area and extent for this minimum was overestimated relative to the observations.However,it captured the general pattern of the sea ice retreat.

A coupled ice-ocean model for the Bohai SeaⅡ.Case study

The coupled ice- ocean model for the Bohai Sea is used for simulating the freezing, melting, and variation of ice cover and the heat bal- ance at the sea- ice, air- ice, and air- sea interfaces of the Bohai Sea during the entire winter in 1998￣1999 and 2000￣2001. The cou- pled model is forced by real time numerical weather prediction fields. The results show that the thermodynamic effects of atmosphere and ocean are very important for the evolvement of ice in the Bohai Sea, especially in the period of ice freezing and melting. Ocean heat flux plays a key role in the thermodynamic coupling. The simulation also presents the different thermodynamic features in the ice covered region and the marginal ice zone. Ice thickness, heat budget at the interface, and surface sea temperature, etc. between the two representative points are discussed.