Guijun Han

Application of the Multigrid Data Assimilation Scheme to the China Seas’ Temperature Forecast

Correlation scales have been used in the traditional scheme of three-dimensional variational data assimilation (3DVAR) to estimate the background (or first guess) error covariance matrix (the B matrix in brief) for the numerical forecast and reanalysis of ocean for decades. However, it is challenging to implement this scheme. On the one hand, determining the correlation scales accurately can be difficult. On the other hand, the positive definite of the B matrix cannot be guaranteed unless the correlation scales are sufficiently small. Xie et al. indicated that a traditional 3DVAR only corrects certain wavelength errors, and its accuracy depends on the accuracy of the B matrix. Generally speaking, the shortwave error cannot be sufficiently corrected until the longwave error is corrected. An inaccurate B matrix may mistake longwave errors as shortwave ones, resulting in erroneous analyses. A new 3DVAR data assimilation scheme, called a multigrid data assimilation scheme, is proposed in this paper for quickly minimizing longwave and shortwave errors successively. By assimilating the sea surface temperature and temperature profile observations into a numerical model of the China Seas, this scheme is applied to a retroactive real-time forecast experiment and favorable results are obtained. Compared to the traditional scheme of 3DVAR, this new scheme has higher forecast accuracy and lower root-meansquare errors. Note that the new scheme demonstrates greatly improved numerical efficiency in the analysis procedure.

Summer surface layer thermal response to surface gravity waves in the Yellow Sea

The Princeton Ocean Model (POM) with generalized coordinate system (POMgcs) is used to study the summer surface-layer thermal response to surface gravity waves in the Yellow Sea (YS). The parameterization schemes of wave breaking developed by Mellor and Blumberg (J Phys Oceanogr 34:693–698, 2004) and Kantha and Clayson (Ocean Model 6:101–124, 2004), respectively, and Stokes production developed by Kantha and Clayson (Ocean Model 6:101–124, 2004) are both included in the Mellor–Yamada turbulence closure model Mellor and Yamada (Rev Geophys 20:851–875, 1982) of POMgcs. Numerical results show that surface gravity waves impact the depth of surface mixed layer of temperature in the YS in summer. The surface mixed layer in the YS cannot be reproduced well and has a visible difference from the observation if the parameterization schemes are not included. A diagnostic analysis of turbulent kinetic energy suggests that both Stokes production and wave breaking play key roles in enhancing the turbulent mixing near the sea surface in the YS. Stokes production seems to have a greater impact throughout the upper mixed layer in the YS in summer than that of wave breaking. In addition, a diagnostic analysis of the momentum balance shows that Coriolis–Stokes forcing has a significant effect on the momentum budget in the upper layer in the YS, and surface gravity waves are able to reduce the velocity of mean flow near the surface and make the mean flow near the surface more homogeneous vertically in the YS.

Improvement of short-term forecasting in the northwest Pacific through assimilating Argo data into initial fields

The impact of assimilating Argo data into an initial field on the short-term forecasting accuracy of temperature and salinity is quantitatively estimated by using a forecasting system of the western North Pacific, on the base of the Princeton ocean model with a generalized coordinate system (POMgcs). This system uses a sequential multigrid three-dimensional variational (3DVAR) analysis scheme to assimilate observation data. Two numerical experiments were conducted with and without Argo temperature and salinity profile data besides conventional temperature and salinity profile data and sea surface height anomaly (SSHa) and sea surface temperature (SST) in the process of assimilating data into the initial fields. The forecast errors are estimated by using independent temperature and salinity profiles during the forecasting period, including the vertical distributions of the horizontally averaged root mean square errors (H-RMSEs) and the horizontal distributions of the vertically averaged mean errors (MEs) and the temporal variation of spatially averaged root mean square errors (S-RMSEs). Comparison between the two experiments shows that the assimilation of Argo data significantly improves the forecast accuracy, with 24% reduction of H-RMSE maximum for the temperature, and the salinity forecasts are improved more obviously, averagely dropping of 50% for H-RMSEs in depth shallower than 300 m. Such improvement is caused by relatively uniform sampling of both temperature and salinity from the Argo drifters in time and space.

A theoretical study of the multigrid three-dimensional variational data assimilation scheme using a simple bilinear interpolation algorithm

In order to solve the so-called “bull-eye” problem caused by using a simple bilinear interpolation as an observational mapping operator in the cost function in the multigrid three-dimensional variational (3DVAR) data assimilation scheme, a smoothing term, equivalent to a penalty term, is introduced into the cost function to serve as a means of troubleshooting. A theoretical analysis is first performed to figure out what on earth results in the issue of “bull-eye”, and then the meaning of such smoothing term is elucidated and the uniqueness of solution of the multigrid 3DVAR with the smoothing term added is discussed through the theoretical deduction for one-dimensional (1D) case, and two idealized data assimilation experiments (one- and two-dimensional (2D) cases). By exploring the relationship between the smoothing term and the recursive filter theoretically and practically, it is revealed why satisfied analysis results can be achieved by using such proposed solution for the issue of the multigrid 3DVAR.

A numerical estimation of the impact of Stokes drift on upper ocean temperature

The impact of Stokes drift on the mixed layer temperature variation was estimated by taking into account an advective heat transport term induced by the Stokes drift in the equation of mixed layer temperature and using the oceanic and wave parameters from a global ocean circulation model (HYCOM) and a wave model (Wave Watch III). The dimensional analysis and quantitative estimation method were conducted to assess the importance of the effect induced by the Stokes drift and to analyze its spatial distribution and seasonal variation characteristics. Results show that the contribution of the Stokes drift to the mixed layer temperature variation at mid-to-high latitudes is comparable with that of the mean current, and a substantial part of mixed layer temperature change is induced by taking the Stokes drift effect into account. Although the advection heat transport induced by the Stokes drift is not the leading term for the mixed layer temperature equation, it cannot be neglected and even becomes critical in some regions for the simulation of the upperocean temperature.

The Reanalysis of Currents and Throughflow Volume Transport in the Taiwan Strait

Currents and throughflow transport in the Taiwan Strait (TS) are systematically analyzed using a new regional ocean reanalysis dataset (China Ocean ReAnalysis, CORA in brief), recently released by National Marine Data and Information Service (NMDIS) in 2009, covering Chinas coastal waters and adjacent seas. The total throughflow volume transport through TS presents a seasonal variability with the annual-mean of 1.49Sv. It reaches maximum in summer with mean value of 2.33Sv, while in winter the mean transport falls to 0.78Sv on average. The estimated throughflow transport is in good agreement with the observational results. The comparisons of the results with historical observations and numerical simulations also show that CORA performs well on reproducing the currents and temperature field in TS region. In addition, CORA not only successfully reproduced the known oceanic phenomena in TS such as temperature front, up-wind flow in fall, winter northward flow blocking in Penghu Channel, and the winter anticyclonic eddy located at north of Changyun Rise (CYR), but it also found a new cyclonic eddy in north TS in winter.

Several Important Issues in Salinity Quality Control of Argo Float

Several floats that deployed by different countries during the last years over the global ocean, which served as a contribution to the Array for real time geostrophic oceanography (Argo), are studied in this paper with respect to their uniqueness in the Delayed Mode Quality Control (DMQC) process. Of these floats, issues were found in the temperature-salinity relationships during the DMQC, which would easily be misinterpreted or overlooked. This study lists and summarizes the possible situations of the floats that might confront in their observation life, together with the reasonable explanations. The authors present some thoughts on evaluating the performance of the conductivity sensors and calibrating the float salinity data.

Impact of sea spray on upper ocean temperature during typhoon passage: simulation with a 1-D turbulent model

At the interface between the lower atmosphere and sea surface, sea spray might significantly influence air-sea heat fluxes and subsequently, modulate upper ocean temperature during a typhoon passage. The effects of sea spray were introduced into the parameterization of sea surface roughness in a 1-D turbulent model, to investigate the effects of sea spray on upper ocean temperature in the Kuroshio Extension area, for the cases of two real typhoons from 2006, Yagi and Soulik. Model output was compared with data from the Kuroshio Extension Observatory (KEO), and Reynolds and AMSRE satellite remote sensing sea surface temperatures. The results indicate drag coefficients that include the spray effect are closer to observations than those without, and that sea spray can enhance the heat fluxes (especially latent heat flux) considerably during a typhoon passage. Consequently, the model results with heat fluxes enhanced by sea spray simulate better the cooling process of the SST and upper-layer temperature profiles. Additionally, results from the simulation of the passage of typhoon Soulik (that passed KEO quickly), which included the sea spray effect, were better than for the simulated passage of typhoon Yagi (that crossed KEO slowly). These promising 1-D results could provide insight into the application of sea spray in general circulation models for typhoon studies.

A Study of Transport and Impact Strength of Fukushima Nuclear Pollutants in the North Pacific Surface

Based on the statistics of surface drifter data of 1979–2011 and the simulation of nuclear pollutant particulate movements simulated using high quality ocean reanalysis surface current dataset, the transport pathways and impact strength of Fukushima nuclear pollutants in the North Pacific have been estimated. The particulates are used to increase the sampling size and enhance the representativeness of statistical results. The trajectories of the drifters and particulates are first examined to identify typical drifting pathways. The results show that there are three types of transport paths for nuclear pollutants at the surface: 1) most pollutant particles move eastward and are carried by the Kuroshio and Kuroshio-extension currents and reach the east side of the North Pacific after about 3.2–3.9 years; 2) some particles travel with the subtropical circulation branch and reach the east coast of China after about 1.6 years according to one drifter trajectory and about 3.6 years according to particulate trajectories; 3) a little of them travel with local, small scale circulations and reach the east coast of China after about 1.3–1.8 years. Based on the particulates, the impact strength of nuclear pollutants at these time scales can be estimated according to the temporal variations of relative concentration combined with the radioactive decay rate. For example, Cesium-137, carried by the strong North Pacific current, mainly accumulates in the eastern North Pacific and its impact strength is 4% of the initial level at the originating Fukushima area after 4 years. Due to local eddies, Cesium-137 in the western North Pacific is 1% of the initial pollutant level after 1.5 years and continuously increases to 3% after 4 years. The vertical movement of radioactive pollutants is not taken into account in the present study, and the estimation accuracy would be improved by considering three-dimensional flows.

Tropical cyclone footprint in the ocean mixed layer observed by Argo in the Northwest Pacific

This study systematically investigated the ocean mixed layer responses to tropical cyclone (TC) using available Argo profiles during the period of 1998–2011 in the northwest Pacific. Results reveal that isothermal layer (IL) deepening and isothermal layer (IL) cooling with evident rightward biases induced by strong TCs are clearer compared to the weak TCs. Likewise, the rightward biases of IL deepening and cooling induced by fast TCs are more obvious than that induced by slow TCs. The upwelling within TCs eye is much stronger for the strong (slow) TCs than weak (fast) TCs. For the strong and slow TCs, the TC-induced rainfall reduces deepening of constant density layer (with its depth called the mixed layer depth, MLD), and in turn increases the barrier layer thickness (BLT). The initial BL prior to TC can restrict IL cooling more markedly under the weak and fast TCs than under the strong and slow TCs. The inertial oscillation is stronger induced by the strong (fast) TCs than by the weak (slow) TCs. In addition, the most pronounced TC-induced mixed layer deepening and IL cooling in July to October climatology occur in the subtropical gyre of the northwest Pacific with enhanced vertical diffusivity. The maximum increase of isothermal layer depth (ILD) and MLD is up to 5 m, with IL cooling up to 0.4°C.